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Preface 



The development of a consistent picture of the processes of decoherence and 
quantum measurement is among the most interesting fundamental problems 
with far-reaching consequences for our understanding of the physical world. 
A satisfactory solution of this problem requires a treatment which is com- 
patible with the theory of relativity, and many diverse approaches to solve or 
circumvent the arising difficulties have been suggested. This volume collects 
the contributions of a workshop on Relativistic Quantum Measurement and 
Decoherence held at the Istituto Italiano per gli Studi Filosofici in Naples, 
April 9-10, 1999. The workshop was intended to continue a previous meeting 
entitled Open Systems and Measurement in Relativistic Quantum Theory, the 
talks of which are also published in the Lecture Notes in Physics Series (Vol. 
526). 

The different attitudes and concepts used to approach the decoherence 
and quantum measurement problem led to lively discussions during the work- 
shop and are reflected in the diversity of the contributions. In the first article 
the measurement problem is introduced and the various levels of compatibility 
with special relativity are critically reviewed. In other contributions the roles 
of non-locality and entanglement in quantum measurement and state vector 
preparation are discussed from a pragmatic quantum-optical and quantum- 
information perspective. In a further article the viewpoint of the consistent 
histories approach is presented and a new criterion is proposed which refines 
the notion of consistency. Also, the phenomenon of decoherence is examined 
from an open system’s point of view and on the basis of superselection rules 
employing group theoretic and algebraic methods. The notions of hard and 
soft superselection rules are addressed, as well as the distinction between real 
and apparent loss of quantum coherence. Furthermore, the emergence of real 
decoherence in quantum electrodynamics is studied through an investigation 
of the reduced dynamics of the matter variables and is traced back to the 
emission of bremsstrahlung. 

It is a pleasure to thank Avv. Gerardo Marotta, the President of the Is- 
tituto Italiano per gli Studi Filosofici, for suggesting and making possible 
an interesting workshop in the fascinating environment of Palazzo Serra di 
Cassano. Furthermore, we would like to express our gratidude to Prof. An- 
tonio Gargano, the General Secretaty of the Istituto Italiano per gli Studi 
Filosofici, for his friendly and efficient local organization. We would also like 
to thank the participants of the workshop. 
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Special Relativity as an Open Question 



David Z. Albert 

Department of Philosophy, Columbia University, New York, USA 



Abstract. There seems to me to be a way of reading some of the trouble we have 
lately been having with the quantum-mechanical measurement problem (not the 
standard way, mind you, and certainly not the only way; but a way that nonethe- 
less be worth exploring) that suggests that there are fairly prosaic physical circum- 
stances under which it might not be entirely beside the point to look around for 
observable violations of the special theory of relativity. The suggestion I have in mind 
is connected with attempts over the past several years to write down a relativistic 
field-theoretic version of the dynamical reduction theory of Ghirardi, Rimini, and 
Weber [Physical Review D34, 470-491 (1986)], or rather it is connected with the 
persistent failure of those attempts, it is connected with the most obvious strategy 
for giving those attempts up. And that (in the end) is what this paper is going to 
be about. 



1 The Measurement Problem 

Let me start out (however) by reminding you of precisely what the quantum- 
mechanical problem of measurement is, and then talk a bit about where 
things stand at present vis-a-vis the general question of the compatibility of 
quantum mechanics with the special theory of relativity, and then I want to 
present the simple, standard, well-understood non-relativistic version of the 
Ghirardi, Rimini, and Weber (GRW) theory [1], and then (at last) I will get 
into the business I referred to above. 

First the measurement problem. Suppose that every system in the world 
invariably evolves in accordance with the linear deterministic quantum-mecha- 
nical equations of motion and suppose that M is a good measuring instrument 
for a certain observable A of a certain physical system S. What it means for 
M to be a “good” measuring instrument for A is just that for all eigenvalues 
tti of A: 



|ready)M|A = ai)s — > jindicates that A = ai)M\A = ai)s, (1) 

where |ready)M is that state of the measuring instrument M in which M 
is prepared to carry out a measurement of A, “ — denotes the evolution 
of the state of M -I- S' during the measurement-interaction between those 
two systems, and [indicates that A = Qi)M is that state of the measuring 
instrument in which, say, its pointer is pointing to the the Oj-position on its 
dial. That is: what it means for M to be a “good” measuring instrument 
for A is just that M invariably indicates the correct value for A in all those 
states of S in which A has any definite value. 



H.-P. Breuer and F. Petruccione (Eds.): Proceedings 1999, LNP 559, pp. 1-13, 2000. 
(c) Springer- Verlag Berlin Heidelberg 2000 
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The problem is that (1), together with the linearity of the equations of 
motion entails that: 

^|ready)M|^ = az)s^^ jindicates that A = ai)MlA = a^}s- (2) 

i i 

And that appears not to be what actually happens in the world. The right- 
hand side of Eq. (2) is (after all) a superposition of various different outcomes 
of the A-measurement - and decidedly not any particular one of them. But 
what actually happens when we measure A on a system S' in a state like 
the one on the left-hand-side of (2) is of course that one or another of those 
particular outcomes, and nothing else, emerges. 

And there are two big ideas about what to do about that problem that 
seem to me to have any chance at all of being on the right track. 

One is to deny that the standard way of thinking about what it means to 
be in a superposition is (as a matter of fact) the right way of thinking about 
it; to deny, for example, that there fails to be any determinate matter of fact, 
when a quantum state like the one here obtains, about where the pointer is 
pointing. 

The idea (to come at it from a slightly different angle) is to construe 
quantum-mechanical wave-functions as less than complete descriptions of the 
world. The idea that something extra needs to be added to the wave-function 
description, something that can broadly be thought of as choosing between 
the two conditions superposed here, something that can be thought of as 
somehow marking one of those two conditions as the unique, actual, outcome 
of the measurement that leads up to it. 

Bohm’s theory is a version of this idea, and so are the various modal 
interpretations of quantum mechanics, and so (more or less) are many-minds 
interpretations of quantum mechanics.^ 

The other idea is to stick with the standard way of thinking about what 
it means to be in a superposition, and to stick with the idea that a quantum- 
mechanical wave-function amounts, all by itself, to a complete description of 
a physical system, and to account for the emergence of determinate outcomes 
of experiments like the one we were talking about before by means of explicit 
violations of the linear deterministic equations of motion, and to try to de- 
velop some precise idea of the circumstance s under which those violations 
occur. 

And there is an enormously long and mostly pointless history of specula- 
tions in the physical literature (speculations which have notoriously hinged on 
distinctions between the “microscopic” and the “macroscopic”, or between 

^ Many-minds interpretations are a bit of a special case, however. The outcomes 
of experiments on those interpretations (although they are perfectly actual) are 
not unique. The more important point, though, is that those interpretations (like 
the others I have just mentioned) solve the measurement problem by construing 
wave- functions as incomplete descriptions of the world. 
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the “reversible” and the “irreversible”, or between the “animate” and the 
“inanimate”, or between “subject” and “object”, or between what does and 
what doesn’t genuinely amount to a “measurement”) about precisely what 
sorts of violations of those equations - what sorts of collapses - are called for 
here; but there has been to date only one fully- worked-out, traditionally sci- 
entific sort of proposal along these lines, which is the one I mentioned at the 
beginning of this paper, the one which was originally discovered by Ghirardi 
and Rimini and Weber, and which has been developed somewhat further by 
Philip Pearle and John Bell. 

There are (of course) other traditions of thinking about the measurement 
problem too. There is the so-called Copenhagen interpretation of quantum 
mechanics, which I shall simply leave aside here, as it does not even pretend 
to amount to a realistic description of the world. And there is the tradition 
that comes from the work of the late Hugh Everett, the so called “many 
worlds” tradition, which is (at first) a thrilling attempt to have one’s cake 
and eat it too, and which (more particularly) is committed both to the propo- 
sition that quantum-mechanical wave-functions are complete descriptions of 
physical systems and to the proposition that those wave-functions invariably 
evolve in accord with the standard linear quantum- mechanical equations of 
motion, and which (alas, for a whole bunch of reasons) seems to me not to 
be a particular candidate either.^ 

And that’s about it. 



2 Degrees of Compatibility with Special Relativity 

Now, the story of the compatibility of these attempts at solving the mea- 
surement problem with the special theory of relativity turns out to be unex- 
pectedly rich. It turns out (more particularly) that compatibility with special 
relativity is the sort of thing that admits of degrees. We will need (as a matter 
of fact) to think about five of them - not (mind you) because only five are 
logically imaginable, but because one or another of those five corresponds 
to every one of the fundamental physical theories that anybody has thus far 
taken seriously. 

Let’s start at the top. 

What it is for a theory to be metaphysically compatible with special rel- 
ativity (which is to say: what it is for a theory to be compatible with special 
relativity in the highest degree) is for it to depict the world as unfolding in 
a four-dimensional Minkowskian space-time. And what it means to speak of 
the world as unfolding within a four-dimensional Minkovskian space-time is 
(i) that everything there is to say about the world can straightforwardly be 

^ Foremost among these reasons is that the many-worlds interpretations seems to 
me not to be able to account for the facts about chance. But that’s a long story, 
and one that’s been told often enough elsewhere. 
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read off of a catalogue of the local physical properties at every one of the con- 
tinuous infinity of positions in a space-time like that, and (ii) that whatever 
lawlike relations there may be between the values of those local properties can 
be written down entirely in the language of a space-time that - that whatever 
lawlike relations there may be between the values of those local properties 
are invariant under Lorentz- transformations. And what it is to pick out some 
particular inertial frame ofreference in the context of the sort of theory we’re 
talking about here - what it is (that is) to adopt the conventions of measure- 
ment that are indigenous to any particular frame of reference in the context 
of the sort of theory we’re talking about here - is just to pick out some par- 
ticular way of organising everything there is to say about the world into a 
story, into a narrative, into a temporal sequen ce of instantaneous global phys- 
ical situations. And every possible world on such a theory will invariably be 
organizable into an infinity of such stories - and those stories will invariably 
be related to one another by Lorentz-transformations. And note that if even 
a single one of those stories is in accord with the laws, then (since the laws 
are invariant under Lorentz-transformations) all of them must be. 

The Lorentz-invariant theories of classical physics (the electrodynamics 
of Maxwell, for example) are metaphysically compatible with special relativ- 
ity; and so (more surprisingly) are a number of radically non-local theories 
(completely hypothetical ones, mind you - ones which in so far as we know 
at present have no application whatever to the actual world) which have 
recently appeared in the literature.^ 

But it happens that not a single one of the existing proposals for making 
sense of quantum mechanics is metaphysically compatible with special rela- 
tivity, and (moreover) it isn’t easy to imagine there ever being one that is. The 
reason is simple: What is absolutely of the essence of the quantum-mechanical 
picture of the world (in so far as we understand it at present), what none 
of the attempts to straighten quantum mechanics out have yet dreamed of 
dispensing with, are wave-functions. And wave-functions just don’t live in 
four-dimensional space-times; wave-functions (that is) are just not the sort 
of objects which can always be uniquely picked out by means of any cata- 
logue of the local properties of the positions of a space-time like that. As 
a general matter, they need bigger ones, which is to say higher- dimensional 
ones, which is to say configurational ones. And that (alas!) is that. 

The next level down (let’s call this one the level of dynamical compati- 
bility with special relativity) is inhabited by pictures on which the physics 
of the world is exhaustively described by something along the lines of a (so- 
called) relativistic quantum field theory - a pure one (mind you) in which 

® Tim Maudlin and Frank Artzenius have both been particularly ingenious in con- 
cocting theories like these, which (notwithstanding their non-locality) are entirely 
formulable in four-dimensional Minkowski space-time. Maudlin’s book Quantum 
Non-Locality and Relativity (Blackwell, 1994) contains extremely elegant discus- 
sions of several such theories. 
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there are no additional variables, and in which the quantum states of the 
world invariably evolve in accord with local, deterministic, Lorentz-invariant 
quantum mechanical equations of motion. These pictures (once again) must 
depict the world as unfolding not in a Minkowskian space-time but in a con- 
figuration one - and the dimensionalities of the configuration space-times in 
question here are (of course) going to be infinite. Other than that, however, 
everything remains more or less as it was above. The configuration space-time 
in question here is built directly out of the Minkowskian one (remember) by 
treating each of the points in Minkowskian space-time (just as one does in 
the classical theory of fields) as an instantaneous bundle of physical degrees 
of freedom. And so what it is to pick out some particular inertial frame of 
reference in the context of this sort of picture is still just to pick out some 
particular way of organizing everything there is to say about the world into 
a temporal sequence of instantaneous global physical situations. And every 
possible world on this sort of a theory will still be organizable into an infinity 
of such stories. And those stories will still be related to one another by means 
of the appropriate generalizations of the Lorentz point-transformations. And 
it will still be the case that if even a single one of those stories is in accord with 
the laws, then (since the laws are invariant under Lorentz-transformations) 
all of them must be. 

The trouble is that there may well not he any such pictures that turn 
out to be worth taking seriously. All we have along these lines at present 
(remember) are the many-worlds pictures (which I fear will turn out not to 
be coherent) and the vaany-minds pictures (which I fear will turn out not to 
be plausible). 

And further down things start to get ugly. 

We have known for more than thirty years now that any proposal for 
making sense of quantum mechanics on which measurements invariably have 
unique and particular and determinate outcomes (which covers all of the 
proposals I know about, or at any rate the ones I know about that are also 
worth thinking about, other than many worlds and many minds) is going to 
have no choice whatever but to turn out to be non-local. 

Now, non-locality is certainly not an obstacle in and of itself even to meta- 
physical compatibility with special relativity. There are now (as I mentioned 
before) a number of explicit examples in the literature of hypothetical dy- 
namical laws which are radically non-local and which are nonetheless cleverly 
cooked up in such a way as to be formulable entirely within Minkowski-space. 
The thing is that none of them can even remotely mimic the empirical pre- 
dictions of quantum mechanics; and that nobody I talk to thinks that we 
have even the slightest reason to hope for one that will. 

What we do have (on the other hand) is a very straightforward trick by 
means of which a wide variety of theories which are radically non-local and 
(moreover) are flatly incompatible with the proposition that the stage on 
which physical history unfolds is Minkowki-space can nonetheless be made 
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fully and trivially Lorentz-invariant; a trick (that is) by means of which a wide 
variety of such theories can be made what you might call formally compatible 
with special relativity. 

The trick [2] is just to let go of the requirement that the physical history 
of the world can be represented in its entirety as a temporal sequence of 
situations. The trick (more particularly) is to let go of the requirement that 
the situation associated with two intersecting space-like hypersurfaces in the 
Minkowski-space must agree with one another about the expectations values 
of local observables at points where the two surfaces coincide. 

Consider (for example) an old-fashioned non-relativistic projection-postu- 
late, which stipulates that the quantum states of physical systems invariably 
evolve in accord with the linear deterministic equations of motion except when 
the system in question is being “measured” ; and that the quantum state of a 
system instantaneously jumps, at the instant the system is measured, into the 
eigenstate of the measured observable corresponding to the outcome of the 
measurement. This is the sort of theory that (as I mentioned above) nobody 
takes seriously anymore, but never mind that; it will serve us well enough, for 
the moment, as an illustration. Here’s how to make this sort of a projection- 
postulate Lorentz-invariant: First, take the linear collapse-free dynamics of 
the measured system - the dynamics which we are generally in the habit of 
writing down as a deterministic connection between the wave-functions on 
two arbitrary equal-time-hyperplanes - and re-write it as a deterministic con- 
nection between the wave-functions on two arbitrary space-like-hypersurfaces, 
as in Fig. 1. Then stipulate that the jumps referred to above occur not (as it 
were) when the equal-time-hyperplane sweeps across the measurement-event, 
but whenever an arbitrary space-like hypersurface undulates across it. 

Suppose (say) that the momentum of a free particle is measured along 
the hypersurface marked t = 0 in Fig. 2, and that later on a measurement 
locates the particle at P. Then our new projection-postulate will stipulate 
among other things) that the wave-function of the particle along hypersurface 
a is an eigenstate of momentum, and that the wave-function of the particle 
along hypersurface b is (very nearly) an eigenstate of position. And none of 
that (and nothing else that this new postulate will have to say) refers in any 
way shape or form to any particular Lorentz frame. And this is pretty. 

But think for a minute about what’s been paid for it. As things stand 
now we have let go not only of Minkowski-space as a realistic description of 
the stage on which the story of the world is enacted, but (in so far as I can 
see) of any conception of that stage whatever. As things stand now (that 
is) we have let go of the idea of the world’s having anything along the lines 
of a narratable story at alll And all this just so as to guarantee that the 
fundamental laws remain exactly invariant under a certain hollowed-out set 
of mathematical transformations, a set which is now of no particular deep 
conceptual interest, a set which is now utterly disconnected from any idea of 
an arena in which the world occurs. 




8 



David Z Albert 



Never mind. Suppose we had somehow managed to resign ourselves to 
that. There would still be trouble. It happens (you see) that notwithstanding 
the enormous energy and technical ingenuity has been expended over the past 
several years in attempting to concoct a version of a more believable theory 
of collapses - a version (say) of GRW theory - on which a trick like this might 
work, even that (paltry as it is) is as yet beyond our grasp. 

And all that (it seems to me) ought to give us pause. 

The next level down (let’s call this one the level of discreet incompatibility 
with special relativity) is inhabited by theories (Bohm’s theory, say, or modal 
theories) on which the special theory of relativity, whatever it means, is unam- 
biguously false; theories (that is) which explicitly violate Lorentz-invariance, 
but which nonetheless manage to refrain from violating it in any of their pre- 
dictions about the outcomes of experiments. These theories (to put it slightly 
differently) all require that there be some legally privileged Lorentz-frame, 
but they all also entail that (as a matter of fundamental principle) no per- 
formable experiment can identify what frame that is. 

And then (at last) there are theories that explicitly violate Lorentz- 
invariance (but presumably only a bit, or only in places we haven’t looked 
yet) even in their observable predictions. It’s that sort of a theory (the sort 
of a theory we’ll refer to as manifestly incompatible with special relativity) 
that I’m going to want to draw your attention to here. 

3 The Theory I Have in Mind 

But one more thing needs doing before we get to that, which is to say some- 
thing about where the theory I have in mind comes from. And where it 
comes from (as I mentioned at the outset) is the non-relativistic spontaneous 
localization theory of Ghirardi, Rimini, Weber, and Pearle. 

GRWP’s idea was that the wave function of an N-particle system 

if{ri,r2,...,rN,t) (3) 

usually evolves in the familiar way - in accordance with the Schrodinger equa- 
tion - but that every now and then (once in something like 10^^/iV seconds), 
at random, but with fixed probability per unit time, the wave function is 
suddenly multiplied by a normalized Gaussian (and the product of those two 
separately normalized functions is multiplied, at that same instant, by an 
overall renormalizing constant). The form of the multiplying Gaussian is: 

Kexp[-{r -rkf/2a^] (4) 

where is chosen at random from the arguments r„, and the width a of 
the Gaussian is of the order of 10“® cm. The probability of this Gaussian 
being centered at any particular point r is stipulated to be proportional to 
the absolute square of the inner product of (3) (evaluated at the instant just 




Special Relativity as an Open Question 



9 



prior to this “jump”) with (4). Then, until the next such “jump”, everything 
proceeds as before, in accordance with the Schrodinger equation. The proba- 
bility of such jumps per particle per second (which is taken to be something 
like 10“^®, as I mentioned above), and the width of the multiplying Gaussians 
(which is taken to be something like 10“^ cm) are new constants of nature. 

That’s the whole theory. No attempt is made to explain the occurrence 
of these “jumps”; that such jumps occur, and occur in precisely the way 
stipulated above, can be thought of as a new fundamental law; a beautiful and 
absolutely explicit law of collapse, wherein there is no talk at a fundamental 
level of “measurements” or “recordings” or “macroscopicness” or anything 
like that. 

Note that for isolated microscopic systems (i.e. systems consisting of small 
numbers of particles) “jumps” will be rare as to be completely unobservable 
in practice; and the width of the multiplying Gaussian has been chosen large 
enough so that the violations of conservation of energy which those jumps 
will necessarily produce will be very small (over reasonable time-intervals), 
even for macroscopic systems. 

Moreover, if it’s the case that every measuring instrument worthy of the 
name has got to include some kind of a pointer, which indicates the outcome 
of the measurement, and if that pointer has got to be a macroscopic physical 
object, and if that pointer has got to assume macroscopically different spatial 
positions in order to indicate different such outcomes (and all of this seems 
plausible enough, at least at first), then the GRW theory can apparently 
guarantee that all measurements have outcomes. Here’s how: Suppose that 
the GRW theory is true. Then, for measuring instruments (M) such as were 
just described, superpositions like 

|H)|M indicates that A) + \B)\M indicates that B) (5) 

(which will invariably be superpositions of macroscopically different localized 
states of some macroscopic physical object) are just the sorts of superposi- 
tions that don’t last long. In a very short time, in only as long as it takes 
for the pointer’s wave-function to get multiplied by one of the GRW Gaus- 
sian (which will be something of the order of 10^®/fV seconds, where N is 
the number of elementary particles in the pointer) one of the terms in (5) 
will disappear, and only the other will propagate. Moreover, the probability 
that one term rather than another survives is (just as standard Quantum 
Mechanics dictates) proportional to the fraction of the norm which it carries. 

And maybe it’s worth mentioning here that there are two reasons why this 
particular way of making experiments have outcomes strikes me at present 
as conspicuously more interesting than others I know about. 

The first has to do with questions of ontological parsimony: We have no 
way whatever of making experiments have outcomes (after all) that does 
without wave-functions. And only many-worlds theories and collapse theo- 
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ries manage to do without anything other than wave- functions^ . And many- 
worlds theories don’t appear to work. 

The second (which strikes me as more important) is that the GRW theory 
affords a means of reducing the probabilities of Statistical Mechanics entirely 
to the probabilities of Quantum Mechanics. It affords a means (that is) of re- 
arranging the foundations of the entirety of physics so as to contain exactly 
one species of chance. And no other way we presently have of making mea- 
surements have outcomes - not Bohm’s theory and not modal theories and 
not many-minds theories and not many-worlds theories and not the Copen- 
hagen interpretation and not quantum logic and not even the other collapse 
theories presently on the market - can do anything like that.^ But let me go 
back to my story. 

4 Approximate Compatibility with Special Relativity 

The trouble (as we’ve seen) is that there can probably not be a version of a 
theory like this which has any sorts of compatibility with special relativity 
that seem worth wanting. 

And the question is what to do about that. 

And one of the things it seems to me one might do is to begin to wonder 
exactly what the all the fuss has been about. One of the things it seems to me 
one might do - given that the theory of relativity is already off the table here 
as a realistic description of the structure of the world - is to begin to wonder 
exactly what the point is of entertaining only those fundamental theories 
which are strictly invariant under Lorentz transformations, or even only those 
fundamental theories whose empirical predictions are strictly invariant under 
Lorentz transformations. 

Why not theories which are are only approximately so? Why not theories 
which violate Lorentz invariance in ways which we would be unlikely to have 
noticed yet? Theories like that, and (more particularly) GRW-like theories 
like that, turn out to be snap to cook up. 

Let’s (finally) think one through. Take (say) standard, Lorentz-invariant, 
relativistic quantum electrodynamics - without a collapse. And add to it some 
non-Lorentz-invariant second-quantized generalization of a collapse-process 
which is designed to reduce - under appropriate circumstances, and in some 
particular preferred frame of reference - to a standard non-relativistic GRW 
Gaussian collapse of the effective wave- function of electrons. And suppose 

^ All other pictures (Bohm, Modal Interpretations, Many-Minds, etc) supplement 
the wave-function with something else; something which we know there to be a 
way of doing without; something which (when you think about it this way) looks 
as if it must somehow be superfluous. 

® This is one of the main topics of a book I have just finished writing, called Time 
and Chance, which is to be published in the fall of 2000 by Harvar d University 
Press. 
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that the frame associated with our laboratory is some frame other than the 
preferred one. And consider what measurements carried out in that labora- 
tory will show. 

This needs to be done with some care. What happens in the lab frame 
is certainly not (for example) that the wave-function gets multiplied by any- 
thing along the lines of a “ Lorentz-transformation” of the non-relativistic 
GRW Gaussian I mentioned a minute ago, for the simple reason that Gaus- 
sians are not the sort of things that are susceptible of having a Lorentz trans- 
formation carried out on them in the first place.® And it is (as a more general 
matter) not to be expected that a theory like this one is going to yield any 
straightforward universal geometrical technique whatever - such as we have 
always had at our disposal, in one form or another, throughout the entire 
modern history of physics - whereby the way the world looks to one observer 
can be read off of the way it looks to some other one, who is in constant 
rectilinear motion relative to the first. The theory we have in front of us at 
the moment is simply not like that. We are (it seems fair to say) in infinitely 
messier waters here. The only absolutely reliable way to proceed on theo- 
ries like this one (unless and until we can argue otherwise) is to deduce how 
things may look to this or that observer by explicitly treating those observers 
and all of their measuring instruments as ordinary physical objects, whose 
states change only and exactly in whatever way it is that they are required 
to change by the microscopic laws of nature, and whose evolutions will pre- 
sumably need to be calculated from the point of view of the unique frame of 
reference in which those laws take on their simplest form. 

That having been said, remember that the violations of Lorentz-invariance 
in this theory arise exclusively in connection with collapses, and that the 
collapses in this theory have been specifically designed so as to have no 
effects whatever, or no effects to speak of, on any of the familiar properties 
or behaviours of everyday localized solid macroscopic objects. And so, in so 
far as we are concerned with things like (say) the length of medium-sized 
wooden dowels, or the rates at which cheap spring-driven wristwatches tick, 
everything is going to proceed, to a very good approximation, as if no such 
violations were occurring at all. 

Let see how far we can run with just that. 

Two very schematic ideas for experiments more or less jump right out at 
you - one of them zeros in on what this theory still has left of the special- 
relativistic length-contraction, and the other on what it still has left of the 
special-relativistic time-dilation. 

The first would go like this: Suppose that the wave-function of a sub- 
atomic particle which is more or less at rest in our lab frame is divided in 
half - suppose, for example, that the wave-function of a neutron whose 2-spin 



The sort of thing you need to start out with, if you want to do a Lorentz trans- 
formation, is not a function of three-space (which is what a Gaussian is) but a 
function of three-space and time. 
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is initially “up” is divided, by means of a Stern-Gerlach magnet, into equal 
y-spin “up” and y-spin “down” components. And suppose that one of those 
halves is placed in box A and that the other half is placed in box B. And 
suppose that those two boxes are fastened on to opposite ends of a little 
wooden dowel. And suppose that they are left in that condition for a certain 
interval - an interval which is to be measured (by the way) in the lab frame, 
and by means of a co-moving cheap mechanical wristwatch. And suppose 
that at the end of that interval the two boxes are brought back together 
and opened, and that we have a look - in the usual way - for the usual 
sort of interference effects. Note (and this is the crucial point here) that the 
length of this dowel, as measured in the preferred frame, will depend radically 
(if the velocity of the lab frame relative to the preferred one is sufficiently 
large) on the dowel’s orientation. If, for example, the dowel is perpendicular 
to the velocity of the lab relative to the preferred frame, it’s length will 
be the same in the preferred frame as in the lab, but if the the dowel is 
parallel to that relative velocity, then it’s length - and hence also the spatial 
separation between A and B - as measured in the preferred frame, will be 
much shorter. And of course the degree to which the GRW collapses wash out 
the interference effects will vary (inversely) with the distance between those 
boxes as measured in the preferred frame. ^ And so it is among the predictions 
of the sort of theory we are entertaining here that if the lab frame is indeed 
moving rapidly with respect to the preferred one, the observed interference 
effects in these sorts of contraptions ought to observably vary as the spatial 
orientation of that device is altered. It is among the consequences of the 
failure of Lorentz-invariance in this theory that (to put it slightly differently) 
in frames other than the preferred one, invariance under spatial rotations 
fails as well. 

The second experiment involves exactly the same contraption, but in this 
case what you do with it is to boost it - particle, dowel, boxes, wristwatch 
and all - in various directions, and to various degrees, but always (so as to 
keep whatever this theory still has in it of the Lorentzian length-contractions 
entirely out of the picture for the moment) perpendicular to the length of the 
dowel. As viewed in the preferred frame, this will yield interference experi- 
ments of different temporal durations, in which different numbers of GRW 
collapses will typically occur, and in which the observed interference effects 
will (in consequence) be washed out to different degrees. 

The sizes of these effects are of course going to depend on things like the 
velocity of the earth relative to the preferred frame (which there can be no 



More particularly: If, in the preferred frame, the separation between the two 
boxes is so small as to be of the order of the width of the GRW Gaussian, the 
washing-out will more or less vanish altogether. 
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way of guessing)®, and the degree to which we are able to boost contraptions 
of the sort I have been describing, and the accuracies with which we are able 
to observe interferences, and so on. 

The size of the effect in the time-dilation experiment is always going to 
vary linearly in where v is the magnitude of whatever boosts we 

find we are able to artificially produce. In the length-contraction experiment, 
on the other hand, the effect will tend to pop in and out a good deal more 
dramatically. If (in that second experiment) the velocity of the contraption 
relative to the preferred frame can somehow be gotten up to the point at 
which is of the order of the width of the GRW Gaussian divided 

by the length of the dowel - either in virtue of the motion of the earth itself, 
or by means of whatever boosts we find we are able to artificially produce, 
or by means of some combination of the two - whatever washing-out there is 
of the interference effects when the length of the dowel is perpendicular to 
its velocity relative to the preferred frame will more or less dis continuously 
vanish when we rotate it. 

Anyway, it seems to me that it might well be worth the trouble to do 
some of the arithmetic I have been alluding to, and to inquire into some of 
our present technical capacities, and to see if any of this might actually be 
worth going out and trying.® 
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® All one can say for certain, I suppose, is that (at the very worst) there must be 
a time in the course of every terrestrial year at which that velocity is at least of 
the order of the velocity of the earth relative to the sun. 

® All of this, of course, leaves aside the question of whether there might be still 
simpler experiments, experiments which might perhaps have already been per- 
formed, on the basis of which the theory we have been talking about here might 
be falsified. It goes without saying that I don’t (as yet) know of any. But that’s 
not saying much. 




Event-Ready Entanglement 



Pieter Kok and Samuel L. Braunstein 
SEECS, University of Wales, Bangor LL57 lUT, UK 



Abstract. We study the creation of polarisation entanglement by means of optical 
entanglement swapping (Zukowski et al, [Phys. Rev. Lett. 71 , 4287 (1993)]). We 
show that this protocol does not allow the creation of maximal ‘event-ready’ en- 
tanglement. Furthermore, we calculate the outgoing state of the swapping protocol 
and stress the fundamental physical difference between states in a Hilbert space and 
in a Fock space. Methods suggested to enhance the entanglement in the outgoing 
state as given by Braunstein and Kimble [Nature 394 , 840 (1998)] generally fail. 



1 Introduction 

Ever since the seminal paper of Einstein, Podolsky and Rosen [1], the concept 
of entanglement has captured the imagination of physicists. The EPR para- 
dox, of which entanglement is the core constituent, points out the non-local 
behaviour of quantum mechanics. This non-locality was quantified by Bell in 
terms of the so-called Bell inequalities [2] and cannot be explained classically. 

Now, at the advent of the quantum information era, entanglement is no 
longer a mere curiosity of a theory which is highly successful in describing 
the natural phenomena. It has become an indispensable resource in quantum 
information protocols such as dense coding, quantum error correction and 
quantum teleportation [3-6]. 

Two quantum systems, parametrised by Xi and X 2 respectively, are called 
entangled when the state ^(xi,X 2 ) describing the total system cannot be 
factorised into states f/’i(a;i) and f/’ 2 (a^ 2 ) of the separate systems: 

'1'{XI,X2) 7^ tpi{xi)tp2{x2) . ( 1 ) 

All the states ^{xi,X 2 ) accessible to two quantum systems form a set S. These 
states are generally entangled. Only in extreme cases 'P{xi,X 2 ) is separable, 
i.e., it can be written as a product of states describing the separate systems. 
The set of separable states form a subset of S with measure zero. 

We arrive at another extremum when the states ^{xi,X 2 ) are maximally 
entangled. The set of maximally entangled states also forms a subset of S 
with measure zero. In order to elaborate on maximal entanglement, we will 
limit our discussion to quantum optics. 

Two photons can be linearly polarised along two orthogonal directions 
X and j/ of a given coordinate system. Every possible state of those two 
photons shared between a pair of modes can be written on the basis of four 
orthonormal states \x, x), [a;, y), \y, x) and \y, y). These basis states generate a 
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four-dimensional Hilbert space. Another possible basis for this space is given 
by the so-called polarisation Bell states: 

\<F^) = (\x,y) ± \y,x))/V2 , 

\<P^) = {\x,x)±\y,y))/V2. (2) 

These states are also orthonormal. They are examples of maximally entangled 
states. The Bell states are not the only maximally entangled states, but they 
are the ones most commonly discussed. For the remainder of this paper we will 
restrict our discussion to the antisymmetric Bell state (in the appendix 
we will explain in more detail why we can do this without loss of generality). 

Suppose we want to conduct an experiment which makes use of polarisa- 
tion entanglement, in particular Ideally, we would like to have a source 

which produces these states at the push of a button. In practice, this might 
be a bit much to ask. A second option is to have a source which only produces 
|>F~) randomly, but flashes a red light when it happens. Such a source would 
create so-called event-ready entanglement: it produces only part of the 
time, but when it does, it tells you so. 

More formally, the outgoing state |V’out|red light flashes) conditioned on the 
red light flashing is said to exhibit event-ready entanglement if it can be 
written as 

|^out|red light flashes ) |iF-) + 0(0 , (3) 

where ^ 1. In the remainder of this paper we will omit the subscript 

‘|red light flashes’ since it is clear that we can only speak of event-ready 
entanglement conditioned on the red light’s flashing. 

Currently, event-ready entanglement has never been produced experimen- 
tally. However, non-maximal entanglement has been created by means of, for 
instance, parametric down-conversion [7]. Rather than a (near) maximally 
entangled state, as in Eq. (3), this process produces states with a large vac- 
uum contribution. Only a minor part consists of an entangled photon state. 
Every time parametric down-conversion is employed, there is only a small 
probability of creating an entangled photon-pair. For the purposes of this 
paper we will call this randomly produced entanglement. 

Parametric down-conversion has been used in several experiments, and 
for some applications randomly produced entanglement therefore seems suf- 
ficient. However, on a theoretical level, maximally entangled states appear as 
primitive notions in many quantum protocols. It is therefore not at all clear 
whether randomly produced entanglement is suitable for all these cases. This 
is one of the main motives in our search for event-ready entanglement, where 
we can ensure that the physical state is maximally entangled. 

In this paper we investigate one particular possibility to create event-ready 
entanglement. It was suggested by Zukowski, Zeilinger, Horne and Ekert [8] 
and Pavicic [9] that entanglement swapping is a suitable candidate. We will 
therefore study this protocol in some detail using quantum optics. 
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Entanglement swapping is essentially the teleportation of one part of an 
entangled pair [3,8,10]. Suppose we have a system of two independent (max- 
imally) entangled photon-pairs in modes a, b and c, d. If we restrict ourselves 
to the Bell states, we have for instance 

\d')abcd = \d^~)ab ® \'I'~)cd ■ (4) 

However, this state can be written on a different basis: 

\d^)abcd = 

+ |<?+)bc . (5) 

If we make a Bell measurement on modes b and c, we can see from Eq. (5) 
that the undetected remaining modes a and d become entangled. For instance, 
when we find modes b and c in a Bell state, the remaining modes a and 
d must be in the |^+) state as well. 

Although maximally entangled states have never been produced experi- 
mentally, entanglement swapping might offer us a solution [9]. An entangled 
state with a large vacuum contribution (as produced by parametric down- 
conversion) can only give us randomly produced entanglement. However, if we 
use two such states and perform entanglement swapping, the Bell detection 
will act as a tell-tale that there were photons in the system. The question is 
whether this Bell detection is enough to ensure that an event-ready entangled 
state appears as a freely propagating wave-function. 

Entanglement swapping has been demonstrated experimentally by Pan et 
al. [10], using parametric down-conversion as the entanglement source. In the 
next section we briefly review the down-conversion mechanism and its role 
in the experiment of Pan et al. Subsequently, in section 3 we study whether 
entanglement swapping can give us event-ready entanglement. 

2 Parametric Down-Conversion and Entanglement 
Swapping 

In this section we review the mechanism of parametric down-conversion and 
the experimental realisation of entanglement swapping. In parametric down- 
conversion a crystal is pumped by a high-intensity laser, which we will treat 
classically (the parametric approximation) . The crystal is special in the sense 
that it has different refractive indices for horizontally and vertically polarised 
light. In the crystal, a photon from the pump is split into two photons with 
half the energy of the pump photon. Furthermore, the two photons have 
orthogonal polarisations. The outgoing modes of the crystal constitute two 
intersecting cones with orthogonal polarisations as depicted in Fig. 1. 

Due to the conservation of momentum, the two produced photons are 
always in opposite modes with respect to the central axis (determined by 




18 



Pietr Kok and Samuel L. Braunstein 




i) 



Fig. 1. A schematic representation of type II parametric down-conversion. A high- 
intensity laser pumps a non-linear crystal. With some probability a photon in the 
pump beam will be split into two photons with orthogonal polarisation | and | <->■) 
along the surface of the two respective cones. Depending on the optical axis of the 
crystal, the two cones are slightly tilted from each other. Selecting the spatial modes 
at the intersection of the two cones yields the outgoing state |0) + -I- O(f^). 

the direction of the pump). In the two spatial modes where the different 
polarisation cones intersect we can no longer infer the polarisation of the 
photons, and as a consequence the two photons become entangled in their 
polarisation. 

However, parametric down-converters do not produce Bell states [8,11,12]. 
They form a class of devices yielding Gaussian evolutions: 

jif) = U{t)\0) = exp[-ifHi/h]\0) , (6) 

with 

^E“l^bSj+H.c., (7) 

'ij 

where "H/ is the interaction Hamiltonian, a| a creation operator and yl^- the 
components of a (complex) symmetric matrix. Here, H.c. stands for the Her- 
mitian conjugate. If A is diagonal the evolution U corresponds to a collection 
of single-mode squeezers [13]. In the case of degenerate type II parametric 
down-conversion used to produce a photon-pair exhibiting polarisation en- 
tanglement, the interaction Hamiltonian in the rotating-wave approximation 
is 

ni = K{albl-albl) + H.c. , ( 8 ) 

with K a parameter which is determined by the strength of the pump and the 
coupling of the electro-magnetic field to the crystal. 

The outgoing state of the down-converter is then 

\1')ab = (l - jo, 0)ab + y)ab ~ \V, x)ab'j 





Event-Ready Entanglement 



19 



Bell detection 




Fig. 2. A schematic representation of the experimental setup for entanglement 
swapping as performed by Pan et al. The pump beam is reverted by a mirror in 
order to create two entangled photon-pairs in different directions (modes a, b, c and 
d). Modes b and c are sent into a beam-splitter (bs). A coincidence in the detectors 
Du and Du at the outgoing modes of the beam-splitter identify a Bell state. 
The undetected modes a and d are now in the Bell state as well. 



- \xy,xy)ab + \y^,x‘^)ab"j + , (9) 

with ^ <C 1, which is a function of k. Here, \x"^) denotes an x-polarised mode in 
a 2-photon Fock state (the case of two y-polarised or an x- and a y-polarised 
photon are treated similarly). 

For the experimental demonstration of entanglement swapping we need 
two independent Bell states. A Bell measurement on one half of either Bell 
state will then entangle the two remaining modes. The photons in these modes 
do not originate from a common source, i.e., they have never interacted. Yet 
they are now entangled. 

In the experiment conducted by Pan et al., instead of having two para- 
metric down-converters, one crystal was pumped twice in opposite directions 
(see Fig. 2). This way, a state which is equivalent to a state originating from 
two independent down-converters was obtained. In order to simplify our dis- 
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The Physical State 




PDCl pdc2 



Fig. 3. A schematic representation of the entanglement swapping setup. Two para- 
metric down-converters (pdc) create states which exhibit polarisation entangle- 
ment. One branch of each source is sent into a beam splitter (bs), after which the 
polarisation beam splitters (pbs) select particular polarisation settings. A coinci- 
dence in detectors Du and Du ideally identify the Bell state. However, since 
there is a possibility that one down-converter produces two photon-pairs while 
the other produces nothing, the detectors Du and D„ no longer constitute a Bell- 
detection, and the freely propagating physical state is no longer a pure Bell 
state. 



cussion, we will treat the experimental setup as if it consists of two separate 
down-converters (see Fig. 3). 

One spatial mode of either down-conversion state is sent into a beam- 
splitter, the output of which is detected. This constitutes the Bell measure- 
ment. In the case where both down-converters create a polarisation entangled 
photon-pair, a coincidence in the photo-detectors and Dy identify the an- 
tisymmetric Bell state |>F“) [14]. The outgoing state should then be the |>F“) 
Bell state as well, thus creating event-ready entanglement. 
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However, there are two problems. First, this is not a complete Bell mea- 
surement [15,16], i.e., it is not possible to identify all the four Bell states 
simultaneously. The consequence is that entanglement swapping occurs a 
quarter of the time (only |>F~) is identified). 

A second, and more serious, problem is that when we study coincidences 
between two down-converters, we need to take higher-order photon-pair pro- 
duction into account [see Eq. (9)]. For instance, one down-converter creates 
a photon-pair with probability |^|^. Two down-converters therefore create 
two photon-pairs with probability |^|^. However, this is roughly equal to the 
probability where one down-converter produces nothing (i.e., the vacuum 
|0)), while the other produces two photon-pairs. In the next section we will 
show that for this reason a coincidence in the detectors and no longer 
uniquely identifies a |<F“) Bell state. 



3 Event-Ready Entanglement 

In this section we first investigate the effect of double-pair production on 
the Bell measurement in the experimental setup depicted in Fig. 3. Subse- 
quently, we study the possibility of event-ready entanglement in the context 
of entanglement swapping. 

There is a possibility that a single down-converter produces a double 
photon-pair. This means that the two photons in the outgoing modes of 
the beam-splitter in Fig. 3 do not necessarily originate from different down- 
converters. We can therefore no longer interpret a detector-coincidence at the 
outgoing modes of the beam-splitter as a projection onto the |<F~) Bell state. 

Another way of looking at this is as follows: consider a two-photon polar- 
isation state. It is a vector in a Hilbert space generated by (for instance) the 
basis vectors \x,x), \x,y), \Vtx) and |y,y). The [S'”) Bell state is a super- 
position of these basis vectors. The key observation is that the two photons 
described in this Hilbert space occupy different spatial modes. When a two- 
photon state entering a 50:50 beam-splitter gives a two-fold coincidence, this 
state is projected onto the |iF~) Bell state. 

This Hilbert space should be clearly distinguished from a (truncated) Fock 
space. In the Fock space two photons can occupy the same spatial mode (see 
for example the state in Eq. (9)). As a consequence, the two input modes of 
a beam-splitter can be the vacuum |0) and a two-photon state (for instance 
\x^)) respectively. In this scenario a detector coincidence at the output of the 
beam-splitter is still possible, but it can not be interpreted as the projection 
of the incoming state |0,a;^) on the |<F“) Bell state (see Fig. 4). 

In the case of the entanglement swapping experiment, two photon-pairs 
are created either by one down-converter alone or both down-converters. This 
means that, conditioned on a detector coincidence, the state entering the 
beam-splitter is a superposition of two single-photon states plus the vacuum 
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a) identifies |^ ) b) no identification 




Fig. 4. A schematic representation of the Bell measurement in the entanglement 
swapping experiment. In Fig. 4a both incoming modes are occupied by a single 
photon. A coincidence in the detectors will then identify a projection onto the 
Bell state. In Fig. 4b, however, one input mode is the vacuum, whereas the 
other is populated by two photons. In this case a detector coincidence cannot be 
interpreted as a projection onto the Bell state. Both instances a) and b) occur 
in the entanglement swapping experiment. 

and a two-photon state. In this setup a detector coincidence therefore does 
not identify the Bell state. 

In Fig. 3 we have added two polarisation beam-splitters in the outgoing 
modes of the beam-splitter. This allows us to condition the outgoing state in 
modes a and d on the different polarisation settings. It should be noted that 
for a Bell detection depicted in Fig. 4a we do not need a polarisation sensitive 
measurement. However, since the detector coincidence is no longer a Bell 
measurement it will be convenient to distinguish between the four different 
polarisation settings at the output modes of the beam-splitter (x,y), 

(y,x) and {y,y)]. 

It turns out that the four outgoing states conditioned on the four different 
polarisation settings have a remarkably simple form: 

l</-(...))ad = (|0,y")-|y",0))/V2, 

l</*(a;.y))ad = (|0,a;y) - \y,x) + \x,y) - \xy,Q))/2, 

|0(j/,x))ad = (|0,a;2/) -b \y,x) - \x,y) - \xy,Q))/2, 

kto,y))ad = (|0,x")-|a:2,0))/V2. (10) 

These states can also be obtained by sending \y, y), \y, x), \x, y) and \x, x) into 
a 50:50 beam-splitter respectively (see Fig. 5). When no distinction between 
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|0, xy) - \y, x) + \x, y) - \xy, 0) |0, y^) - |y^ 0) 





\y) \y) 



Fig. 5. The four (nnnormalised) outgoing states from a 50:50 beam-splitter condi- 
tioned on the four input states \x,x), \x,y), \y,x) and \y,y). The outgoing states 
correspond to the four possible outgoing states of entanglement swapping. 



the four possible polarisation settings in the two-fold detector-coincidence is 
made [10,8], the state of the two remaining (undetected) modes (the physical 
state) will be a mixture p of the four states in Eqs. (10). 

Using a technical mathematical criterion based on the partial transpose 
of the outgoing state p [17,18], it can be shown that p is entangled {p satisfies 
this entanglement criterion). However, it can not be used for event-ready 
detections of polarisation entanglement since the states in Eqs. (10) are not 
of the form of Eq. (3). More specifically, these states are nowhere near the 
maximally entangled states we wanted to create. 

We now ask the question whether entanglement swapping can still be used 
to create event-ready entanglement. Observe that the experiment conducted 
by Pan et al. closely resembles the experimental realisation of quantum tele- 
portation by Bouwmeester et al. [4]. There, one of the states produced by 
parametric down-conversion was used to prepare a single-photon state which 
was then teleported. However, due to the double pair-production, the tele- 
ported state was a mixture of the teleported photon and the vacuum. 

In Ref. [11], Braunstein and Kimble presented three possible ways to im- 
prove the experimental setup of quantum teleportation in order to minimise 
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the unwanted double pairs created in a single down-converter. Since the ex- 
perimental setup considered here closely resembles the setup of the telepor- 
tation experiment, one might expect that the remedy given by Braunstein 
and Kimble will be effective here as well. 

However, this is not the case. The first approach was to make sure that in 
one of the outgoing modes (the state-preparation mode) there was only one 
photon present. In order to achieve this, a detector cascade was suggested 
which, upon a two-fold detector coincidence, would reveal a two-photon state. 
Since in the entanglement swapping experiment there is no conditional mea- 
surement of the outgoing modes (apart from the Bell measurement), this 
approach doesn’t work here. 

The second approach is to differentiate between the coupling of the two 
down-converters by lowering the intensity of the pump of one of them. How- 
ever, since the setup of the entanglement swapping experiment is symmetric, 
it can be easily seen that the (unnormalised) outgoing states become 



!</>(.,.)) « ?||o,2/2)-c?l2/^o)+o(a, 

\<l^ix,y)) oc ^2 |0>a;y) -^ 1^2 (ly.a;) - |a;, 2/)) -^1 |x?/,0) -kO(C^), 

\<l>iv,x)) oc Cl |0,a;2/) -kCi6 (|2/,x) - |a;,2/)) |xy,0) -kO(C^), 

!</>(,,,)) cx C 2 |0,a:2)-C||a:^0) + O(a. (11) 



Varying Ci and C 2 will not allow us to create any of the states which have the 
form of Eq. (3). 

The third approach presented by Braunstein and Kimble involved the 
use of quantum non-demolition (qnd) measurements. In the case of entan- 
glement swapping, we should be able to identify whether there are photons 
in the outgoing modes with such qnd detectors. If we can tell that both 
outgoing modes are populated by a photon without destroying the non-local 
correlations (i.e., without gaining information about the polarisations of the 
photons), we can create event-ready entanglement. However, such qnd de- 
tectors correspond to technology not yet available for optical photons (for 
QND detections of microwave photons see Ref. [19]). 

Can we turn any of the states in Eq. (10) into the form of Eq. (3) by 
other means? Additional photon sources would take us beyond the entangle- 
ment swapping protocol and we will not consider them here. Alternatively, 
we might be able to use a linear interferometer to obtain event-ready entan- 
glement. First, observe that instead of taking Eqs. (10) as the input of the 
interferometer, we can use \x,x), \x,y), \y,x) or \y,y), since these states are 
obtained from Eqs. (10) by means of a (unitary) beam-splitter operation. If 
the linear interferometer described above does not exist for these separable 
states, neither will it for the outgoing states of the entanglement swapping 
protocol, since these states only differ by a unitary transformation. 
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Next, suppose we have a linear interferometer described by the unitary 
matrix U [20], which transforms the creation operators of the electro-magnetic 
field according to 

cii y ^ hj , 
j 

3 

where the Uij are the components of U and i,j enumerate both the modes 
and polarisations. There is no mixing between the creation and annihilation 
operators, because photons do not interact with each other. 

In order to create event-ready entanglement, the creation operators which 
give rise to the states obtained by the entanglement swapping protocol (for 
instance \x,y) = aj,6j|0)) should be transformed according to 

aUl (cUl - cldl)/V2 . (13) 

Relabel the modes Ux, ay, bx and by as a\ to 04, respectively. Without loss of 
generality (and leaving the normalisation aside for the moment) we can then 
write 

a\a\ b\b\ - b\b\ . (14) 

Substituting Eq. (12) into Eq. (14) generates ten equations for eight variables 

uy. 

^ 11^21 — U12U22 — '^13'^23 — '*^ 14^24 — d , 

^11^23 + ■*^13'*^2l) “ (^H'^24 + ^^14^^21) “ d, 

(U12U23 + ^134422) = (U12U24 + U44U22) = d, 

(ul^U*22 + ul2U*2,) =1, 

+ =-l. (15) 

It can be easily verified that there are no solutions for the u*j which satisfy 
these ten equations simultaneously. This means that there is no passive inter- 
ferometer which transforms the states resulting from entanglement swapping 
into event-ready entangled states. 

4 Conclusions 

We speak of event-ready entanglement when upon the production of a max- 
imally polarisation entangled state the source flashes a red light (or gives 
a similar macroscopic indication). We have studied the possibility of event- 
ready entanglement preparation by means of entanglement swapping as per- 
formed by Pan et al. [10]. This was suggested by Zukowski, Zeilinger, Horne 
and Ekert [8] and Pavicic [9]. 
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In entanglement swapping, we perform a Bell measurement on two parts 
of two (maximally) entangled states, leaving the two undetected parts en- 
tangled. In general, these parts have never interacted. In the experiment of 
Pan et al, the entangled states are produced by means of parametric down- 
conversion. Due to higher order corrections in the down-converters, the phys- 
ical state leaving the entanglement swapping apparatus is a random mixture 
of four states. These states correspond to the four possible polarisation set- 
tings in the Bell measurement. They are equivalent to the outgoing state of a 
50:50 beam-splitter conditioned on the four possible input states \x, x), \x, y), 
|j/, a;) and \y,y). However, these states are not the ones we were looking for. 

With respect to a related experiment involving quantum teleportation [4] , 
it has been described [11] how to enhance the fidelity of the outgoing state. 
However, these methods fail here, as well as the application of a linear inter- 
ferometer with passive elements. We need at least a quantum non-demolition 
measurement or a quantum computer of some kind to turn the outgoing states 
of the entanglement swapping experiment into event-ready entanglement. 

Appendix: Transformation of Maximally Entangled 
States 

In this appendix we will show that each maximally entangled two-system 
state can be transformed into any other by a local unitary transformation 
on one subsystem alone. For example, every two-photon polarisation Bell 
state can be transformed into any other maximally polarisation entangled 
two-photon state by a linear optical transformation on one mode. 

We will treat this in a formal way by considering an arbitrary maximally 
entangled state of two A^-level systems in the Schmidt decomposition: 

1 ^ 

• ( 16 ) 

i.e., a state with equal amplitudes on all possible branches. Here {ju-i)} and 
{jmi)} form two orthonormal bases of the subsystems. We can obtain any 
maximally entangled state by applying the (bi-local) unitary transformation 
Ux®U 2 - This means that each maximally entangled state can be transformed 
into any other by a pair of local unitary transformations on each of the 
subsystems. We will now show that any two maximally entangled states jV') 
and |^/>') are connected by a local unitary transformation on one subsystem 
alone: 

\i)) = U®tW) . (17) 

First, we show that any transformation Ui®U 2 on a, particular maximally 
entangled state \(j)) can be written as H ® 11 [())), where V = C/iC/J. Take 

1 ^ 



( 18 ) 
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For any U we have 

u®u*\4)) = 10) . 

Proof. Using the completeness relation 

N 



\nk){nk\ = ]1 






on both subsystems we have 

N 



|nfc,m/)(nfc,m/| (C/0 1) (10 t/*) I 






By writing out \(j)) explicitly according to Eq. (18) we obtain 



N 



C/(g)C/*|0) = ^ UkjUfj\nk,mi) 

^ 3,k,l=l 

1 ^ 

= 7 ^ E {UU'^)ki\nk,mi) 






= - 7 = E 

fc,i=l 



(19) 



( 20 ) 



(21) 



(22) 



which is equal to |0). □ 

Next, using the relation (19) we will show that every transformation U\ 0 
U 2 acting on |0) is equivalent to a transformation U ® 11 acting on |0), where 
V = C/iC/J. 

Proof. From Eq. (19) we obtain the equation 

(l®[/^)(t/®C/*)|0) = (C/®ll)|0), (23) 

which immediately gives us 

U 0 il|0) = 1 0 U^|0) and 

[/'^(g)]l|0) = 1(g) C/|0) . (24) 



From Eq. (24) we obtain 



Ui 0 C/ 2 I 0 ) = (Ui 0 1) (1 0 U 2 ) 10) 
= (Ui 0 1) (Uj 0 1) 10) 
= UiU^ (g) 1|0) . 



( 25 ) 
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Similarly, 

Ui (g) U2\(t>) = 1 (g) c/f U2\(j)) . (26) 

We therefore obtain that Ui (g) U2\4>) is equal to F (g) t\4>) with V = C/iC/J, 
and similarly that it is equal to 11 (g) y'|( 5 i) with V' = U^U2- □ 

Since every maximally entangled stated can be obtained by applying U\ (g) 
U2 to |(/)), two maximally entangled states |^) and \il)') can be transformed 
into any other by choosing 



\ii) = [/ (g) 11 1^) 

W) = V®t\<j>) , (27) 

which gives 

1^) = c/yt ^ ]1 1 ^') . (28) 

Thus each maximally entangled two-system state can be obtained from any 
other by means of a local unitary transformation on one subsystem alone. 

Returning to a pair of maximally polarisation entangled photons, we know 
that any unitary transformation on a single optical mode consisting of sin- 
gle photons can be viewed as a combination of a polarisation rotation and 
a polarisation dependent phase shift of that mode. We can easily perform 
such operations with linear optical elements. If we can create one maximally 
entangled state, we can create any maximally entangled state. It is therefore 
sufficient to restrict our discussion to, for example, the Bell state. 
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Abstract. The processes of radiation damping and decoherence in Quantum Elec- 
trodynamics are studied from an open system’s point of view. Employing functional 
techniques of field theory, the degrees of freedom of the radiation field are eliminated 
to obtain the influence phase functional which describes the reduced dynamics of 
the matter variables. The general theory is applied to the dynamics of a single 
electron in the radiation field. From a study of the wave packet dynamics a quanti- 
tative measure for the degree of decoherence, the decoherence function, is deduced. 
The latter is shown to describe the emergence of decoherence through the emission 
of bremsstrahlung caused by the relative motion of interfering wave packets. It is 
argued that this mechanism is the most fundamental process in Quantum Elec- 
trodynamics leading to the destruction of coherence, since it dominates for short 
times and because it is at work even in the electromagnetic field vacuum at zero 
temperature. It turns out that decoherence trough bremsstrahlung is very small for 
single electrons but extremely large for superpositions of many-particle states. 



1 Introduction 

Decoherence may be defined as the (partial) destruction of quantum coher- 
ence through the interaction of a quantum mechanical system with its sur- 
roundings. In the theoretical analysis decoherence can be studied with the 
help of simple microscopic models which describe, for example, the interaction 
of a quantum mechanical system with a collection of an infinite number of 
harmonic oscillators, representing the environmental degrees of freedom [1,2]. 
In an open system’s approach to decoherence one derives dynamic equations 
for the reduced density matrix [3] which yields the state of the system of 
interest as it is obtained from an average over the degrees of freedom of the 
environment and the resulting loss of information on the entangled state of 
the combined total system. The strong suppression of coherence can then be 
explained by showing that the reduced density matrix equation leads to an 
extremely rapid transitions of a coherent superposition to an incoherent sta- 
tistical mixture [4,5]. For certain superpositions the associated decoherence 
time scale is often found to be smaller than the corresponding relaxation or 
damping time by many orders of magnitude. This is a signature for the fun- 
damental distinction between the notions of decoherence and of dissipation. 
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A series of interesting experimental investigations of decoherence have been 
performed as, for example, experiments on Schrodinger cat states of a cavity 
field mode [6] and on single trapped ions in a controllable environment [7]. 

If one considers the coherence of charged matter, it is the electromagnetic 
field which plays the role of the environment. It is the purpose of this paper to 
study the emergence of decoherence processes in Quantum Electrodynamics 
(QED) from an open system’s point of view, that is by an elimination of 
the degrees of freedom of the radiation field. An appropriate technique to 
achieve this goal is the use of functional methods from field theory. In section 
2 we combine these methods with a super-operator approach to derive an 
exact, relativistic representation for the reduced density matrix of the matter 
degrees of freedom. This representation involves an influence phase functional 
that completely describes the influence of the electromagnetic radiation field 
on the matter dynamics. The influence phase functional may be viewed as 
a super-operator representation of the Feynman-Vernon influence phase [1] 
which is usually obtained with the help of path integral techniques. 

In section 3 we treat the problem of a single electron in the radiation 
field within the non-relativistic approximation. Starting from the influence 
phase functional, we formulate the reduced electron motion in terms of a 
path integral which involves an effective action functional. The corresponding 
classical equations of motion are demonstrated to yield the Abraham-Lorentz 
equation describing the radiation damping of the electron motion. In addition, 
the influence phase is shown to lead to a decoherence function which provides 
a measure for the degree of decoherence. 

The general theory will be illustrated with the help of two examples, 
namely a free electron (section 4) and an electron moving in a harmonic po- 
tential (section 5) . For both cases an analytical expression for the decoherence 
function is found, which describes how the radiation field affects the electron 
coherence. 

We shall use the obtained expressions to investigate in detail the time- 
evolution of Gaussian wave packets. We study the influence of the radia- 
tion field on the interference pattern which results from the collision of two 
moving wave packets of a coherent superposition. It turns out that the ba- 
sic mechanism leading to the decoherence of matter waves is the emission 
of bremsstrahlung through the moving wave packets. The resultant picture 
of decoherence is shown to yield expressions for the decoherence time and 
length scales which differ substantially from the conventional estimates de- 
rived from the prominent Caldeira-Leggett master equation. In particular, it 
will be shown that a superposition of two wave packets with zero velocity 
does not decohere and, thus, the usual picture of decoherence as a decay of 
the off-diagonal peaks in the corresponding density matrix does not apply to 
decoherence through bremsstrahlung. 

We investigate in section 6 the possibility of the destruction of coherence 
of the superposition of many-particle states. It will be argued that, while the 
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decoherence effect is small for single electrons at non-relativistic speed, it is 
drastically amplified for certain superpositions of many-particle states. 
Finally, we draw our conclusions in section 7. 



2 Reduced Density Matrix of the Matter Degrees of 
Freedom 

Our aim is to eliminate the variables of the electromagnetic radiation field 
to obtain an exact representation for the reduced density matrix pm of the 
matter degrees of freedom. The starting point will be the following formal 
equation which relates the density matrix Pm(tf) of the matter at some final 
time tf to the density matrix p(ti) of the combined matter-field system at 
some initial time ti, 



Pm{tf) = try 



|r<_ exp 




The Liouville super-operator C{x) is defined as 



( 1 ) 



£{x)p = -i[n{x),p], 



(2) 



where 'H(x) denotes the Hamiltonian density. Space-time coordinates are 
written as x^ = (x^,x) = (t,x), where the speed of light c is set equal 
to 1. All fields are taken to be in the interaction picture and T<_ indicates 
the chronological time-ordering of the interaction picture fields, while tr f de- 
notes the trace over the variables of the radiation field. Setting h = c = 1 we 
shall use here Heaviside-Lorentz units such that the fine structure constant 
is given by 



a = 



Airhc 



1 

137' 



(3) 



To be specific we choose the Coulomb gauge in the following which means 
that the Hamiltonian density takes the form [8-10] 



Hix) = Hcix) + Titrix). 



(4) 



Here, 



'H-trix) = j^{x)Af,{x) (5) 

represents the density of the interaction of the matter current density j^(x) 
with the transversal radiation field. 



A>^{x) = {0,A{x)), V-A{x) = 0, 



(6) 



= ^f{x)A°{x) 



1 r 3 f{x°,x)f{x°,y) 

2j ^ 4:Tr\x-y\ 



and 



(7) 
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is the Coulomb energy density such that 

Hc(x^) = ^Jd^xJ d^y 



f{x°,x)f{x^,y) 
4 tt\x — y\ 



(8) 



is the instantaneous Coulomb energy. Note that we use here the convention 
that the electron charge e is included in the current density j^(x) of the 
matter. 

Our first step is a decomposition of chronological time-ordering operator 
T<_ into a time-ordering operator T^_ for the matter current and a time- 
ordering operator for the electromagnetic field, 

T^=Tl_T^. (9) 



This enables one to write Eq. (1) as 



Pm{tf) = Tl_ try T^exp 



d'^x{Cc{x) + Ct^{x)) 



P{U) 



(10) 



where we have introduced the Liouville super-operators for the densities of 
the Coulomb field and of the transversal field, 



Cc{x)p= -\[Uc{.x),p], Ctr{x)p = -i[j^{x)A^{x),p\. (11) 



The currents commute under the time-ordering T|_. We may therefore 
treat them formally as commuting c-number fields under the time-ordering 
symbol. Since the super-operator Cc{x) only contains matter variables, the 
corresponding contribution can be pulled out of the trace. Hence, we have 



Pm{tf) = Tl_ exp 



rtf 



d'^xCc{x) 



tr/ < exp 



ftf 



d^’^xCtrix) 



p{u) 



( 12 ) 

We now proceed by eliminating the time-ordering of the H- fields. With 
the help of the Wick-theorem we get 



T^exp 



rtf 



d'^xCtrix) 



(13) 



exp 



^ d‘^x'[Ctr{x),Ctr{x')]9{t-t') 



exp 



r^f 



d'^xCtrix) 



In order to determine the commutator of the Liouville super-operators we 
invoke the Jacobi identity which yields for an arbitrary test density p, 



[Ctr{x),Ctr{x')]p = Ctr{x)Ctr{x')p - Ct,{x') Ctr{x) p 

= -[UtAx), [HtAx'), p]] + [UtAx'), [HtAx), p\] 

= -[[ntAx),'HtAx')],p\- (14) 
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The commutator of the transversal energy densities may be simplified to read 






since the contribution involving the commutator of the currents vanishes by 
virtue of the time-ordering operator T^. Thus, it follows from Eqs. (14) and 
(15) that the commutator of the Liouville super-operators may be written as 

[Ctv{x),Ctr{x')]p= -[A^{x),A^{x')][f{x)f{x'),p]. (16) 

It is useful to introduce current super-operators J+{x) and J-{x) by means 
of 

Jl{x)p = f{x)p, Jt{x)p = pf{x). (17) 

Thus, J+{x) is defined to be the current density acting from the left, while 
J-{x) acts from the right on an arbitrary density. With the help of these 
definitions we may write the commutator of the Liouville super-operators as 



[£tr(a;),/ltr(a;')] = ~ [A^l.{x),A^,{x')]J^{x)J'^{x') 
+ \A^{x),A^{x')]J‘^{x)J'L{x'). 

Inserting this result into Eq. (13), we can write Eq. (12) as 



Pm{tf) = Tl_ ( exp 



Ut. 

1 



d'^xCc{x) 
Af 



d^x^ d"^x'6{t — t')[Af^{x) , A^{x')]J!^{x) J'^{x') 

+ ^ d^‘a;'6<(t - t')[^/i(a;), ^!y(a;')]</-(a;)-^-(a;') 



•try < exp 



d^xCf[{x) 



p{u) 



(18) 



This is an exact formal representation for the reduced density matrix of the 
matter variables. Note that the time-ordering of the radiation degrees of 
freedom has been removed and that they enter Eq. (18) only through the 
functional 

Af 



W[J+, J-] = try < exp 



d*xCu{x)p{U) 



(19) 



since the commutator of the A-fields is a c-number function. 



3 The Influence Phase Functional of QED 

The functional (19) involves an average over the field variables with respect 
to the initial state p(h) of the combined matter-field system. It therefore 
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contains all correlations in the initial state of the total system. Here, we 
are interested in the destruction of coherence. Our central goal is thus to 
investigate how correlations are built up through the interaction between 
matter and radiation field. We therefore consider now an initial state of low 
entropy which is given by a product state of the form 

p{k) = PmiU)® Pf, (20) 

where Pm(ti) is the density matrix of the matter at the initial time and 
the density matrix of the radiation field describes an equilibrium state at 
temperature T, 

P/ = ^exp(-/3iJ/). (21) 

Here, Hf denotes the Hamiltonian of the free radiation field and the quantity 
Zf = tr f[exp{—PHf)] is the partition function with f3 = l/ksT. In the 
following we shall denote by 

{0)f = tTf{0pf} (22) 

the average of some quantity O with respect to the thermal equilibrium state 

( 21 ). 

The influence of the special choice (20) for the initial condition can be 
eliminated by pushing ti — >■ — oo and by switching on the interaction adiabat- 
ically. This is the usual procedure used in Quantum Field Theory in order to 
define asymptotic states and the S'-matrix. The matter and the field variables 
are then described as m-fields, obeying free field equations with renormalized 
mass. These fields generate physical one-particle states from the interacting 
ground state. 

For an arbitrary initial condition p(ti) the functional 1F[J+, J_] can be 
determined, for example, by means of a cumulant expansion. Since the initial 
state (20) is Gaussian with respect to the field variables and since the Liouville 
super-operator Ctr{x) is linear in the radiation field, the cumulant expansion 
terminates after the second order term. In addition, a linear term does not 
appear in the expansion because of {A^{x)) f = 0. Thus we immediately 
obtain 

1F[J+, J_] = exp PmiU). (23) 

Inserting the definition for the Liouville super-operator Ctr{x) into the expo- 
nent of this expression one finds after some algebra, 

^ d‘^x'{Ctr{x)Ctr{x'))fPm 

= ~\J^ {['Htr{x),['Htr{x'),Pm^ Pf]]} 




Radiation Damping and Decoherence in Quantum Electrodynamics 



37 



= -\J^ d^x' [{A:,{x')A^{x))fJ!l{x)Jl{x') 

+{A^ix)A,{x'))frix)r{x') 
-{A^{x')Af,{x)) fJ^{x)J'l{x') 
-{Af^{x)A^{x'))fJ^{x)J'l{x')] Pm- 

On using this result together with Eq. (23), Eq. (18) can be cast into the 
form, 

Pmitf) =T^(^exp d‘‘x£c(a;) (24) 

+ ^y J {e{t-t')[Af,{x),A^{x')] 

+{A^{x')Af,{x))f)j^{x)Jl{x') 

+ {0{t - t')[Afj_{x),A„{x')] 

-{A,{x)A,{x'))f)r{x)r_{x') 

+{A^{x')A^{x))fJ^{x)J'L{x') 

+ {Af,{x)Ai,{x')) fJ^{x)J'^{x')] ^Pm(ti). 

At this stage it is useful to introduce a new notation for the correlation 
functions of the electromagnetic field, namely the Feynman propagator and 
its complex conjugated (T_> denotes the anti-chronological time-ordering), 

iDf{x - x')fj,i, = (T<_(A^(x)A,,(x')))/ 

= - t')[A^{x),A^{x')] + (Aj,(x')A^(x))/, 

iD*p{x - x')^^ = ~{T^{A^{x)Ai,{x'))) f 

= 9{t - t')[A^{x),A^{x')] - {A^{x)A^{x'))f, (25) 

as well as the two-point correlation functions 

D+{x - x')^u = {A^{x)A^,{x'))f, 

L)_(x - x')^^ = {A^{x')A^{x))f. (26) 

As is easily verified these functions are related through 

-iDf{x - x')fj,u + iD*f{x - x')fj,^ + D+{x - x')^u + D^{x - x')^i/ = 0. (27) 

With the help of this notation the density matrix of the matter can now be 
written as follows. 



Pm{tf) — TI_ 



exp 



d^x£c(x) 



(28) 




38 



Heinz-Peter Breuer and Francesco Petruccione 



+ ^ J d'^x J d^x' {-iDp^x - x')f_u,J^{x)J'^{x') 

+iDp{x - x')^yJ^{x)J'^{x) 

+D-{x - x')^^Jp{x)J'i{x') 

+ D+{x - x')pi,J^{x)Jp{x')]]) Pm{ti)- 

This equation provides an exact representation for the matter density matrix 
which takes on the desired form: It involves the electromagnetic field variables 
only through the various two-point correlation functions introduced above. 
One observes that the dynamics of the matter variables is given by a time- 
ordered exponential function whose exponent is a bilinear functional of the 
current super-operators J±{x). Formally we may write Eq. (28) as 

= Tl_ exp (i^[J+, J-]) Pm{U), (29) 



where we have introduced an influence phase functional 



i^[J+,J_]= / d^xCc{x) 




(30) 



X [-\Df{x - x)^^Jp{x)J'f{x) + iDp{x - x')fj,„Jfl{x)Jf{x') 
+D_{x - x')fj,„Jp{x)Jf{x') + D+{x - x')fj,„J^{x)Jf_{x')} . 



It should be remarked that the influence phase ^[J+, J-] is both a functional 
of the quantities J± (x) and a super-operator which acts in the space of den- 
sity matrices of the matter degrees of freedom. There are several alternative 
methods which could be used to arrive at an expression of the form (30) 
as, for example, path integral techniques [1] or Schwinger’s closed time-path 
method [11]. The expression (30) for the influence phase functional has been 
given in Ref. [12] without the Coulomb term and for the special case of zero 
temperature. In our derivation we have combined super-operator techniques 
with methods from field theory, which seems to be the most direct way to 
obtain a representation of the reduced density matrix. 

For the study of decoherence phenomena another equivalent formula for 
the influence phase functional will be useful. To this end we define the com- 
mutator function 



D{x - x')^y = i[A;x(x), Aflx')] 

= \{D+{x - x')f^^ - D_{x - x')f_,flj (31) 

and the anti-commutator function 



Di{x-x')p„ = {{Ap{x),Ai,{x')})f 

— d-|-(x X )piy D—{x X )piy 



(32) 
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Of course, the previously introduced correlation functions may be expressed 
in terms of D{x — x')^^ and D\{x — x')^^, 

D+{x-x')^i, = ^Di{x - x')f^„ - ^D{x - x')f^,,, (33) 

D_{x-x')f^,, = ^Di{x - x')f^i, + ^D{x - x')f^,,, (34) 

iDp{x - x')f^„ = ^Di{x - x')fj.u - ^sign(t - t')D{x - x')^^, (35) 

-iD*p{x - x')^i, = ]^Di{x - x')fj,y + ^sign(t - t')D{x - x')^„. (36) 

Correspondingly, we define a commutator super-operator Jc{x) and an anti- 
commutator super-operator Ja{x) by means of 

Ji^{x)p ^ [r{x),p], Ji^{x)p ^ {f{x),p}, (37) 

which are related to the previously introduced super-operators J±{x) by 

Jl!{x) = Jt{x) - r{x), Jl^ix) = J^(x) + r{x). (38) 

In terms of these quantities the influence phase functional may now be written 
as 

i<P[Jc,Ja] = J d''^xCc{x) (39) 

-^-D,{x-x'),,Jl!ix)j:ix')Y 

This form of the influence phase functional will be particularly useful later 
on. It represents the influence of the radiation field on the matter dynamics 
in terms of the two fundamental 2-point correlation functions D(x — x') and 
Di{x — x'). Note that the double space-time integral in Eq. (39) is already 
a time-ordered integral since the integration over t' = x''^ extends over the 
time interval from ti to t = x^ . 

For a physical discussion of these results it may be instructive to compare 
Eq. (28) with the structure of a Markovian quantum master equation in 
Lindblad form [3], 




where Hm generates the coherent evolution and the Ai denote a set of opera- 
tors, the Lindblad operators, labeled by some index i. One observes that the 
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terms of the influence phase functional involving the current super-operators 
in the combinations J+J_ and J_J+ correspond to the gain terms in the 
Lindblad equation having the form These terms may be interpreted 

as describing the back action on the reduced system of the matter degrees 
of freedom induced by “real” processes in which photons are absorbed or 
emitted. The presence of these terms leads to a transformation of pure states 
into statistical mixtures. Namely, if we disregard the terms containing the 
combinations J+J_ and J_ J+ the remaining expression takes the form 

Pm{tf) « (41) 



where 

— i^ d'^x'Hcix) (42) 

It It ~ ■ 

Eq. (41) shows that the contributions involving the Feynman propagators 
and the combinations J+J+ and J_ J_ of super-operators preserve the purity 
of states [12]. Recall that all correlations functions have been defined in terms 
of the transversal radiation held. We may turn to the covariant form of the 
correlation functions if we replace at the same time the current density by 
its transversal component The expression (42) is then seen to contain the 
vacuum-to- vacuum amplitude A[j] of the electromagnetic held in the presence 
of a classical, transversal current density j([.(a;) [13], 



U{tf,ti) = Tf_exp 



A[j] = exp 




J d^^x'Dpix - . 



(43) 



With the help of the decomposition (35) of the Feynman propagator into a 
real and an imaginary part we And 



A[j] = exp 






(44) 



The vacuum-to-vacuum amplitude is thus represented in terms of a complex 
action functional with the real part 



= ^Jd‘^xJ d‘^x'sign{t - t')D{x - x')f,^j^^{x)j!^^{x'), (45) 



and with the imaginary part 



The imaginary part yields the probability that no photon is emitted by 
the current j([.. 



lR[j]]2 = exp(-25(2)y 



(47) 
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In covariant form we have 



D{x 



x') 



— 



1 

27T 



sign(t — t')S[{x 






(48) 



and, hence, 



S'(i) = j d‘‘x j d‘^x' 5[{x - x'f]f^{x)j'^^{x'). (49) 

This is the classical Feynman- Wheeler action. It describes the classical mo- 
tion of a system of charged particles by means of a non-local action which 
arises after the elimination of the degrees of freedom of the electromagnetic 
radiation field. In the following sections we will demonstrate that it is just 
the imaginary part which leads to the destruction of coherence of the 
matter degrees of freedom. 

4 The Interaction of a Single Electron with the 
Radiation Field 

In this section we shall apply the foregoing general theory to the case of a 
single electron interacting with the radiation field where we confine ourselves 
to the non-relativistic approximation. It will be seen that this simple case 
already contains the basic physical mechanism leading to decoherence. 

4.1 Representation of the Electron Density Matrix in the 
Non-Relativistic Approximation 

The starting point will be the representation (29) for the reduced matter 
density with expression (39) for the influence phase functional <?. It must be 
remembered that the correlation functions D{x — x')^^ and Di{x — x')^i, have 
been defined in terms of the transversal radiation field using Coulomb gauge 
and that they thus involve projections onto the transversal component. In 
fact, we have the replacements, 

D{x - x')^^ — ;> D{x - x')^j = - D{x - x') 

for the commutator functions, and 

Di{x - x')^t, — Di{x - x')ij = + (^Sij - Di{x - x') 

for the anti-commutator function, where 

n, I 



[exp {—\k{x — x')) — exp (ik(x — a;'))] , 
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and 

Di{x 



d^k 

2(27r)^w 



[exp {—ik{x 



x')) + exp (ik{x — x'))] coth (f3u!j2) , 



with the notation k^ = {uj,k) = (|fe|,fc) for the components of the wave 
vector. It should be noted that the commutator function is independent of the 
temperature, while the anti-commutator function does depend on T through 
the factor coth(/3w/2) = 1-1- 2N{ui), where N{uj) is the average number of 
photons in a mode with frequency u. Hence, invoking the non-relativistic 
(dipole) approximation we may replace 



D{x-x')i 



= SijD{t-t') = 6i. 



dujJ{u;)smuj{t — t'), (50) 



and 



Di{x - x')ij — )> Di(t - t')ij = 6ijDi(t - t') (51) 

roo 

= Sij / dwJ(w) coth (/3w/2) cos w(t — t'), 

Jo 

where we have introduced the spectral density 

J(w) = — ^a;6>(12 - w), (52) 

with some ultraviolet cutoff ^2 (see below) . It is important to stress here that 
the spectral density increases with the first power of the frequency uj. Had 
we used dipole coupling -ex-E of the electron coordinate x to the electric 
field strength E, the corresponding spectral density would be proportional 
to the third power of the frequency. This means that the coupling to the 
radiation field in the dipole approximation may be described as a special 
case of the famous Caldeira-Leggett model [2] and that in the language of 
the theory of quantum Brownian motion [14] the radiation field constitutes a 
super-Ohmic environment [15,16]. Note also that we now include the factor 
into the definition of the correlation function. Within the non-relativistic 
approximation we may thus replace the current density by 



j{t,x) — )> ’^S{x-x{t)) + 6{x-x{t))'^, 



(53) 



where p{t) and x{t) denote the momentum and position operator of the 
electron in the interaction picture with respect to the Hamiltonian 



= 



P 

2m 



■V(x) 



( 54 ) 



for the electron, V{x) being some external potential. 
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We are thus led to the following non-relativistic approximation of Eq. (29), 



Pm{tf) = exp 



Jt. [2 mm 

2 mm 



(55) 






This equation represents the density matrix (neglecting the spin degree of 
freedom) for a single electron interacting with the radiation field at tempera- 
ture T. In accordance with the definitions (37) and (38) Pc is a commutator 
super-operator and Pa an anti-commutator super-operator. In the theory of 
quantum Brownian motion the function D{t — t') is called the dissipation 
kernel, whereas Di{t — t') is referred to as noise kernel. 



4.2 The Path Integral Representation 

The reduced density matrix given in Eq. (55) admits an equivalent path 
integral representation [14] which may be written as follows, 

Pm(^/,^J,t/) — ^ j ^ f 1 ^ f } I f 1 Pm(.^i , ti^ , (56) 

with the propagator function 

J{xf,x'f,tf,x„x[,U) = / D£cD£c'exp{i (S'm[x] - Sm[x']) + i<^[x,x']} . 

(57) 

This is a double path integral which is to be extended over all paths x{f) and 
x'(t) with the boundary conditions 

x'{ti) = x[, x'{tf) = x'f, x{U) = Xi, x{tf) = Xf. (58) 
Sm\x] denotes the action functional for the electron. 




while the influence phase functional becomes, 

\<l>[x,x'] = dt' — t') (x{t) — x'(t)) (x{t') + x'{t')) 

-^Di{t-t') {x{t)-x{t)) (®(t') (t'))| -(60) 

We define the new variables 



q = X — x' , 



r = -{x + x'), 



(61) 
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and set, for simplicity, the initial time equal to zero, U = 0. We may then 
write Eq. (56) as 

Pm{rf,qf,tf)=J d^ri J d^q^J{rf,qf,tf;ri,q,)pm{r^,qi,0). (62) 

The propagator function 

J{rf,qf,tf;ri,qi) = J Dr J Dqexp{U[r, q]} (63) 

is a double path integral over all path r{t), q{t) satisfying the boundary 
conditions, 

r(0)=r„ r{tf)=rf, q{Q) = q^, q{tf) = qj. (64) 

The weight factor for the paths r{t), q{t) is defined in terms of an effective 
action A functional, 

f 1 1 

A[r,q]= / dt [mrq — V{r + -q) + V{r — -q) 

Jo \ 2 2 

+ f dt f dt'9{t — t')D{t — t')q{t)r{t') 

Jo Jo 

Di{t - t')q{t)q{t'). (65) 

4 Jo Jo 

The first variation of A is found to be 

SA=-J^ dti^Sq{t) mr{t) + ^Vr{V{r+^q) + V{r-^q)) 

+ dt'D(t-t')r(t') + dt'Di{t-t')q{t') 

+5r{t) mq{t) + 2Vg(E(r + ^q) + E(r - ^q)) 

|> (66) 

which leads to the classical equations of motion, 

mr{t) + ^Vr(E(r + iq) + y(r- ^q)) + dt' D{t - t')r{t') 

(67) 

and 

mq{t) + 2Vq{V{r+^q) + V{r-^q)) + ^ dt' D{t' - t)q{t') = 0. (68) 
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4.3 The Abraham-Lorentz Equation 

The real part of the equation of motion (67), which is obtained by setting the 
right-hand side equal to zero, yields the famous Abraham-Lorentz equation 
for the electron [17]. It describes the radiation damping through the damping 
kernel D{t — t') [15]. To see this we write the real part of Eq. (67) as 

d /■* 

mr{t) + — dt'D{t - (69) 

di Jo 

where Eext(l) denotes an external force derived from the potential V. The 
damping kernel can be written (see Eqs. (50) and (52)) 

d 

D{t — t')= / dw— ^wsina;(t — t') = — / duj cos oj{t — t') 

Jo oTT Stt dt Jq 

_ d sin d , 



where we have introduced the function 



/«)=”■“. 

To be specific the UV-cutoff 17 is taken to be 

hf} = mc^, 




(70) 

(71) 


which implies that 

mc^ c 

= A^’ 




(72) 


where 

Ac = — 




(73) 


me 

is the Compton wavelength. For an electron we 


have 


Ac ~ 3.8 X 10“^^m and 17 s 


i 0.78 X lO^^s"^ 


(74) 



The term of the equation of motion (69) involving the damping kernel 
can be written as follows. 






- [ dt'/(t - t')r{t') + /(O)f(t) - f{t)r{0). 
. Jo 



r{t') (75) 



d 

dt 



For times t such that ilt 1, i.e. 1 10 ^^s, we may replace 
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and approximate f{t) « 0, while Eq. (70) yields /(O) = ^2. Thus we obtain, 

ct „2 d r 



which finally leads to the equation of motion, 

e^f2\ ^ . _ , , 

+ = -P^ext(t), 



Stt^ ) 



Gtt 



(77) 



(78) 



where v = f is the velocity. This is the famous Abraham-Lorentz equation 
[17]. The term proportional to the third derivative of r(t) describes the damp- 
ing of the electron motion through the emitted radiation. This term does not 
depend on the cutoff frequency, while the cutoff-dependent term yields a 
renormalization of the electron mass. 



niii = m + Am = m + . 



(79) 



It is important to note that the electro-magnetic mass Am diverges linearly 
with the cutoff. The equation of motion (78) can be obtained heuristically by 
means of the Larmor formula for the power radiated by an accelerated charge. 
More rigorously, it has been derived by Abraham and by Lorentz from the 
conservation law for the field momentum, assuming a spherically symmet- 
ric charge distribution and that the momentum is of purely electromagnetic 
origin [17]. 

For the cutoff 17 chosen above we get 



Am 



me^ 

Stt^ 



Stt 



(80) 



and, hence. 



Am 



= —a « 0.0031. 



m 3 tt 

The decomposition (79) of the mass is, however, unphysical, since the electron 
is never observed without its self-field and the associated field momentum. In 
other words, we have to identify the renormalized mass mji with the observed 
physical mass which enables us to write Eq. (78) as 



mR [h(t) - Toh(t)] = Fext(t). 



(81) 



Here, the radiation damping term has been written in terms of a characteristic 
radiation time scale tq given by 

To = ^ « 0.6 X 10“^^s, (82) 

birmR 3 

where denotes the classical electron radius, 

g2 

Te = = aXc ~ 2.8 X 10 ^^m. (83) 

AirmR 
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It is well-known that Eq. (81), being a classical equation of motion for the 
electron, leads to the problem of exponentially increasing runaway solutions. 
Namely, for i^ext = 0 we have 



V — tqV = 0. (84) 

In addition to the trivial solution of a constant velocity, v = const, one also 
finds the solution 

i,{t) = ■w(0)exp(t/ro), 

describing an exponential growth of the acceleration for 'h(O) yf 0. In order 
to exclude these solutions one imposes the boundary condition 

i}{t) — >• 0 for t — oo, 

if -Fext also vanishes in this limit. This boundary condition can be imple- 
mented by rewriting Eq. (81) as an integro-differential equation 

nOO 

mRr{t) = / dsexp(-s)Fext(i + Tos). (85) 

Jo 

On differentiating Eq. (85) with respect to time, it is easily verified that one 
is led back to Eq. (81). However, for Fext = 0 it follows immediately from 
Eq. (85) that V = const, such that runaway solutions are excluded. 

On the other hand, Eq. (85) shows that the acceleration depends upon 
the future value of the force. Hence, the electron reacts to signals lying a time 
of order tq in the future, which is the phenomenon of pre-acceleration. This 
phenomenon should, however, not be taken too seriously, since the description 
is only classical. The time scale Tq corresponds to a length scale Vg which 
is smaller than the Compton wavelength Ac by a factor of a, such that a 
quantum mechanical treatment of the problem is required. 



4.4 Construction of the Decoherence Function 

In this subsection we derive the explicit form of the propagator function (63) 
for the reduced electron density matrix in the case of quadratic potentials, 

V{x) = ^niRLolx'^ . (86) 

Our aim is to introduce and to determine the decoherence function which 
provides a quantitative measure for the degree of decoherence. On using 



V{r + q/2) + V{r- q/2) = mRUjlr"^ + mRiJlq^j^ 

-V (r + q/2) + V(r- q/2) = -mRU^lr ■ q. 



(87) 
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the classical equations of motion take the form 



'mu 

rriR 



r{t)-\-u;l / d5exp(-s)r(t + tqs) 
^0 
pOO 

/ dsexp(-5)g(t - tqs) 
Jo 






= 0 . 



(89) 



Note, that Eq. (89) is the backward equation of the real part of Eq. (88). More 
precisely, if q{t) solves Eq. (89), then r{t) = q{tf — t) is a solution of (88) 
with the right-hand side set equal to zero. 

The above equations of motion lead to the following renormalized action 
functional 



A[r,q] = 



dt rriR 
rtf 



pOO 

r{t)q{t) - Lolq{t) / ds exp(-s)r(t -|- tqs) 

Jo 

dt / dt' Di{t — t')q{t)q{t'). 

Jo 



(90) 



In the following we shall use this renormalised action functional instead of 
the action given in Eq. (65). By variation with respect to q(t) we immediately 
obtain Eq. (88), whereas the variation with respect to r{t) yields: 



dt mu 



pOO 

q{t)Sr{t) + ojI / dsexp(— s)q(t)(5r(t -I- tqs) 
Jo 



= 0 , 



which implies 
rtf 



pf 

dtq{t)5r{t) + / ds i dtexp(— s)q(t)i5r(t -I- tqs) 

Jo Jo 

rtf poo ptf-\-TQS 

/ dtq{t)6r{t) + ujq / ds dtexp{—s)q{t — Tos)Sr{t) 

Jo Jo J Tns 



= 0 . 



(91) 



In the last time integral we may extend the integration over the time interval 
from 0 to t/. This is legitimate since tq is the radiation time scale: By setting 
this variation of the action equal to zero we thus neglect times of the order 
of the pre-acceleration time, which directly leads to the equation of motion 
(89). 

Since the action functional is quadratic the propagator function can be 
determined exactly by evaluating the action along the classical solution and 
by taking into account Gaussian fluctuations around the classical paths. We 
therefore assume that r{t) and q{t) are solutions of the classical equations 
of motion (88) and (89) with boundary conditions (64). The effective action 
along these solutions may be written as 



Aci[r, q] = - riqi] 
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rtf 



dtmRq{t) 



roo 

r{t)+Lol / dsexp(-s)r(i + Tos) 

Jo 



1 

'4 



dt j dt' Di{t — t')q{t)q{t'), 
3 Jo 



(92) 



or, equivalently, 



Ai[r-,q] = ^ ^ dtq{t)^J dt' Di{t - t’)q{t’) 

• pt ^ pt ^ 

+ ^ dt dt'Di{t - t')q{t)q{t'). (93) 

4 Jo Jo 

Eq. (88) shows that the solution r{t) is, in general, complex due to the cou- 
pling to q{t) via the noise kernel Di{t — t'). Consider the decomposition of 
r(t) into real and imaginary part, 

r{t) = r^^^t) +ir^‘^'^{t), (94) 



where is a solution of the real part of Eq. (88), while solves its 
imaginary part. 



rriR 



d2) 



it) 



-UJn 



ds exp(— s)r {t + tqs) 






(95) 

We now demonstrate that, in order to determine the action along the 
classical paths, it suffices to find the homogeneous solution and to insert 
it in the action functional [14]. In other words we have 



■Ac\[r^^\q] = Ac\[r,q], 



(96) 



where 



Aci[r^'^\q] = mR[ry'qf - r^'q,] 






^(1). 



+ - dt dt'Di{t-t')q{t)q{t'). 
^ Jo Jo 

To proof this statement we first deduce from Eq. (95) that 



(97) 



rtf rtf 



dt q{t) 



dt 



dt'Zli(t — t')q{t') 



pt f pOO 

= —iniR / dt q{t) r*'^^(t) -h Wg / dsexp(— s)r*^^^(t -I- rgs) 

Jo L Jo 

r-(2) -(2) 1 

= -miR[r) >qf - r\ >qj\ 

ptf r 1*00 

—iniR / dt r'^"^\t)q(t) + uJq / dsexp(— s)r^^^(t -I- Tgs)<7(t) 

Jo . Jo 



.( 98 ) 
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The term within the square brackets is seen to vanish if one employs Eq. (89) 
and the same arguments that were used to derive the equation of motion 
from the variation (91) of the action functional. Furthermore, we made use 
of r(^)(0) = = 0 which means that the real part of the solution 

satisfies the given boundary conditions. Hence we find 



from which we finally obtain with the help of (93), 



■J2), 



(99) 



Aci[r, q] = mR[rfqf - r^q,] - imR[rf'’qf - 

+ |/ dt J dt'Di{t-t')q{t)q{t') 






(100) 



This completes the proof of the above statement. 

Summarizing, the procedure to determine the propagator function for the 
electron can now be given as follows. One first solves the equations of motion 

pOO 

/ dsexp{—s)r{t-\-Tos)=0, ( 101 ) 

Jo 

poo 

q{t)+uJo / dsexp(-s)q(t - tqs) = 0, (102) 

Jo 

together with the boundary conditions (64). With the help of these solutions 
one then evaluates the classical action, 

■J^c\[r,q] = mR[rfqf - riqi] + ^ f dt [ dt' Di{t - t')q{t)q{t'), (103) 

4 Jo Jo 

which immediately yields the propagator function 



J{'i'f,qf,tf,r,,qi) = iVexp {iAi[r, q]} 

= iVexp{imfl(r/q/ - Viq,) + r{qf,qi,tf)} ■ (104) 

Here, is a normalization factor which is determined from the normalization 
condition 

J d^rfJ{rf,qj = 0,tf;ri,q^) = 6{qi). (105) 

The function r{qf,qi,tf) introduced in Eq. (104) will be referred to as the 
decoherence function. It is given in terms of the noise kernel Di(t — t') as 

r(Qf^QiAf) = -^ J dt J dt'Di{t-t')q{t)q{t'). 



(106) 
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Explicitly we find with the help of Eq. (51), 
r*/ 



r = -1 /■' 

4./„ 



Plf POO 

dt / dt' / du>J{uj)coth{Pu!/2)cosuj{t — t')q{t)q{t'). (107) 
Jo Jo 



The double time-integral can be written as 

f^f f^f f^f 

Re / dt dt' eyip[iuj{t — t')]q{t)q{t') = / dteyip{iujt)q{t) 

Jo Jo Jo 

Hence, the decoherence function takes the form 

1 

r{qf,qi,tf) = -- J dwJ(w)coth(/3w/2) |g(w)|^ 
where we have introduced 

rtf 

Q{uj) = / dt exp{\ujt)q{t) . 



(108) 



(109) 



( 110 ) 



It can be seen from the above expressions that T is a non-positive func- 
tion. The decoherence function will be demonstrated below to describe the 
reduction of electron coherence through the influence of the radiation held. 



5 Decoherence Through the Emission of 
Bremsstrahlung 

As an example we shall investigate in this section the most simple case, 
namely that of a free electron coupled to the radiation held. This case is 
of particular interest since it allows an exact analytical determination of 
the decoherence function and already yields a clear physical picture for the 
decoherence mechanism. Having determined the decoherence function, we 
proceed with an investigation of its influence on the propagation of electronic 
wave packets. 

5.1 Determination of the Decoherence Function 

We set Wo = 0 to describe the free electron. The equations of motion (101) 
and (102) with the boundary conditions (64) can easily be solved to yield 

r{t) = ri+ q{t) = qi+ ^\ (Ill) 

tf 

Making use of Eq. (104) and determining the normalization factor from 
Eq. (105) we thus get the propagator function, 

J{rf,qf,tf;r„q,) = 
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As must have been expected J is invariant under space translations since it 
depends only on the difference rf — ri. Furthermore, one easily recognizes 
that the contribution 

G{rf-r,,qf-q„tf)=(^^^'^ exp |i^(r/ - - qi)| (113) 

is simply the propagator function for the density matrix of a free electron 
with mass mp> for a vanishing coupling to the radiation field. We can thus 
write the electron density matrix as follows, 



Pm{rf,qf,tf) = J dS J d^qiG{rf -Ti,qf - qi,tf) 

X exp{^{qf,q^,tf)}pm{r^,q^,0), (114) 



which exhibits that the decoherence function F describes the influence of the 
radiation field on the electron motion. 

We proceed with an explicit calculation of the decoherence function. It 
follows from Eqs. (110) and (111) that 



Ff 



Q{uj) = / dtexp(ia;t) 

Jo tf 



qf - qi _ exp(iwt/) - 1 



where 



- <li)- 



Therefore, the decoherence function is found to be 



r = - 









, 1 — cos Lot f , , ^ 
dw coth(/3w/2). 



LO 



(115) 

(116) 

(117) 



where we have used expression (52) for the spectral density J{to). The deco- 
herence function may be decomposed into a vacuum contribution Tvac and a 
thermal contribution Fth, 

T = rvac+/^th, (118) 

where ^ ^ ^ 

1 — cos Lot f 



vac g^2 / 



dw- 



10 



LO 



and 



^ , 1 — cos Lot f , , 

rth = -^irz^ I dw ^ ^ [coth(/3w/2) - 1] . 



67T^ 



'0 



LO 



(119) 



(120) 



The frequency integral appearing in the vacuum contribution can be eval- 
uated in the following way. Substituting x = Lot f we get 



dw 



1 — cos Lot f 



LO 



, 1-cosa: , ^ ^ 

dx =lnf2tf + G + O 

X 



lo 



( [2tf ) 



( 121 ) 
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where C « 0.577 is Euler’s constant [18]. For 121/ ^ 1 we obtain asymptoti- 
cally 

= ( 122 ) 

To determine the thermal contribution Tth we first write Eq. (120) as 
follows, 



Tth = - 



g2y,2 rtf rn 



rf 

/ dt dw [coth(/3(^/2) — 1] sinwl = — (123) 
do do Stt 



Introducing the integration variable x = fdoj we can cast the double integral 
/ into the form 

rtf rl30 

I= — I dt dx [coth(a;/2) — 1] sin (lx//3) . 
do 



Here, we have (3fl = hQ/ksT and, using the cutoff hf2 = mc^, we get 



For temperatures T obeying 



I3f2 = 



ksT <C mc^ 



the upper limit of the x-integral may be shifted from (3f2 to oo. Condition 
(124) states that 

rnksT ^ ’ 

which means that the thermal wavelength Ath = h/ \/2mkBT is much larger 
than the Compton wavelength, 

Ath > Ac. (125) 



Thermal and Compton wavelength are of equal size at a temperature of 
about 10® Kelvin. Condition (124) therefore means that T <C 10® K. Under 
this condition we now obtain 



1 rf 



pOO 

dt / dx [coth(x/2) — 1] sin (tx/P) 
do 



1 r*f 

pJo 



, r , { P 

dl TTCOth — 

\P ) t 



sinh {irtf/P) 
TTtf/P 
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where we have employed the formula 



/.oo ^ 

/ dx [coth(x/2) — 1] siuTX = 7rcoth(7rr) . 

Jo T 



(127) 



The quantity 

tb = - = « 2.4 • 10-12 g/T[K] (128) 

7T TTfcs J 

represents the correlation time of the thermal radiation field. Putting these 
results together we get the following expression for the thermal contribution 
to the decoherence function, 



Tth 



/ sinh(t//TB) \ (g/ - g,)^ 
6^2 1, tf/TB ) t) 



(129) 



Adding this expression to the vacuum contribution (122) and introducing 
a = /Att he and further factors of c, we can finally write the expression for 
the decoherence function as 



, 2of 



pn f2tf + ln 
Alternatively, we may write 

r{qf,qt,tf) = - 



sinh(f/ / tb) \ 



(g/ - q^? 



J\ {ctfY 

jqf - 
2A(t/)2 ’ 



where the quantity L{tf) defined by 






ho Alt f + In 



sinh(t//TB) \ 

tf/T-B J\ 






(130) 



(131) 



(132) 



may be interpreted as a time-dependent coherence length. 

The vacuum contribution Tvac to the decoherence function (130) appar- 
ently diverges with the logarithm of the cutoff 17. This is, however, an artificial 
divergence which can be seen as follows. The decoherence function is defined 
in terms of the Fourier transform Q{uj) of q{t), see Eqs. (109) and (110). 
Evaluating Q{lo) as in Eq. (115) we assume that the velocity is zero prior to 
the initial time t = 0, that it suddenly jumps to the value given by Eq. (116), 
and that it again jumps to zero at time tf. This implies a force having the 
shape of two (5-function pulses around t = 0 and t = tf. Such a force acts 
over two infinitely small time intervals and leads to sharp edges in the clas- 
sical path. More realistically one has to consider a finite time scale Tp for the 
action of the force which must be still large compared to the radiation time 
scale tq. We may interpret the time scale Tp as a preparation time since it 
represents the time required to prepare the initial state of a moving electron. 
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A natural, physical cutoff frequency of the order 17 ~ 1/rp is thus introduced 
by the preparation time scale Tp and we may set 

Qtf = ^ (133) 

Tp 

in the following. It should be noted that the weak logarithmic dependence 
on 17 shows that the precise value of the preparation time scale Tp is rather 
irrelevant. The important point is that the preparation time introduces a 
new time scale which removes the dependence on the cutoff. The vacuum 
decoherence function can thus be written, 



Tv, 



(g/ ~ 

37T \tp) {ctfY ’ 



(134) 



showing that it vanishes for large times essentially as 

The thermal contribution Tth is determined by the thermal correlation 
time tb- For T — > 0 we have tb — > oo, and this contribution vanishes. For 
large times tf ^ tb the thermal decoherence function may be approximated 

by 



2a tf {qf - q^f 
3tt Tb (c</)2 



(135) 



which shows that Tth vanishes as Thus, for short times the vacuum 
contribution dominates, whereas the thermal contribution is dominant for 
large times. Both contributions Ihac and Tth are plotted separately in Fig. 1 
which clearly shows the crossover between the two regions of time. 

Eq. (132) implies that the vacuum coherence length is roughly of the order 



T(ty)vac ^ c * tj. 



(136) 



To see this let us assume a typical preparation time scale of the order Tp ~ 
10“^^s. If we take tf to be of the order of Is we find that ln(ty/rp) ^ 48. 
In the rather extreme case tf ^ 10^^, which is of the order of the age of the 
universe, we get ln(t//rp) ~ 87. On using 37r/4a « 322 and Eq. (132) for 
T = 0 one is led to the estimate (136). 



5.2 Wave Packet Propagation 



Having obtained an expression for the decoherence function F we now proceed 
with a detailed discussion of its physical significance. For this purpose it will 
be helpful to investigate first how F affects the time-evolution of an electronic 
wave packet. We consider the initial wave function at time t = 0, 






(^) ' 



{x- 

4ct 



a) 



2 



iko{x - a) , 



(137) 
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Fig. 1. The vacuum contribution TTac and the thermal contribution /th of the de- 
coherence function F (Eq. (130)). For a fixed value |q/ — Qi\ = 0.1 • ctb, the two 
contributions are plotted against the time t/ which is measured in units of the ther- 
mal correlation time tb- The temperature was chosen to be T = IK. One observes 
the decrease of both contributions for increasing time, demonstrating the vanishing 
of decoherence effects for long times. The thermal contribution Tth vanishes as , 
while the vacuum contribution TTac decays essentially as leading to a crossover 
between two regimes dominated by the vacuum and by the thermal contribution, 
respectively. 



describing a Gaussian wave packet centered at a; = a with width (Jq. With 
the help of Eqs. (113), (114) and (131) we get the position space probability 
density at the final time tf, 



Pm{rf,tf) = Pm{rf,qf = 0,t/) 

3 



(138) 



= / d^ri / d^(7* 



/ rriR y 

\2TTtf) 



exp 



iTOfl , ^ Q 

-{Tf - r,)q, - 



tf 



2L(t;)2j 



xV’o(ri + ^g*)V’S(n - ^Qt)- 



The Gaussian integrals may easily be evaluated with the result, 



Pmi'f'fi tf) — 



2TTa{tfY^ 



3/2 



exp 



jrf - b)^ 

2cr(t/)2 



(139) 
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where 



and 



b = a — 



kptf 

rriR 



a{tfr^a^,+ 












(140) 



(141) 



This shows that the wave packet propagates very much like that of a free 
Schrodinger particle with physical mass mR. The centre b of the proba- 
bility density moves with velocity —kp/niR, while its spreading, given by 
Eq. (141), is similar to the spreading which is obtained from the 

free Schrodinger equation. 



0"(^/)free ~ ^0 + 






(142) 



If we write 



cr(t/)^ — <7o + 



y-2 



4m|jcrj 



1-h 



my 



(143) 



we observe that the decoherence function affects the probability density only 
though the width <r{tf) and leads to an increase of the spreading. In view 
of the estimate (136) the correction term in Eq. (143) is, however, small for 
times satisfying 

L(t/) ^ c • » (To. (144) 



This means that the influence of the radiation field can safely be neglected 
for times which are large compared to the time it takes a light signal to travel 
the width of the wave packet. 




Fig. 2. Sketch of the interference experiment used to determine the decoherence 
factor. Two Gaussian wave packets with initial separation 2a approach each other 
with opposite velocities of equal magnitude v = ko/mu. 



Let us now study the evolution of a superposition of two Gaussian wave 
packets separated by a distance 2a. This case has been studied already by 
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Barone and Caldeira [15] who find, however, a different result. We assume 
that the packets have equal widths uo and that they are centered initially at 
X = ±a = ±(a, 0, 0). The packets are supposed to approach each other with 
the speed v = k^jraB. > 0 (see Fig. 2). For simplicity the motion is assumed 
to occur along the a;-axis. Thus we have the initial state 



ipoix) = Ai 



A 9 



27T(Tq^ 



27rcr2y 



3/4 


(a:-a)2 ' 

2 ^kg[x a) 


j exp 


, 3/4 
j exp 


(x + a)^ , , 

2 “1“ ikq[X + a) 

4^0 



(145) 



where kg = (kg, 0,0) and Ai, A 2 are complex amplitudes. Our aim is to 
determine the interference pattern that arises in the moment of collision of 
the two packets at a? = 0. Using again Eqs. (113), (114) and (131) and doing 
the Gaussian integrals we find 



Pm{i' f , tf) 



V27ra(t/)2j 

X {jxlip + \ A2\^ + 2ReAiA*2e-xp[tp{r f)]) . 



(146) 



We recognize a Gaussian envelope centered at ry = 0 with width cr(tf), an 
incoherent sum j^i p + |xl 2 p, and an interference term proportional to AiA^. 
The interference term involves a complex phase given by 

2a? 

= -2ikgrf{l-e)~ (147) 



The term —2\kgrf describes the usual interference pattern as it occurs for a 
free Schrodinger particle, while the contribution 2ikgrfS leads to a modifica- 
tion of the period of the pattern. The final time tf = amn/kg is the collision 
time and 



a kg 

tf m,R 



(148) 



is the speed of the wave packets. The factor e is given by 



'm-RL{tf)^a{tf)‘^ 4cr2 tj j 



(149) 



Obviously we always have 0 < £ < 1. Furthermore, for the situation consid- 
ered in Eq. (144) we have e <C 1. Thus, we get 



2a^ 



= -2ifeory 



L{tfr 



(150) 
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The last expression clearly reveals that the real part of the phase ‘fir f) 
describes decoherence, namely a reduction of the interference contrast de- 
scribed by the factor 



D 





r (2ar 1 




exp 


[ 2L(ty)2j 


= exp 



distance^ 

2 (coherence length)^ ’ 



(151) 



which multiplies the interference term. As was to be expected from the general 
formula for F, the decoherence factor D is determined by the ratio of the 
distance of the two wave packets to the coherence length. 

Alternatively, we can write the decoherence factor in terms of the velocity 
(148) of the wave packets. In the vacuum case we then get 



^vac — exp 



37T \TpJ 




(152) 



This clearly demonstrates that it is the motion of the wave packets which 
is responsible for the reduction of of the interference contrast: If one sets 
into relative motion the two components of the superposition in order to 
check locally their capability to interfere, a decoherence effect is caused by 
the creation of a radiation field. As can be seen from Eq. (117) the spectrum 
of the radiation field emitted through the moving charge is proportional to 
1 /uj which is a typical signature for the emission of bremsstrahlung. Thus 
we observe that the physical origin for the loss of coherence described by the 
decoherence function is the creation of bremsstrahlung. 

It is important to recognize that the frequency integral of Eq. (117) con- 
verges for w — >■ 0, see Eq. (121). The decoherence function F is thus infrared 
convergent which is obviously due to the fact that we consider here a pro- 
cess on a finite time scale tf. This means that we have a natural infrared 
cutoff of the order of l7min ~ l/^/> addition to the natural ultraviolet 
cutoff ~ 1/A> introduced earlier. The important conclusion is that the 
decoherence function is therefore infrared as well as ultraviolet convergent. 

It might be instructive, finally, to compare our results with the corre- 
sponding expressions which are derived from the famous Caldeira-Leggett 
master equation in the high-temperature limit (see, e.g. [3]). From the latter 
one finds the following expression for the coherence length 



^(^/)cL — 



^th 



(153) 



where 7 is the relaxation rate. This is to be compared with the expressions 
(132) for the coherence length. For large temperatures we have the following 
dominant time and temperature dependence, 

L{tff ~ Y ^(^/)cL ~ (154) 
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Hence, while both expressions for the coherence length are proportional to 
the inverse temperature, the time dependence is completely different. Namely, 
for tf — > oo we have 

L{tff — ^ oo and L{tf)cL — ^ 0, (155) 

and, therefore, complete coherence in the case of bremsstrahlung and total 
destruction of coherence in the Caldeira-Leggett case. 



6 The Harmonically Bound Electron in the Radiation 
Field 



As a further illustration let us investigate briefly the case of an electron in the 
radiation field moving in a harmonic external potential. Another approach 
to this problem may be found in [19], where the authors arrive, however, at 
the conclusion that there is no decoherence effect in the vacuum case. 

We take wq > 0 and solve the equation of motion (101) with the help of 
the ansatz 

r(t) = roexp(zt), (156) 

where, for simplicity, we consider the motion to be one-dimensional. Substi- 
tuting this ansatz into (101) one is led to a cubic equation for z, 

— tqz^ + Wq = Q. (157) 



For vanishing coupling to the radiation field (tq = 0) the solutions are lo- 
cated at z± = ±iwo) describing the free motion of a harmonic oscillator with 
frequency loq. 

For To > 0 the cubic equation has three roots, one is real and the other 
two are complex conjugated to each other. The real root corresponds to the 
runaway solution and must be discarded. Let us assume that the period of 
the oscillator is large compared to the radiation time. 

To < — . (158) 

Wo 



Because of To ^ 10“^^s this assumption is well satisfied even in the regime 
of optical frequencies. We may thus determine the complex roots to lowest 
order in loqTq, 



z± = ±iwo - ^tqujq. 



(159) 



The purely imaginary roots ±iwo of the undisturbed harmonic oscillator are 
thus shifted into the negative half plane under the influence of the radiation 
field. The negative real part describes the radiative damping. In fact, we see 
that r{t) decays as exp(— 7 ^/ 2 ), where 



7 = toWq = -a 



fkvQ 

mjiC^ 



(160) 
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is the damping constant for radiation damping [17]. In the following we con- 
sider times tf of the order of magnitude of one period ojotf ~ 1 - Because of 
jtf = (tooTo) {ujot f) we then have 7^/ ~ tqWo 1. In this case the damping 
can be neglected and we may use the free solution in order to determine the 
decoherence function. 

Let us consider again the case of a superposition of two Gaussian wave 
packets in the harmonic potential. The packets are initially separated by a 
distance 2a and approach each other with opposite velocities of equal mag- 
nitude such that they collide after a quarter of a period, tf = 7t/2wo. The 
corresponding free solution q(t) is therefore given by 



= Qi cos Wot + <7/ sin ujot. 



(161) 



To describe the situation we have in mind we take qi = 2a (initial separation 
of the wave packets) and 9 / = 0 (to get the probability density). Hence, we 
have 

q{t) = —2auJo sin uot, (162) 

and we evaluate the Fourier transform. 



Q{uj) 




dt exp{iujt)q{t) 



aojQ 



exp(i[a; -I- oJo]tf) ~ 1 
CO + LOo 



exp(i[u; — coo]tf) — 1 
10 — IOq 



This yields the decoherence function 



r = r{qf = 0,9*,f/) 

e^(awo)^ r 1 — cos(w -I- Wo)G 



1 — COs(w — iO[j)tf 
{to - Wo)^ 



coth 



(163) 




We discuss the case of zero temperature. The frequency integral in Eq. 
(163) then approaches asymptotically the value 2 In Qtf which leads to the 
following expression for the decoherence factor, 



D 



vac 



exp Tvac = exp 



37T \TpJ 




(164) 



The interesting point to note here is that this equation is the same as Eq. (152) 
for the free electron, with the only difference that the square (u/c)^ of the 
velocity, which was constant in the previous case, must now be replaced with 
its time averaged value ((u/c)^). 



7 Destruction of Coherence of Many-Particle States 

For a single electron the vacuum decoherence factor (152) turns out to be very 
close to 1 , as can be illustrated by means of the following numerical example. 
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We take Tp to be of the order of 10“^^s and t/ of the order of Is. Using a 
velocity v which is already as large as 1/10 of the speed of light, one finds 
that Uvac ^ 10“^, corresponding to a reduction of the interference contrast of 
about 1%. This demonstrates that the electromagnetic field vacuum is quite 
ineffective in destroying the coherence of single electrons. 

For a superposition of many-particle states the above picture can lead, 
however, to a dramatic increase of the decoherence effect. Consider the su- 
perposition 

\^) = IV'i) + IV'2) (165) 

of two well-localized, spatially separated fV-particle states If/i) and \ip 2 )- We 
have seen that decoherence results from the imaginary part of the influence 
phase functional <P[Jc, Ja], that is from the last term on the right-hand side 
of Eq. (39) involving the anti-commutator function Di{x — of the elec- 
tromagnetic field. Thus, it is the functional 

r[Jc] = -^j^ d'^a;^ d‘^x'Di{x-x')p,^J^{x)J';!^{x'), (166) 

which is responsible for decoherence. This shows that the decoherence func- 
tion for fV-electron states scales with the square of the particle number. 
Thus we conclude that for the case of the superposition (165) the decoherence 
function must be multiplied by a factor of that is the decoherence factor 
for Wparticle states takes the form, 



D 



N 



exp 









37T ' 




\c) 



(167) 



This scaling with the particle number obviously leads to a dramatic increase 
of decoherence for the superposition of Wparticle states. To give an example 
we take N = 6-10^^, corresponding to 1 mol, and ask for the maximal velocity 
V leading to a 1% suppression of interference. With the help of (167) we find 
that V ~ 10“^®m/s. This means that, in order to perform an interference 
experiment with 1 mol electrons with only 1% decoherence, a velocity of at 
most 10“^®m/s may be used. For a distance of Im this implies, for example, 
that the experiment would take 3 x 10® years! 



8 Conclusions 

In this paper the equations governing a basic decoherence mechanism oc- 
curring in QED have been developed, namely the suppression of coherence 
through the emission of bremsstrahlung. The latter is created whenever two 
spatially separated wave packets of a coherent superposition are moved to one 
place, which is indispensable if one intends to check locally their capability to 
interfere. We have seen that the decoherence effect through the electromag- 
netic radiation field is extremely small for single, non-relativistic electrons. 
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The decoherence mechanism is thus very ineffective on the Compton length 
scale. An important conclusion is that decoherence does not lead to a local- 
ization of the particle on arbitrarily small length scales and that no problems 
with associated UV-divergences arise here. 

The decoherence mechanism through bremsstrahlung exhibits a highly 
non-Markovian character. As a result the usual picture of decoherence as a 
decay of the off-diagonals in the reduced density matrix does not apply. In 
fact, consider a superposition of two wave packets with zero velocity. The 
expression (131) for the decoherence function together with the estimate 
L{tf)vac ~ c-tf for the vacuum coherence length T(t/)vac show that de- 
coherence effects are negligible for times tf which are large in comparison to 
the time it takes light to travel the distance between the wave packets. The 
off-diagonal terms of the reduced density matrix for the electron do therefore 
not decay at all, which shows the profound difference between the decoher- 
ence mechanism through bremsstrahlung and other decoherence mechanisms 
(see, e. g. [20]). 

A result of particular interest from a fundamental point of view is that 
coherence can already be destroyed by the presence of the electromagnetic 
field vacuum if superpositions of many-particle states are considered. An im- 
portant conclusion which can be drawn from this picture of decoherence in 
QED refers to various alternative approaches to decoherence and the closely 
related measurement problem of quantum mechanics: In recent years sev- 
eral attempts have been made to modify the Schrodinger equation by the 
addition of stochastic terms with the aim to explain the non-existence of 
macroscopic superpositions through some kind of macrorealism. Namely, the 
random terms in the Schrodinger equation lead to a spontaneous destruc- 
tion of superpositions in such a way that macroscopic objects are practically 
always in definite localized states. Such approaches obviously require the in- 
troduction of previously unknown physical constants. In the stochastic theory 
of Ghirardi, Pearle and Rimini [21], for example, a single particle microscopic 
jump rate of about has to be introduced such that decoherence is 

extremely weak for single particles but acts sufficiently strong for many par- 
ticle assemblies. It is interesting to observe that the decoherence effect caused 
by the presence of the quantum field vacuum yields a similar time scale in a 
completely natural way without the introduction of new physical parameters. 
Thus, QED indeed provides a consistent picture of decoherence and it seems 
unnecessary to propose new ad hoc theories for this purpose. 

It must be emphasized that the above picture of decoherence in QED 
has been derived from the well-established basic postulates of quantum me- 
chanics and quantum field theory. It therefore does not, of course, constitute 
a logical disprove of alternative approaches. However, it does represent an 
example for a basic decoherence mechanism in a microscopic quantum field 
theory. In particular, it provides a unified explanation of decoherence which 
does not suffer from problems with renormalization (as they occur, e.g. in 
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alternative theories [22]) and which does not exclude a priori the existence of 
macroscopic quantum coherence. Only under certain well-defined conditions 
regarding time scales, relative velocities and the structure of the state vector, 
it is true that decoherence becomes important. Thus, decoherence is traced 
back to a dynamical effect and not to a modification of the basic principles 
of quantum mechanics. 

In this paper we have discussed in detail only the non-relativistic approxi- 
mation of the reduced electron dynamics. For a treatment of the full relativis- 
tic theory, including a Lorentz invariant characterization of the decoherence 
induced by the vacuum field, one can start from the formal development given 
in section 2. An investigation along these lines could also be of great interest 
for the study of measurement processes in the relativistic domain [23] . 
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Abstract. It is well known that the dynamical mechanism of decoherence may 
cause apparent superselection rules, like that of molecular chirality. These ‘environ- 
ment-induced’ or ‘soft’ superselection rules may be contrasted with ‘hard’ super- 
selection rules, like that of electric charge, whose existence is usually rigorously 
demonstrated by means of certain symmetry principles. We address the question of 
whether this distinction between ‘hard’ and ‘soft’ is well founded and argue that, 
despite first appearance, it might not be. For this we first review in detail some of 
the basic structural properties of the spaces of states and observables in order to 
establish a fairly precise notion of superselection rules. We then discuss two exam- 
ples: 1.) the Bargmann superselection rule for overall mass in ordinary quantum 
mechanics, and 2.) the superselection rule for charge in quantum electrodynamics. 



1 Introduction 

To explain the (apparent) absence of interferences between macroscopically 
interpretable states - like states describing spatially localized objects - is the 
central task for any attempt to resolve the measurement problem. First at- 
tempts in this direction just imposed additional rules, like that of the Copen- 
hagen school, who defined a measurement device as a system whose state- 
space is classical, in the sense that the superposition principle is fully broken: 
superpositions between any two states simply do not exist. In a more modern 
language this may be expressed by saying that any two states of such a sys- 
tem are disjoint, i.e., separated by a superselection rule (see below). Proper 
quantum mechanical systems, which in isolation do obey the superposition 
principle, can then inherit superselection rules when coupled to such classical 
measurement devices. 

Whereas there can be no doubt that the notion of classicality, as we 
understand it here, is mathematically appropriately encoded in the notions of 
disjointness and supers election rules, there still remains the physical question 
how these structures come to be imposed. In particular, if one believes that 
fundamentally all matter is described by some quantum theory, there is no 
room for an independent classical world. Classicality should be a feature that 
is emerging in accordance with, and not in violation of, the basic rules of 
quantum mechanics. This is the initial credo of those who believe in the 
program of decoherence [16], which aims to explain classicality by means of 
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taking into account dynamical interactions with ambient systems, like the 
ubiquitous natural environment of the situation in question. (Note: It is not 
claimed to resolve the full measurement problem.) This leads to the notion 
of ‘environment-induced superselection rules’ [40]. 

Fundamental to the concept of dynamical decoherence is the notion of 
‘delocalization’ [24]. The intuitive idea behind this is that through some dy- 
namical process certain state characteristics (‘phase relations’), which were 
locally accessible at one time, cease to be locally accessible in the course 
of the dynamical evolution. Hence locally certain superpositions cannot be 
verified anymore and an apparent obstruction to the superposition principle 
results. Such mechanisms are considered responsible for the above mentioned 
environment-induced superselection rule, of which a famous physical exam- 
ple is that of molecular chirality (see e.g. [38] and references therein). It has 
been established in many calculations of realistic situations that such dy- 
namical processes of delocalization can be extremely effective over short time 
scales. But it is also intuitively clear that, mathematically speaking, it will 
never be strict in any finite time. Hence one will have to deal with notions of 
approximate- respectively asymptotic (for t — >■ oo) superselection rules and 
disjointness of states [31,26], which needs some mathematical care. 

Since for finite times such dynamical superselection rules are only approx- 
imately valid, they are sometimes called ‘soft’. In contrast, ‘hard’ superse- 
lection rules are those which are rigorously established mathematical results 
within the kinematical framework of the theory, usually based on symmetry 
principles (see section 3 below), or on first principles of local QFT, like in the 
proof for the superselection rule for electric charge [33]. Such presentations 
seem to suggest that there is no room left for a dynamical origin of ‘hard’ 
superselection rules, and that hence these two notions of superselection rules 
are really distinct. However, we wish to argue that at least some of the ex- 
isting proofs for ‘hard’ superselection rules give a false impression, and that 
quite to the contrary they actually need some dynamical input in order to be 
physically convincing. We will look at the case of Bargmann’s superselection 
rule for total mass in ordinary quantum mechanics (which is clearly more of 
an academic example) and that of charge in QED. The discussion of the latter 
will be heuristic insofar as we will pretend that QED is nothing but quan- 
tum mechanics (in the Schrodinger representation) of the infinite-dimensional 
(constrained) Hamiltonian system given by classical electrodynamics. For a 
brief but general orientation on the subject of superselection rules and the 
relevant references we refer to Wightman’s survey [39]. 

Let us stress again that crucial to the ideas presented here is of course 
that ‘delocalized’ does not at all mean ‘destroyed’, and that hence the loss 
of quantum coherence is only an apparent one. This distinction might be 
considered irrelevant FAPP (for all practical purposes) but it is important 
in attempts to understand apparent losses of quantum coherence within the 
standard dynamical framework of quantum mechanics. 
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As used here, the term ‘local’ usually refers to locality in the (classical) 
configuration space Q of the system, where we think of quantum states in the 
Schrodinger representation, i.e., as L^-functions on Q. Every parametrization 
of Q then defines a partition into ‘degrees of freedom’. Locality in Q is a more 
general concept than locality in ordinary physical space, although the latter 
forms a particular and physically important special case. Moreover, on a 
slightly more abstract level, one realizes that the most general description of 
why decoherence appears to occur is that only a restricted set of so-called 
physical observables are at ones disposal, and that with respect to those the 
relevant ‘phase relations’ apparently fade out of existence. It is sometimes 
convenient to express this by saying that decoherence occurs only with respect 
(or relative) to a ‘choice’ of observables [27]. Clearly this ‘choice’ is not meant 
to be completely free, since it has to be compatible with the dynamical laws 
and the physically realizable couplings (compare [24]). (In this respect the 
situation bears certain similarities to that of ‘relevant’ and ‘irrelevant’ degrees 
of freedom in statistical mechanics.) But to fully control those is a formidable 
task - to put it mildly. In any case it will be necessary to assume some a priori 
characterizations of what mathematical objects correspond to observables, 
and to do this in such a general fashion that one can effectively include 
superselection rules. This will be done in the next section. 

2 Elementary Concepts 

In this section we wish to convey a feeling for some of the concepts underlying 
the notion of superselection rules. We will take some care and time to do 
this, since many misconceptions can (and do!) arise from careless uses of 
these concepts. To gain intuition it is sometimes useful to dispense with some 
technicalities associated with infinite dimensions and continuous spectra and 
just look at finite dimensional situations; we will follow this strategy where 
indicated. We use the following, generally valid notations: TL denotes a Hilbert 
space, B{'H) the algebra of bounded operators on H. The antilinear operation 
of taking the hermitean conjugate is denoted by * (rather than f), which 
makes B{'H) a *-algebra. Given a set {Av} where A G A (= some index set), 
then by {A^}' we denote the ‘commutant’ of {Aa} in B{'H), defined by 

{Aa}' := {B G B{H) I HAa = AaH, VA G A}. (1) 

Note that if the set {Aa} is left invariant under the *-map (in this case we 
call the set ‘self-adjoint’), then {Aa}' is a *-subalgebra of B{'H). Also, the 
definition (1) immediately implies that 

QA'. (2) 

2.1 Superselection Rules 

There are many different ways to give a structural definition of superselection 
rules. Some stress the notion of state others the notion of observable. Often 
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this dichotomy seems to result in very different attitudes towards the fun- 
damental significance of superselection rules. This really seems artificial in a 
quantum mechanical context. In quantum field theory, i.e., if the underlying 
classical system has infinitely many degrees of freedom, the situation seems 
more asymmetric. This is partly due to the mathematical difficulties to de- 
fine the full analog of the Schrodinger representation, i.e., to just construct 
the Hilbert space of states as space over the classical configuration space. 
In this paper we will partly ignore this mathematical difficulty and proceed 
heuristically by assuming that such a Schrodinger representation (of QED) 
exists to some level of rigour. 

In traditional quantum mechanics, which stresses the notion of state, a 
system is fundamentally characterized by a Hilbert space, TL, the vectors 
of which represent (pure) states. We say ‘represent’ because this labeling 
by states through vectors is redundant: non-zero vectors which differ by an 
overall complex number label the same state, so that states can also be labeled 
by rays. We will use VH to denote the space of rays in "H. In many cases of 
interest this Hilbert space is of course just identified with the space of L^- 
functions over the classical configuration space. Now, following the original 
definition given by [35], we say that a superselection rule operates on 
"H, if not all rays represent pure states, but only those which lie entirely in 
certain mutually orthogonal subspaces Hi C H, where 

n = ^m. (3) 

i 

The only rays which correspond to pure states are those in the disjoint union 

\JVH,. (4) 

I 

Since no vector which lies skew to the partition (3) can, by assumption, rep- 
resent a pure state, the superposition principle must be restricted to the Hi- 
Moreover, since observables map pure states to pure states, they must leave 
the Hi invariant and hence all matrix-elements of observables between vectors 
from different sectors vanish. The Hi are called coherent sectors if the observ- 
ables act irreducibly on them, i.e., if no further decomposition is possible; this 
is usually implied if a decomposition (3) is written down. States which lie in 
different coherent sectors are called disjoint. Note that disjointness of states 
is essentially also a statement about observables, since it means orthogonal- 
ity of the original states and the respective states created from those with 
all observables. The existence of disjoint states is the characteristic feature 
of superselection rules. 

From this we see that a partition (3) into coherent sectors implies that the 
set of physical observables is strictly smaller than the set of all self-adjoint 
(w.l.o.g. bounded) operators on H. It can be characterized by saying that 
observables are those self adjoint operators on H which commute with the 
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orthogonal projectors Pi ^ 'Hi. So the Pi are themselves observables and 
generate the center (see below) of the algebra of observables. 

This suggests a ‘dual’, more algebraic way to look at superselection rules, 
which starts with the algebra of observables O. Then superselection rules are 
said to occur if the algebra of observables, O, - which we think of as being 
given by bounded operators on some Hilbert space H^ - has a non-trivial 
center 0°. Recall that 

0‘^ -.= {A&0\AB = BA, \/B G O}. (5) 

Suppose 0° is generated by self-adjoint elements {C^, /x = 1, 2, ..} which have 
simultaneous eigenspaces Hi, then the Hi’s are just the coherent sectors. 
Indeed, as already remarked, matrix elements of operators from O between 
states from different coherent sectors (i.e. differing in the eigenvalue of at 
least one C^) necessarily vanish. Thus if (j>i and 4>2 are two non-zero vectors 
from Hi and Hj with i yf j, their superposition (j) := (pi + p 2 defines a state 
whose density matrix p := (=orthogonal projector onto the ray generated 
by (p) satisfies 



tr(pH) — tr((Aipi -I- A 2 P 2 )^)> VH G O, (6) 

where Ai ,2 = ll</>i, 2 |P/||<('|P and pi ,2 = ^Vi, 2 - This means that p is a non-pure 
state of O, since it can be written as a non-trivial convex combination of other 
density matrices. Hence we come back to the statements expressed by (3) and 
(4). Also note the following: in quantum mechanics the decomposition of a 
non-pure density matrix as a convex combination of pure density matrices - 
the so-called extremal decomposition - is generically not unique, thus pre- 
venting the (ignorance-) interpretation as statistical “mixtures”.^ However, 
for the special density matrices of the form p = \(p){(p\, where \p) G H, the 

^ In Algebraic Quantum Mechanics one associates to each quantum system an 
abstract C*-algebra, C, which is thought of as being the mathematical object that 
fully characterizes the system in isolation, i.e. its intrinsic or ‘ontic’ properties. 
But this is not yet what we call the algebra of observables. This latter algebra 
is not uniquely determined by the former. It is obtained by studying faithful 
representations of C in some Hilbert space H, such that C can be identified with 
some subalgebra of B{H) (the bounded operators on H). This is usually done 
by choosing a reference state (positive linear functional) on C and performing 
the GNS construction. Then C inherits a norm which is used to close C (as 
topological space) in B{H). It is this resulting algebra which corresponds to our 
O. Technically speaking it is a von Neumann algebra which properly contains an 
embedded copy of C. The added observables (those in O — C) do not describe 
intrinsic but contextual properties. For example, it may happen that O has non- 
trivial center whereas C doesn’t. In this case the superselection rules described 
by O are contextual. See [31] for a more extended discussion of this point. 

^ Hence the term ‘mixture’ for a non-pure state is misleading since we cannot tell 
the components and hence have no ensemble interpretation. For this reason we 
will say ‘non-pure state’ rather than ‘mixture’. 
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extremal decomposition is unique and given hy 4> = where 4>i is 

the orthogonal projection of into T-Li, the orthogonal projector onto 
(pi’s ray, and Xi = This is the relevance of superselection rules 

for the measurement problem: to produce unique extremal decompositions 
- and hence statistical ‘mixtures’ in the proper sense of the word - into an 
ensemble of pure states. There is a long list of papers dealing with the mathe- 
matical problem of how superselection sectors can arise dynamically; see e.g. 
[19,30,2,26] and the more general discussions in [27,31]. 



2.2 Dirac’s Requirement 

Dirac was the first who spelled out certain rules concerning the spaces of 
states and observables [6]. He defined the notion of compatible (i.e., simul- 
taneously performable) observations, which mathematically are represented 
by a set of commuting observables, and the notion of a complete set of such 
observables, which is meant to say that there is precisely one state for each 
set of simultaneous “eigenvalues”. Starting from the hypothesis that states 
are faithfully represented by rays, Dirac deduced that a complete set of such 
mutually compatible observables existed. But this only makes sense if all the 
observables in question have purely discrete spectra. 

In the general case one has to proceed differently: We heuristically define 
Dirac’s requirement as the statement, that there exists at least one complete 
set of mutually compatible observables and show how it can be rephrased 
mathematically so that it applies to all cases. In doing this we essentially 
follow Jauch’s exposition [22]. To develop a feeling for what is involved, we 
will first describe some of the consequences of Dirac’s requirement in the most 
simplest case: a finite dimensional Hilbert space. We will use this insight to 
rephrase it in such a way to stay generally valid in infinite dimensions. 



Gaining intuition in finite dimensions. So let % be an n-dimensional 
complex Hilbert-space, then BifH) is the algebra of complex n x n matrices. 
Physical observables are represented by hermitean matrices in BifH), but 
we will explicitly not assume the converse, namely that all hermitean matri- 
ces correspond to physical observables. Rather we assume that the physical 
observables are somehow given to us by some set S of hermitean matrices. 
This set does not form an algebra, since taking products and complex linear 
combinations does not preserve hermiticity. But for mathematical reasons it 
would be convenient to have such an algebraic structure, and just work with 
the algebra O generated by this set, called the algebra of observables. [Note 
the usual abuse of language, since only the hermitean elements in O are ob- 
servables.] But for this replacement of 5 by G to be allowed S must have been 
a set of hermitean matrices which is uniquely determined by O, for otherwise 
we can not reconstruct the set S from O. To grant us this mathematical con- 
venience we assume that S was already maximal, i.e. that S already contains 
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all the hermitean matrices that it generates. But we stress that there seems 
to be no obvious reason why in a particular practical situation the set of 
physically realizable observables should be maximal in this sense. 

We may choose a set {Oi, . . . Om} of hermitean generators of O. Then O 
may be thought of as the set of all complex polynomials in these (generally 
non-commuting) matrices. But note that we need not consider higher powers 
than (n — 1) of each Oi, since each complex n x n matrix O is a zero of its 
own characteristic polynomial po, i.e. satisfies po{0) = 0, by the theorem of 
Cayley-Hamilton. Since this polynomial is of order n, O” can be re-expressed 
by a polynomial in O of order at most (n — 1). For example, the *-algebra 
generated by a single hermitean matrix O can be identified with the set of 
all polynomials of degree at most (n — 1) and whose multiplication law is as 
usual, followed by the procedure of reducing all powers n and higher of O via 
po{0) = 0. 

Now let {^ 1 , • • • , Am} ='■ {Ai} be a complete set of mutually commuting 
observables. It is not difficult to show that there exists an observable A and 
polynomials pi, i = such that Ai = Pi{A) (see [20] for a simple 

proof). This actually means that the algebra generated by {Ai} is just the n- 
dimensional algebra of polynomials of degree at most n — 1 in ^ (see below for 
justification), which we call A. This algebra is abelian, which is equivalently 
expressed by saying that A is contained in its commutant (compare (1)): 



AC A' 



‘A is abelian’ 



( 7 ) 



Now comes the requirement of completeness. In terms of A it is easy to 
see that it is equivalent to the condition that A has a simple spectrum (i.e. 
the eigenvalues are pairwise distinct). This has the following consequence: 
Let B be an observable that commutes with A, then B is also a function 
of A, i.e., Pb{B) = A for some polynomial pb- The proof is simple: We 
simultaneously diagonalize A and B with eigenvalues and j3a, a = 1, • • • , n. 
We wish to find a polynomial of degree n — 1 such that pB{ota) = Pa- Writing 
Pb{x) = -I- • • • -I- oo, this leads to a system of n linear equations 

{a\ := power of Oa) 



n— 1 

ttoOfa = Pa, for a = 1, ... , n, (8) 

6=0 

for the n unknowns (ag, • • • , On-i). Its determinant is of course just the Van- 
dermonde determinant for the n tuple (oi, • • • , a„): 



det{a[]} = J]^(aa - «&), 

a<b 



(9) 



which is non-zero if and only if (=iff) A’s spectrum is simple. This implies 
that every observable that commutes with A is already contained in A. (It 
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follows from this that the algebra generated by {Ai} is equal to, and not 
just a subalgebra of, the algebra generated by {A}, as stated above.) Since a 
*-algebra is generated by its self-adjoint elements (observables), A cannot be 
properly enlarged as abelian *-algebra by adding more commuting generators. 
In other words, A is maximal. Since A' is a *-algebra, this can be equivalently 
expressed by 



A' C A 



‘A is maximal’ 



(10) 



Equations (7) and (10) together are equivalent to Dirac’s condition, which 
can now be stated in the following form, first given by Jauch [22]: the algebra 
of observables O contains a maximal abelian *-subalgebra A C O, i.e.. 



Dirac’s requirement, 1®* version: 3 A^ O satisfying A = A' 



( 11 ) 



This may seem as if Dirac’s requirement could be expressed in purely 
algebraic terms. But this is deceptive, since the very notion of ‘commutant’ 
(compare (1)) makes reference to the Hilbert space % through B{TL). Without 
further qualification the term ‘maximal’ always means maximal in B{Ti.).^ 
This reference to % can be further clarified by yet another equivalent 
statement of Dirac’s requirement. Since A consists of polynomials in the 
observable A, which has a simple spectrum, the following is true: there exists 
a vector \g) & such that for any vector (j> G T~L there exists a polynomial 
such that 

= !<(’)■ ( 12 ) 

Such a vector |g) is called a generating or cyclic vector for ^ in "H. The 
proof is again very simple: let {^i, • • • , </>„} be the pairwise distinct, non-zero 
eigenvectors of A (with any normalization); then choose 



= ( 13 ) 

i=l 



Equation (12) now defines again a system of n linear equations for the n 
coefficients a„_i,---,ao of the polynomial p,f,, whose determinant is again 
the Vandermonde determinant (9) for the n eigenvalues of A. 

Conversely, if A had an eigenvalue, say a\, with eigenspace "Hi of two or 
higher dimensions, then such a cyclic \g) cannot exist. To see this, suppose it 
did, and let \(pj^) £ Hi be orthogonal to the projection of jg) into Hi- Then 
{(j>j-\p{A)g) = 0 for all polynomials p. Thus \4>i) is unreachable, contradicting 
our initial assumption. Hence a simple spectrum of A is equivalent to the 
existence of a cyclic vector. 

® The condition for an abelian A £ O to he maximal in O wonld he A = A C\0. 
Snch abelian subalgebras always exist, in contrast to those A £ O which satisfy 
the stronger condition to be maximal in the ambient algebra B{H). 
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The general case. In infinite dimensions we have to care a little more 
about the topology on the space of observables, since here there are many 
inequivalent ways to generalize the finite dimensional case. The natural choice 
is the so-called ‘weak topology’, which is characterized by declaring that a 
sequence {Ai} of observables converges to the observable A if the sequence 
{(j)\Ai\ip) of complex numbers converges to {(plAl'tp) for all \(j)), If/') G "H. Hence 
one also requires that the algebra of observables is weakly closed (i.e., closed 
in the weak topology). Such a weakly closed *-subalgebra of B{'H) is called 
a W*- or von-Neumann-algebra (we shall use the first name for brevity). 

A crucial and extremely convenient point is, that the weak topology is 
fully encoded in the operation of taking the commutant (see (1)), in the 
following sense: Let {Aa} be any subset of then {Aa}' is automatically 

weakly closed (see [22] p 716 for a simple proof) and hence a lT*-algebra. 
Moreover, the weak closure of a *-algebra A C B{Ti) is just given by A” 
(the commutant of the commutant). Hence we can characterize a IT*-algebra 
purely in terms of commutants: A is W* iff A = A” . 

This allows to easily generalize the notion of ‘algebra generated by observ- 
ables’: Let {Oa} be a set of self-adjoint elements in B{%), then O := {0\}'' 
is called the {W*-) algebra generated by this set. This definition is natural 
since {Oa}" is easily seen to be the smallest IT*-algebra containing {Oa}> for 
if {Oa} C B C O for some IT*-algebra B, then taking the commutant twice 
yields B = O.^ 

Now we see that Dirac’s requirement in the form (11) directly translates to 
the general case if all algebras involved (i.e. A and O) are understood as W*- 
algebras. Now we also know what a ‘complete set of (bounded) commuting 
observables’ is, namely a set {Aa} C B{'H) whose generated lT*-algebra 
A := {Aa}" is maximal abelian: A = A'. This latter condition is again 
equivalent to the existence of a cyclic vector \g) for A, where in infinite 
dimensions the definition of cyclic is that |A|(/)} is dense in (rather than 
equal to) T~L. It is also still true that there is an observable A such that all 
A\ are functions (in an appropriate sense, not just polynomials of course) 
of A [34]. But since A’s spectrum may be (partially) continuous, there is no 
direct interpretation of a ‘simple’ spectrum as in finite dimensions. Rather, 
one now defines simplicity of the spectrum of A by the existence of a cyclic 
vector for A = {A}". 

Now we come to our final reformulation of Dirac’s condition. Namely, 
looking at (11), we may ask whether we could not reformulate the existence 
of such a maximal abelian A purely in terms of the algebra of observables O 
alone. This is indeed possible. We have AQO^ O' C A' = AQO, hence 
O' C O. Since O = O" the last condition is equivalent to saying that O' is 
abelian {O' C O"), or to saying that O' is the center O" of O, since by (1) and 

^ Note: for any M C B{'H) definition (1) immediately yields M C M" and hence 
M' D M'" (by (2)). But also M' C M"' (by replacing M — >• M'); therefore 
M' = M'" for any M C B{H). 
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(5) the center can be written as = O (1 O' . Now, conversely, it was shown 
in [23] that an abelian O' implies the existence of a maximal abelian A Q O. 
Hence we have the following alternative formulation of Dirac’s requirement, 
first spelled out, independently of (11), by Wightman [37], who called it the 
‘hypothesis of commutative superselection rules’: 



Dirac’s requirement, 2'’'^ version: O' is abelian 



(14) 



There are several interesting ways to interpret this condition. From its 
derivation we know that it is equivalent to the existence of a maximal abelian 
A C O. But we can in fact make an apparently stronger statement, which also 
relates to the earlier footnote 3, namely: (14) is equivalent to the condition, 
that any abelian A ^ O that is maximal in O, i.e. satisfies A = A' C\0, 
also maximal in B{TL).^ 



2.3 Dirac’s Condition and Gauge Symmetries 

Another way to understand (14) is via its limitations on gauge-symmetries. 
To see this, we mention that any lF*-algebra is generated by its unitary 
elements. Hence O' is generated by a set {U\} of unitary operators. Each 
U\ commutes with all observables and therefore generates a one-parameter 
group of gauge-transformations. Condition (14) is then equivalent to saying 
that the total gauge group, which is generated by all U\, is abelian. Note also 
that an abelian O' implies that the gauge-algebra, {Ux}" = O', is contained 
in the observables. O' C O" = O, so that O' = O''. From this one can infer 
the following central statement: 



Dirac’s requirement implies that gauge- and 
sectorial structures are fully determined by the 
center O" of the algebra of observables O. 



(15) 



To see in what sense this is true we remark that for lF*-algebras we can 
simultaneously diagonalize all observables in O". That means that we can 
write Ti in an essentially unique way as direct integral over the real line of 
Hilbert spaces 'H(X) using some (Lebesgue-Stieltjes-) measure cr: 

n= / da(A)-H(A). (16) 

JR 

Operators in O respect this decomposition in the sense that each O G O acts 
on V. componentwise via some bounded operator 0{\) on 'H(X). If O G O" 

® Proof: We need to show that {O' abelian) {A = A' C\ O ^ A — A'). ‘=^’: O' 
abelian implies O' C O" — O and AGO implies O' C A' , so that O' G A' C\ O. 
Hence A = A' GO implies O' G A, which implies A' G O" = O, and hence 
A — A' . ‘<i=’: {A = A' G O ^ A — A') is equivalent to A' G O, which implies 
O' G A' = A and hence that O' is abelian. 
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then each 0(A) is a multiple ^(A) G C of the unit operator. Moreover, the 
set of all {0(A)} induced from O for each fixed A acts irreducibly on "H(A).® 
Hence, provided that Dirac’s requirement is satisfied, (16) is the generally 
valid version of (3). The notion of disjointness now acquires an intuitive 
meaning: two states I'f'i) and ltf' 2 ) are separated by a superselection rule (are 
disjoint), iff their component-state-functions A — >■ |'i/'i(A)) and A — >■ |'i/' 2 (A)) 
have disjoint support on R (up to measure-zero sets). Note that by spec- 
tral decomposition the superselection observables can be decomposed into 
the projectors in which for (16) are all given by multiplications with 
characteristic functions x(A) for cr-measurable sets in R. 



Non-abelian gauge groups We have seen that the fulfillment of Dirac’s 
requirement allows to give a full structural characterisation for the spaces of 
(pure) states and observables. How general is this result? Does it exclude cases 
of physical interest? At first glance this seems indeed to be the case: just con- 
sider a situations with non-abelian gauge groups; for example, the quantum 
mechanical system of n > 2 identical spinless particles with n-particle Hilbert 
space "H = T^(R^") on which the permutation group G = S„ of n objects acts 
in the obvious way by unitary operators U(g). That these particles are iden- 
tical means that observables must commute with each U{g). Without further 
restrictions on observables one would thus define O := {U{g), g G G}'. Hence 
O' is the W*-algebra generated by all U{g), which is clearly non-abelian, thus 
violating (14). But does this generally imply that general particle statistics 
cannot be described in a quantum-mechanical setting which fulfills Dirac’s 
requirement? The answer to this question is ‘no’. Let us explain why. 

If we decompose % according to the unitary, irreducible representations 
of G we obtain ([10,14]) 

p(n) 

H = ^H,, (17) 

i=l 

where i labels the p{n) inequivalent, unitary, irreducible representations Di 
of G of dimension di. Each 'Hi has the structure Hi = 'f^'Hi, where G acts 

irreducibly via Di on and trivially on Hi whereas O acts irreducibly via 
some *-representation on Hi and trivially on . Hi and Hj are inequivalent 
if t yf j. Hence we see that Hi furnishes an irreducible representation for O, iff 
di = i.e., for the Bose and Fermi sectors only. Pure states from these sectors 

are just the rays in the corresponding Hi- In contrast, for di > 1, given a 
non-zero vector \<j)) G Hi, all non-zero vectors in the di-dimensional subspace 
(g) 10) c Hi define the same pure state, i.e., the same expectation-value- 
functional on O. Furthermore, a vector in Hi = ‘Si'Hi which is not a pure 

® It is this irreducibility statement which depends crucially on the fulfillment of 
Dirac’s requirement. In general, the 0(A)’s will act irreducibly on H{X) for each 
A, iff 0“ is maximal abelian in O' , i.e., iff O" = (0“)' n O' . But we already saw 
that (14) also implies 0‘^ = O' so that this is fulfilled. 
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tensor product defines a non-pure state, since the restriction of O G O to 
T~Li is of the form 1 0 O, which means that a vector in TLi defines a state 
given by the reduced density matrix obtained by tracing over the left (i.e. 
C^*) state space. From elementary quantum mechanics we know that the 
resulting state is pure, iff the vector in Hi was a pure tensor product (i.e. 
of rank one). Hence in those Hi where di > I not all vectors correspond to 
pure states, and those which do represent pure states in a redundant fashion 
by higher dimensional subspaces, sometimes called ‘generalized rays’ in the 
older literature on parastatistics [28]. 

However, the factors C^* are completely redundant as far as physical 
information is concerned, which is already fully encoded in the irreducible 
representations of O on Hi; no further physical information is contained 
in di-fold repetitions of tt^. Hence we can define a new, truncated Hilbert 
space 

p(n) 

H:=0H*. (18) 

i=l 

This procedure has also been called ‘elimination of the generalized ray’ in 
the older literature on parastatistics [18] - see also [14] for a more recent 
discussion of this point. Since every pure state in H is also contained in H, just 
without repetition, these two sets are called ‘phenomenological equivalent’ in 
the literature on QFT (e.g. in chapter 6.1.C of [4]). The point is that pure 
states are now faithfully labelled by rays in the Hi and that O' - where the 
commutant is now taken in B(H) rather than B(H) — is generated by 1 and 
the p{n) (commuting!) projectors into the Hi's. Hence Dirac’s requirement is 
satisfied. But clearly the original gauge group has no action on H anymore, 
but there is also no physical reason why one should keep it.^ It served to 
define O, but then only its irreducible representations tt^ are of interest. 
Only a residual action of the center of G still exists, but the gauge group 
generated by the projectors into the Hi consists in fact of the continuous 
group of p{n) copies of C/(l), one global phase change for each sector. Its 
meaning is simply to induce the separation into the different sectors {Hi, tt^), 
and that in accordance with Dirac’s requirement. 

To sum up, we have seen that even if a theory is initially formulated via 
non-abelian gauge groups, we can give it a physically equivalent formula- 
tion that has at most a residual abelian gauge group left and hence obeys 
Dirac’s requirement. Hence the ‘obvious’ counterexamples to Dirac’s require- 
ment turn out to be harmless. This is generally true in quantum mechanics. 



^ In [10] Dirac’s requirement together with the requirement that the physical 
Hilbert space must carry an action of the gauge group has been used to “prove” 
the absence of parastatistics. In our opinion there seems to be no physical reason 
to accept the second requirement and hence the “proof” ; compare [18] and [14] . 
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but in quantum field theory there are genuine possibilities to violate Dirac’s 
condition which we will ignore here.® 

3 Superselection Rules via Symmetry Requirements 

The requirement that a certain group must act on the set of all physical 
states is often the (kinematical) source of superselection rules. Here I wish to 
explain the structure of this argument. 

Note first that in quantum mechanics we identify the states of a closed 
system with rays and not with vectors which represent them (in a redundant 
fashion). It is therefore not necessary to require that a symmetry group G 
acts on the Hilbert space H, but rather it is sufficient that it acts on VH, 
the space of rays, via so-called ray-representations. Mathematically this is a 
non-trivial relaxation since not every ray-representation of a symmetry group 
G (i.e. preserving the ray products) lifts to a unitary action of G on H. What 
may go wrong is not that for a given g G G we cannot find a unitary (or 
anti-unitary) operator Ug on H; that is assured by Wigner’s theorem (see 
[3] for a proof). Rather, what may fail to be possible is that we can choose 
the Ug’s in such a way that we have an action, i.e., that Ug^Llgj = Ug-^g^- 
As is well known, this is precisely what happens for the implementation of 
the Galilei group in ordinary quantum mechanics. Without the admission 
of ray representations we would not be able to say that ordinary quantum 
mechanics is Galilei invariant. 

To be more precise, to have a ray-representation means that for each 
g G G there is a unitary® transformation Ug which, instead of the usual 
representation property, are only required to satisfy the weaker condition 

Ug,Ug, =eXp(i^(6fl,6f2))Ugjg, (l9) 

for some function ^ : G x G — >■ K, called multiplier exponent, satisfying^® 

e(l,ff)=e( 5 ,l) = 0 , ( 20 ) 

^(51,52) - ^{91,9293) + ^{9192,93) - ^{92,93) = 0 . ( 21 ) 

The second of these conditions is a direct consequence of associativity: 

U.,(Ug,Ug3) = (Ug,UgJUg3. 

® An abelian O' implies that G is a von Neumann algebra of type I (see [7] , chap- 
ter 8) whereas truly infinite systems in QFT are often described by type III 
algebras. 

® For simplicity we ignore anti-unitary transformations. They cannot arise if, for 
example, G is connected. 

The following conditions might seem a little too strong, since it would be sufficient 
to require the equalities in (20) and (21) only mod 27t; this also applies to (22). 
But for our application in section 4 it is more convenient to work with strict 
equalities, which in fact implies no loss of generality; compare [32]. 
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Obviously these maps project to an action of G on VH. Any other lift of this 
action on VH onto H is given by a redefinition — >■ Ug := exp(z7(5))Ug, 

for some function 7 : G — >■ M with 7(1) = 0, resulting in new multiplier 
exponents 

^'(51.52) =C(5i, 52) +7(51) -7(3192) +7(52), (22) 

which again satisfy (20) and (21). The ray representations U and are then 
said to be equivalent, since the projected actions on VH are the same. We 
shall also say that two multiplier exponents are equivalent if they satisfy 
(22) for some 7. 

We shall now see how the existence of inequivalent multiplier exponents, 
together with the requirement that the group should act on the space of 
physical states, may clash with the superposition principle and thus give rise 
to superselection rules. For this we start from two Hilbert spaces H' and T~L" 
and actions of a symmetry group G on VH' and VH", i.e., ray representations 
U' and U" on H' and H" up to equivalences (22). We consider % — %' ® H" 
and ask under what conditions does there exist an action of G on VH which 
restricts to the given actions on the subsets VH' and VH" . Equivalently: 
when is U = © U" a ray representation of G on "H for some choice of 

ray-representations and U" within their equivalence class? To answer this 
question, we consider 

U.,Ug, = (U'g^©U")(U'g,©U") 

= exp(i$'(gi,52))Ug^g,^ ®exp(C"(9i,92))U" g^ (23) 

and note that this can be written in the form (19), for some choice of 
within their equivalence class, iff the phase factors can be made to coin- 
cide, that is, iff and are equivalent. This shows that there exists a 
ray-representation on H which restricts to the given equivalence classes of 
given ray representations on H' and H' , iff the multiplier exponents of the 
latter are equivalent. Hence, if the multiplier exponents and are not 
equivalent, the action of G cannot be extended beyond the disjoint union 
VH' UVH". Conversely, if we require that the space of physical states must 
support an action of G, then non-trivial superpositions of states in H' and 
H" must be excluded from the space of (pure) physical states. 

This argument shows that if we insist of implementing G as symmetry 
group, superselection rules are sometimes unavoidable. A formal trick to avoid 
them would be not to require G, but a slightly larger group, G, to act on the 
space of physical states. G is chosen to be the group whose elements we label 
by (0,g), where 6* G R, and the multiplication law is 

9192 = {0i,9i){S2,92) = {01 +02 + f {91, 92), 9192)- (24) 

It is easy to check that the elements of the form {9, 1) lie in the center of G 
and form a normal subgroup = M which we call Z. Hence GjZ = G but G 
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need not be a subgroup of G. G is a central R extension^^ of G (see e.g. 
[32]). Now a ray-representation U of G on "H defines a proper representation 
[/ of G on "H by setting 

:= exp(i6»)Ug. (25) 

Then G is properly represented on %' and %" and hence also on 'H = ■ 

The above phenomenon is mirrored here by the fact that Z acts trivially on 
VT-L' and VTi!' but non-trivially on VT-L, and the superselection structure 
comes about by requiring physical states to be fixed points of Z's action. 

4 Bargmann’s Superselection Rule 

An often mentioned textbook example where a particular implementation 
of a symmetry group allegedly clashes with the superposition principle, such 
that a superselection rule results, is Galilei invariant quantum mechanics (e.g. 
[9]; see also Wightman’s review [39]). We will discuss this example in detail 
for the general multi-particle case. (Textbook discussions usually restrict to 
one particle, which, due to Galilei invariance, must necessarily be free.) It will 
serve as a test case to illustrate the argument of the previous chapter and also 
to formulate our critique. Its physical significance is limited by the fact that 
the particular feature of the Galilei group that is responsible for the existence 
of the mass superselection rule ceases to exist if we replace the Galilei group 
by the Poincare group (i.e. it is unstable under ‘deformations’). But this is not 
important for our argument. Let now G be the Galilei group, an element 
of which is parameterized by (i?, v, a, b), with R a rotation matrix in 5*0(3), 
V the boost velocity, a the spatial translation, and b the time translation. Its 
laws of multiplication and inversion are respectively given by 

9i92 = {Ri,Vi,ai,bi){R2,V2, a2,b2) 

= (i?ii ?2 , Vi + Ri ■ V 2 , ai + Ri ■ tt 2 + Vib 2 , &i -k & 2 ), (26) 

g~^ = {R,v,a,b)~^ = (R~^, —R~^ ■ V , —R~^ ■ {a — vb) , —b). (27) 

We consider the Schrodinger equation for a system of n particles of positions 
Xi, masses rrii, mutual distances := |]a:j — ®j|| which interact via a Galilei- 
invariant potential V{{rij}), so that the Hamilton operator becomes H = 
—h^ + V. The Hilbert space is H = L^(M^",d^Xi • • -d^Xn). 

G acts on the space {configurations} x {times} = as follows: Let g = 

(i?, V, a, b), then g{{xi}, t) := {{R-Xi+vt+a\ , t+b). Hence G has the obvious 

Had we defined the multiplier exponents mod 27 t (compare footnote 10) then we 
would have obtained a 17(1) extension, which would suffice so far. But in the 
next section we will definitively need the R extension as symmetry group of the 
extended classical model discussed there. 

In General Relativity, where the total mass can be expressed as a surface integral 
at ‘infinity’, the issue of mass superselection comes up again; see e.g. [15] and [8]. 
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left action on complex- valued functions on However, 

these transformations do not map solutions of the Schrodinger equations into 
solutions. But, as is well known, this can be achieved by introducing an 
R^"^i-dependent phase factor (see e.g. [13] for a general derivation). We set 
M = the total mass and J'c = tor the center-of-mass. 

Then the modified transformation, Tg, which maps solutions (i.e. curves in 
7i) to solutions, is given by 

:=exp(iM[v- (r^-a) - ^v‘^{t-b)]) (28) 

However, due to the modification, these transformations have lost the prop- 
erty to define an action of G, that is, we do not have Tg^ o Tg^ = Tg^g^. 
Rather, a straightforward calculation using (26) and (27) leads to 

Tgi oTg2 = exp(i^(5i,52))Tg,g3, (29) 

with non-trivial multiplier exponent 

^{ 91 , 92 ) = ^{vi- Ri- tt 2 + ^vlb 2 ). (30) 

Although each Tg is a mapping of curves in "H, it also defines a unitary 
transformation on T~L itself. This is so because the equations of motion define 
a bijection between solution curves and initial conditions at, say, t = 0, which 
allows to translate the map Tg into a unitary map on R, which we call Ug. 
It is given by 

= exp {^M[v ■ {tc - a) -I- \v‘^b]) exp(iH5)V’({R“^(®i-a+i’6)}), 

( 31 ) 

and furnishes a ray-representation whose multiplier exponents are given by 
(30). It is easy to see that the multiplier exponents are non-trivial, i.e., not 
removable by a redefinition (22). The quickest way to see this is as follows: 
suppose to the contrary that they were trivial and that hence (22) holds 
with ^' = 0. Trivially, this equation will continue to hold after restriction to 
any subgroup Gq C G. We choose for Go the abelian subgroup generated by 
boosts and space translations, so that the combination 7(gi) — 7 ( 5 i(? 2 ) + 7 (ff 2 ) 
becomes symmetric in 51,(72 G Gq. But the exponent (30) stays obviously 
asymmetric after restriction to Gq. Hence no cancellation can take place, 
which contradicts our initial assumption. 

The same trick immediately shows that the multiplier exponents are in- 
equivalent for different total masses M . Hence, by the general argument given 
in the previous chapter, if %' and %" correspond to Hilbert spaces of states 
with different overall masses M' and M" , then the requirement that the 
Galilei group should act on the set of physical states excludes superpositions 
of states of different overall mass. This is Bargmann’s superselection rule. 

I criticize these arguments for the following reason: The dynamical frame- 
work that we consider here treats ‘mass’ as parameter(s) which serves to 
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specify the system. States for different overall masses are states of differ- 
ent dynamical systems, to which the superposition principle does not even 
potentially apply. In order to investigate a possible violation of the super- 
position principle, we must find a dynamical framework in which states of 
different overall mass are states of the same system; in other words, where 
mass is a dynamical variable. But if we enlarge our system to one where mass 
is dynamical, it is not at all obvious that the Galilei group will survive as 
symmetry group. We will now see that in fact it does not, at least for the 
simple dynamical extension which we now discuss. 

The most simple extension of the classical model is to maintain the Hamil- 
tonian, but now regarded as function on an extended, 6n -I- 2n - dimensional 
phase space with extra ‘momenta’ mi and conjugate generalized ‘positions’ 
\i- Since the A^’s do not appear in the Hamiltonian, the m^’s are constants 
of motion. Hence the equations of motion for the XiS and their conjugate 
momenta Pi are unchanged (upon inserting the integration constants mi) and 
those of the new positions are 






dmi 2m‘f ’ 



(32) 



which, upon inserting the solutions {xi{t),Pi(t)}, are solved by quadrature. 

Now, the point is that the new Hamiltonian equations of motion do not 
allow the Galilei group as symmetries anymore. But they do allow the R- 
extension G as symmetries [13]. Its multiplication law is given by (24), with ^ 
as in (30). The action of G on the extended space of {configurations} x {times} 
is now given by 



g{{x,},{X,},t) = {9,R,v,a,b){{xff,{Xff,t) 

= {{Rxi -\-vt-\- a} , {Xi - {-^9 V ■ R - Xi , t -I- 6). (33) 

With (24) and (30) it is easy to verify that this defines indeed an action. 
Hence it also defines an action on curves in the new Hilbert space % := 
given by 

j gip := Ip o , (34) 

which already maps solutions of the new Schrodinger equation to solutions, 
without invoking non-trivial phase factors. This is seen as follows: Let 

<l'{{xi},{Xi},t) G U 

and (!>{{xi},{mi},t) its Fourier transform in the (Xi,mi) arguments: 



-n/2 



d”m exp 



. n 

2=1 



'P{{x,}, {A*},t) = {2nh) 



<P{{xi},{mi},t). 

( 35 ) 
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For each set of masses {nn} the function := 

satisfies the original Schrodinger equation. Since (34) does not mix different 

sets of {nii} it induces a map for each such set: 

: = exp [i6» + j^M {v ■ (r^ - a) - ^v'^{t-b))] 

X (36) 

Via the Fourier transform (35) we represent H as direct integral of "H{mj’s, 
each of which isomorphic to our old % = L^(M^", d^a^i • • • d^a;„), and on each 
of which (36) defines a unitary representation U oi G the form (25) with Ug 
the ray-representation (31). This shows how the much simpler transformation 
law (34) contains the more complicated one (28) upon writing "H as a direct 
integral of vector spaces 

In the new framework the overall mass, M, is a dynamical variable, and 
it would make sense to state a superselection rule with respect to it. But now 
G rather than G is the dynamical symmetry group, which acts by a proper 
unitary representation on %, so that the requirement that the dynamical 
symmetry group should act on the space of physical states will now not lead 
to any superselection rule. Rather, the new and more physical interpretation 
of a possible superselection rule for M would be that we cannot localize 
the system in the coordinate conjugate to overall mass, which we call A, 
i.e., that only the relative new positions \ — Xj are observable. (This is 
so because M generates translations of equal amount in all Ai.) But this 
would now be a contingent physical property rather than a mathematical 
necessity. Note also that in our dynamical setup it is inconsistent to just state 
that M generates gauge symmetries, i.e. that A corresponds to a physically 
non existent degree of freedom. For example, a motion in real time along 
A requires a non- vanishing action (for non- vanishing M), due to the term 
f dtMAva the expression for the action. 

If decoherence were to explain the (ficticious) mass superselection rule, it 
would be due to a dynamical instability (as explained in [24]) of those states 
which are more or less localized in A. Mathematically this effect would be 
modelled by effectively removing the projectors onto yl-subintervalls from the 
algebra of observables, thereby putting M (i.e. its projectors) into the center 
of O. Such a non-trivial center should therefore be thought of as resulting 
from an approximation-dependent idealisation. 



A system {(Ai, fhi}) of canonical coordinates including M = 5()i rrii is e.g. Ai : = 
Ai, rhi = M and Xi = Xi — Xi, rhi = rm for i = 2...n. Then T = Ai. 
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5 Charge Superselection Rule 



In the previous case I said that superselection rules should be stated within a 
dynamical framework including as dynamical degree of freedom the direction 
generated by the superselected quantity. What is this degree of freedom in 
the case of a superselected electric charge and how does it naturally appear 
within the dynamical setup? What is its relation to the Coulomb field whose 
role in charge-decoherence has been suggested in [15]? In the following discus- 
sion I wish to investigate into these questions by looking at the Hamiltonian 
formulation of Maxwell’s equation and the associated canonical quantization. 

In Minkowski space, with preferred coordinates {x^ = z)} (labo- 

ratory rest frame), we consider the spatially finite region Z = {{t,x,y,z) : 
x"^ + y'^ + z"^ < R^}. E denotes the intersection of Z with a slice t = const, 
and dS =: Sr its boundary (the laboratory walls). Suppose we wish to solve 
Maxwell’s equations within Z, allowing for charged solutions. It is well known 
that in order for charged configurations to be stationary points of the action, 
the standard action functional has to be supplemented by certain surface 
terms (see e.g. [11]) which involve new fields on the boundary, which we call 
A and /, and which represent a pair of canonically conjugate variables in 
the Hamiltonian sense. On the laboratory walls, dS, we put the boundary 
conditions that the normal component of the current and the tangential com- 
ponents of the magnetic field vanish. Then the appropriate boundary term 
for the action reads 



dtduj{X + <p)f, 



(37) 



where (j) is the scalar potential and ckv the measure on the spatial boundary 2- 
sphere rescaled to unit radius. Adding this to the standard action functional 
and expressing all fields on the spatial boundary by their multipole moments 
(so that integrals /g^dwi?^, dw = measure on unit sphere, become 
one arrives at a Hamiltonian function 



H= / [1(£;2 + (Vx A)2) + <^(p- V-£;)-A-j] + V</>z^(Az^-/i™). 

Im 

(38) 

Here the pairs of canonically conjugate variables are {A{x),—E{x)) and 
(Aim, fim), and Eim are the multipole components of n ■ E, 



Elm ■= [ dujR^ Yim.n ■ E, (39) 

JdS 

where n is the normal to dE. The scalar potential 4> has to be considered as 
Lagrange multiplier. With the given boundary conditions the Hamiltonian is 
differentiable with respect to all the canonical variables^"* and leads to the 



14 



This would not be true without the additional surface term (37). Without it 
one does not simply obtain the wrong Hamiltonian equations of motions, but 
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following equations of motion 



A 

-E 



^Im 






m 



r Tj 

=j-Vx(VxA) 
dH 



dflr, 



— 4^lm j 



dH 

dXlm 



= 0 . 



(40) 

(41) 

(42) 

(43) 



These are supplemented by the equations which one obtains by varying with 
respect to the scalar potential (f), which, as already said, is considered as 
Lagrange multiplier. Varying first with respect to 4>{x) (i.e. within E) and 
then with respect to (pim (i-e. on the boundary dE), one obtains 



G(x) : = V • E{x) — p{x) = 0, (44) 

Gi^-. =Eim-flm = 0. (45) 

These equations are constraints (containing no time derivatives) which, once 
imposed on initial conditions, continue to hold due to the equations of motion. 

15 

This ends our discussion of the classical theory. The point was to show that 
it leaves no ambiguity as to what its dynamical degrees of freedom are, and 
that we had to include the variables Xim along with their conjugate momenta 
fim in order to gain consistency with the existence of charged configurations. 
The physical interpretation of the Xim’s is not obvious. Equation (42) merely 
relates their time derivative to the scalar potential’s multipole moments on 
the boundary, which are clearly highly non-local quantities. The interpreta- 
tion of the fim’s follow from (45) and the definition of Eim, i-e. they are the 
multipole moments of the electric flux distribution (p{n) := R^n ■ E{R^n). 
In particular, for / = 0 = m we have 

/oo = (47r)"5 Q, (46) 

none at all! Concerning the Langrangean formalism one should be aware that the 
Euler-Lagrange equations may formally admit solutions (e.g. with long-ranged 
(charged) fields) which are outside the class of functions which one used in the 
variational principle of the action (e.g. rapid fall-off). Such solutions are not 
stationary points of the action and their admittance is in conflict with the vari- 
ational principle unless the expression for the action is modified by appropriate 
boundary terms. 

Equation (41) together with charge conservation, p -|- V • j = 0, shows that (44) 
is preserved in time, and (43) together with the boundary condition that n • j 
and n x (V x A) vanish on dE show that (45) is preserved in time. 
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where Q is the total charge of the system. Hence we see that the total charge 
generates motions in Aqo • But this means that the degree of freedom labelled 
by Aqo truly exists (in the sense of the theory) . For example, a motion along 
Aqo will cost a non- vanishing amount of action oc Q(Ao“' — A dec- 

laration that Ago really labels only a gauge degree of freedom is incompatible 
with the inclusion of charged states. Similar considerations apply of course to 
the other values of l,m. But note that this conclusion is independent of the 
radius R of the spatial boundary 2-sphere dS. In particular, it continues to 
hold in the limit i? — >■ oo. We will not consistently get rid of physical degrees 
of freedom that way, even if we agree that realistic physical measurements 
will only detect field values in bounded regions of space-time. See [12] for 
more discussion on this point and the distinction between proper symmetries 
and gauge symmetries. 

It should be obvious how these last remarks apply to the statement of 
a charge superselection rule. Without entering the technical issues (see e.g. 
[33]), its basic ingredient is Gauss’ law (for operator- valued quantities), lo- 
cality of the electric field and causality. That Q commutes with all (quasi-) 
local observables then follows simply from writing Q as surface integral of 
the local flux operator R^n ■ E, and the observation that the surface may 
be taken to lie in the causal complement of any bounded space-time region. 
Causality then implies commutativity with any local observable. 

In a heuristic Schrodinger picture formulation of QED one represents 
states W by functions of the configuration variables A{x) and Xim- The mo- 
mentum operators are obtained as usual: 



-E{x) - 


. 6 

5A{x) 


(47) 


flm 


■ ^ 
^dXim ' 


(48) 



In particular, the constraint (45) implies the statement that on physical states 
'R we have^® 

QtF = -iVdTT-^tf' . (49) 

oXoo 

This shows that a charge superselection rule is equivalent to the statement 
that we cannot localize the system in its Agg degree of freedom. Removing by 
hand the multiplication operator Agg (i.e. the projectors onto Agg-intervals) 
from our observables clearly makes Q a central element in the remaining 
algebra of observables. But what is the physical justification for this removal? 
Certainly, it is valid FAPP if one restricts to local observations in space-time. 
To state that this is a fundamental restriction, and not only an approximate 

Clearly all sorts of points are simply sketched over here. For example, charge 
quantization presumably means that Ago should be taken with a compact range, 
which in turn will modify (48) and (49) . But this is irrelevant to the point stressed 
here. 
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one, is equivalent to saying that for some fundamental reason we cannot 
have access to some of the existing degrees of freedom, which seems at odds 
with the dynamical setup. Rather, there should be a dynamical reason for 
why localizations in Aqo seem FAPP out of reach. The idea of decoherence 
would be that localizations in Aqo are highly unstable against dynamical 
decoherence. 

We have mainly focussed on the charge superselection operator /oo, al- 
though the foregoing considerations make it clear that by the same argument 
any two different asymptotic flux distributions also define different supers- 
election sectors of the theory. Do we expect these additional superselection 
rules to be physically real? First note that for I > 0 the fim are not directly 
related to the multipole moments of the charge distributions, as the latter 
fall-off faster than ^ and are hence not detectable on the sphere at infin- 
ity. Conversely, the higher multipole moments fim are not measurable (in 
terms of electromagnetic fields) within any finite region of space-time, unlike 
the charge, which is tight to massive particles; any finite sphere enclosing all 
sources has the same total flux. But the fim can be related to the kinematical 
state of a particle through the retarded Coulomb field. In fact, given a par- 
ticle with constant momentum p, charge e and mass m, one obtains for the 
electric flux distribution at time t on a sphere centered at the instantaneous 
(i.e. at time t) particle positiond^ 



V3p(n) 



em^ [p^ + m^] 2 
47t [(p . n)^ -I- m?] i 



(51) 



Hence different incoming momenta would induce different flux distributions 
and therefore lie in different sectors. Given that these sectors exist this means 

Formula (51) requires a little more explanation: for a particle with general tra- 
jectory z{t) let t' be the retarded time for the space-time point (x,t), i.e., 
F = t — 11® — 2(F) II (c = 1 in our units). Now we can use the well known 
formula for the retarded electric field (e.g. (14.14) in [21]) and compute the flux 
distribution on a sphere which lies in the space of constant time t, where it is 
centered at the retarded position z{t') of the particle. This flux distribution can 
be expressed as function of the retarded momentum p' := p{t') and the retarded 
direction n' := [x — 2 :(t')]/||® — 2 (F)|| as follows {E' \/p'^ + m^): 



V'p'in') 



em^ 1 

47t [E' — p' ■ n']2 



(50) 



If the particle moves with constant velocity v := z, the expression for the retarded 
Coulomb held can be rewritten in terms of the instantaneous position z{t) by 
using z{t) — z{t') -b u||® — z{t')\\. With respect to this center it is purely radial. 
Then one calculates the hux distribution on a sphere which again lies in the space 
of constant time t, but now centered at z{t) rather than z{t'). This function can 
be expressed in terms of the instantaneous direction n := [® — z(t)]/||® — z{t)\\ 
and the instantaneous momentum p pit). One obtains (51). 
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that different incoming momenta cannot be coherently superposed and no 
incoming localized states be formed, unless one also adds the appropriate 
incoming infrared photons to just cancel the difference of the asymptotic flux 
distributions. This is achieved by imposing the ‘infrared coherence condition 
of Zwanziger [41]^® the effect of which is to ‘dress’ the charged particles with 
infrared photons which just subtract their retarded Coulomb fields at large 
spatial distances. Hence coherent superpositions of particles with different 
momenta can only be formed if they are dressed by the right amount of 
incoming infrared photons. 

As a technical aside we remark that this can be done without violating the 
Gupta-Bleuler transversality condition k'^a^{k)\m) = 0 in the zero-frequency 
limit, precisely because of the surface term (37) [11]. This resolved an old is- 
sue concerning the compatibility of the infrared coherence condition on one 
hand, and the Gupta-Bleuler transversality condition on the other [17,42]. 
From what we said earlier concerning the consistency of the variational prin- 
ciple in the presence of charged states, such an apparent clash of these two 
conditions had to be expected: without the surface variables one cannot main- 
tain gauge invariance at spatial infinity (i.e. in the infrared limit) and at the 
same time include charged states. In the charged sectors the longitudinal 
infrared photons correspond to real physical degrees of freedom and it will 
naturally lead to inconsistencies if one tries to eliminate them by imposing 
the Gupta-Bleuler transversality condition also in the infrared limit. How- 
ever, a gauge symmetry in the infrared limit can be maintained if one adds 
the asymptotic degrees of freedom in the form of surface terms. 

These remarks illustrate how the rich superselection structure associ- 
ated with different asymptotic flux distributions fim renders the problem 
of characterizing state spaces in QED for charged sectors fairly complicated. 
This problem has been studied within various formalisms including algebraic 
QFT [5] and lattice approximations, where the algebra of observables can be 
explicitly presented [25]. However, all this takes for granted the existence of 
the superselection rules, whereas we would like to see whether they really 
arise from some physical impossibility to localize the system in the degrees 
of freedom labelled by Xim- What physics should prevent us from forming 
incoming localized wave packets of charged undressed (in the sense above) 
particles, which would produce coherent superpositions of asymptotic flux 
distributions from the sectors with I > 1? This cries out for a decoherence 
mechanism to provide a satisfying physical explanation. The case of charge 
superselection is, however, more elusive, since localizations in Aqo do not have 
an obvious physical interpretation. Gompare the controversy between [1,29] 
on one side and [36] on the other. 



18 



Basically it says that the incoming scattering states should be eigenstates to the 
photon annihilation operators a)j’(fc) in the zero-frequency limit. 
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Abstract. Classical mechanics and standard Copenhagen quantum mechanics re- 
spect subspace implicatious. For example, if a particle is coufined in a particular 
region R of space, then in these theories we can deduce that it is confined in regions 
containing R. However, subspace implications are generally violated by versions of 
quantum theory that assign probabilities to histories, such as the consistent his- 
tories approach. I define here a new criterion, ordered consistency, which refines 
the criterion of consistency and has the property that inferences made by ordered 
consistent sets do not violate subspace relations. This raises the question: do the op- 
erators defining our observations form an ordered consistent history? If so, ordered 
consistency defines a version of quantum theory with greater predictive power than 
the consistent histories formalism. If not, and our observations are defined by a non- 
ordered consistent quantum history, then subspace implications are not generally 
valid. 



1 Introduction 

We take it for granted that we can infer quantitatively less precise statements 
from our observations. For example, if we know that an atom is confined in 
some region R of space, we believe we are free to assume for calculational 
purposes only that it lies in some larger region containing R. Our under- 
standing of the world and our interpretation of everyday experience tacitly 
rely on subspace implications of this general type: if a physical quantity is 
known to lie within a range R\, then it lies in all ranges R 2 Z) R\. 

In classical mechanics, subspace inferences follow from the correspondence 
between physical states and points in phase space: if the state of a system 
lies in a subset 5*1 of phase space, it lies in all subsets S 2 D -S'!. Similarly, 
in Copenhagen quantum theory, they hold since if the state of a quantum 
system lies in a subspace Hi of Hilbert space, it lies in all subspaces H 2 Z) Hi. 

However, neither classical mechanics nor (presumably) Copenhagen quan- 
tum theory is fundamentally correct. If the basic principles of quantum theory 
apply to the universe as a whole, then a post-Copenhagen interpretation of 
quantum theory seems to be needed, and any justification of the subspace im- 
plications must ultimately be given in terms of that interpretation. Though it 
may seem hard to imagine how to make sense of nature without allowing sub- 
space inferences, there are versions of quantum theory in which they do not 
hold. In particular, this is true of recent attempts to develop history-based 
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formulations of quantum theory [1-5], which rely on the notion of consistent 
or decoherent sets of histories. 

This paper suggests a way in which quantum theory can plausibly be 
interpreted via statements about histories, without violating subspace impli- 
cations — the motivation being that both history-based interpretations and 
subspace implications seem natural. The interpretation which results, a re- 
finement of the consistent histories approach based on ordered consistent sets 
of histories, is certainly not the ultimate answer to the problems of quantum 
theory. It does not solve the measurement problem, for example. But perhaps 
it is a step in the right direction, or at least in an interesting direction. It 
also helps to make precise the question as to what it would mean if subspace 
implications actually were violated in nature, which the last part of the paper 
examines. 

The language of quantum histories may not necessarily be the right way to 
interpret the quantum theory of closed systems. Bohm theory [6] and dynam- 
ical collapse theories [7] are particularly interesting alternatives, for example. 
Quantum histories still play a role in these theories as usually interpreted, 
but it is certainly possible to interpret them so as to ascribe no ontologi- 
cal status to the evolving quantum state, in which case quantum histories 
have no ontological status either, and the orderings amongst them become 
physically irrelevant. 

In this paper, which focusses on the relation between quantum histories 
and orderings, I take it for granted — as an interesting hypothesis, not a 
dogma — that quantum histories have some ontological status. I will focus 
on the consistent histories approach. This is not to suggest that there are no 
other interesting quantum histories approaches. Other ideas are outlined in 
Ref. [8], Section VII of Ref. [9], and Ref. [10], for example. But the discussion 
of this paper is best carried out within some specific formulation of quantum 
theory, and consistent histories is the most widely studied example. As an 
approach to closed system quantum histories, it has many interesting fea- 
tures. Its main drawback is that it does not solve the measurement problem, 
which in the language of consistent histories takes the form of a set selec- 
tion problem. The search for ways to refine the definition of consistency in 
order to solve, or at least reduce, the set selection problem is an independent 
motivation for the definition of ordered consistency proposed below. 



2 Partial Ordering of Quantum Histories 

We consider versions of quantum theory that assign probabilities to histories 
— i.e., to collections of events. There are several reasonably standard ways 
of representing an event in quantum theory. The simplest is as a projection 
operator, labelled by a particular time, corresponding to a statement about 
an observable in non-relativistic quantum mechanics. Thus, the statement 
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that a particle was in the interval I at time t is represented by the projection 

/ \x){x\d^x. (1) 

Jx^l 

The representation of events as projections can be generalised to quantum ef- 
fects [10-12], defined by operators A such that A and I— A are both positive.^ 
Events can also be defined, at least formally, in the path integral version of 
quantum theory by dividing the set of paths into exclusive subsets. For exam- 
ple, by considering the appropriate integrals one can attach probabilities to 
the events that a particle did, or did not, enter a particular region of space- 
time [13]. Further generalisations have been discussed by Isham, Linden and 
collaborators [5,14,15]. 

Each of these representations has a natural partial ordering. For projec- 
tions, we take A > B A the range of i? is a subspace of that of A. This 
corresponds to the logical implication of Copenhagen quantum theory al- 
ready mentioned: if the state of a quantum system lies in the range of B, 
then it necessarily lies in the range of A. For quantum effects, we take A> B 
if and only if (A — B) is positive. For events defined by the position space 
path integral, we take the statement that a particle entered R\ to imply the 
statement that it entered i ?2 if Ri C i? 2 . 

We define a quantum history as a collection of quantum events and ex- 
tend the partial ordering to histories in the natural way. For example, if 
we take quantum histories in non-relativistic quantum mechanics to be de- 
fined by sequences of projections at different fixed times, we compare two 
quantum histories H = {Ai,ti, ■ ■ ■ ', An;tn} and H' = . . . ; t„,} 

as follows. First add to each history the identity operator at every time 
at which it contains no proposition and the other does, so obtaining rela- 
belled representations of the histories as H = {Ai, Ti; . . . ; A^; T^} and H' = 
{A'l.Ti', . . . ] A'f^, Tm}- Then define the partial ordering by taking H > H' ii 
and only if Ai > A'^ for all i and H > H' \i H > H' and H ^ H' . This 
ordering was first considered in the consistent histories literature by Isham 
and Linden [14], whose work, and its relation to the ideas of this paper, are 
discussed in the Appendix. 

3 Consistent Histories 

The consistent sets of histories for a closed quantum system are defined in 
terms of the space of states T~L, the initial density matrix p, and the hamil- 
tonian H. In the simplest version of the consistent histories formulation of 
non-relativistic quantum mechanics, sets of histories correspond to sets of pro- 
jective decompositions. In order to be able to give a physical interpretation 
of any of the consistent sets, we need also to assume that standard observ- 
ables, such as position, momentum and spin, are given. Individual quantum 

^ Different definitions can be found: this, the simplest, is adequate for our purposes. 
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events are defined by members of projective decompositions of the identity 
into orthogonal hermitian projections a = {-P*}, with 

^ P* = 1 and P*P^' = %P* . (2) 

i 

Each such decomposition defines a complete and exclusive list of events at 
some fixed time, and a time label is thus generally attached to the decompo- 
sitions and the projections: the labels are omitted here, since the properties 
of interest here depend only on the time ordering of events. 

Suppose now we have a set of decompositions S = {cti, . . . , cr„}. Then the 
histories given by choosing one projection from each of the decompositions 
aj in all possible ways define an exhaustive and exclusive set of alternatives. 
We follow Gell-Mann and Hartle’s definitions, and say that 5 is a consistent 
(or medium decoherent) set of histories if 

Tr(P;- . . . PlVPf^ . . . = 6,,,, . . . ...in). (3) 

When S is consistent, p{ii . . .in) is the probability of the history H = 
{P{^ , . . . We will want later to discuss the properties of individual 

histories without reference to any fixed consistent set, and we define a con- 
sistent history to be a history which belongs to some consistent set S. Finally, 
we say the set 

is a consistent extension of a consistent set of histories S = {ai, . . . , cr„} by 
the set of projections r = {Q* : i = 1, . . . , m} if t is a projective decomposi- 
tion and S' is consistent. 

Suppose now that we have a collection of data defined by the history 

H = {P{\...,Pl-} (5) 

which has non-zero probability and belongs to the consistent set S. This his- 
tory might, for example, describe the results of a series of experiments or 
the observations made by an observer. Given a choice of consistent extension 
S' oi S, we can make probabilistic inferences conditioned on the history H . 
For example, if S' has the above form, the histories extending H in S' are 
H' = . . . , Q*, . . . , and the proposition Q' has condi- 

tional probability 

p{Q'\H)=p{H')/p{H). (6) 

These probabilities and conditional probabilities are the same in every con- 
sistent set which includes H' and (hence) H . However, when we want to em- 
phasize that the calculation can be carried out in some particular set S, we 
will attach a suffix. For example, we might write Ps'{Q'\H) = Ps'{H') /ps{H) 
for the above equation. 

The formalism itself gives no way of choosing any particular consistent 
extension. In the view of the original developers of the consistent histories 
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approach, the different S' are to be thought of as ways of producing possible 
pictures of the past and future physics of the system which, though generally 
incompatible, are all equally valid. More formally, they can be seen as incom- 
patible logical structures which allow different classes of inferences from the 
data [16,17]. 

It is this freedom in the choice of consistent extension which, it has been 
argued elsewhere [9,18-23,10] gives rise to the most serious problems in the 
consistent histories approach. Standard probabilistic predictions and deter- 
ministic retrodictions can be reproduced in the consistent histories formalism 
by making an ad hoc choice of consistent set , but cannot be derived from the 
formalism itself. In fact, it is almost never possible to make any unambigu- 
ous predictions or retrodictions: there are almost always an infinite number 
of incompatible consistent extensions of the set containing a given history 
dataset [9,20]. 

The problem is not simply that the formalism supplies descriptions of 
physics which are complementary in the standard sense, although that in itself 
is sufficient to ensure that the formalism is only very weakly predictive. Even 
on the assumption that we will continue to observe quasiclassical physics, 
no known interpretation of the formalism allows us to derive the predictions 
of classical mechanics and Copenhagen quantum theory [18]. Hence the set 
selection problem: probabilistic predictions can only be made conditional on 
a choice of a consistent set, yet the consistent histories formalism gives no 
way of singling out any particular set or sets as physically interesting. 

One possible solution to the set selection problem would be an axiom 
which identifies a unique physically interesting set, or perhaps a class of 
such sets, from the initial state and the dynamics. Another would be the 
identification of a physically natural measure on the space of consistent sets, 
according to which the physically relevant consistent set is randomly chosen. 
Possible set selection criteria have been investigated [24,25], but no generally 
workable criterion has emerged. ^ 



4 Consistent Sets and Contrary Inferences: A Brief 
Review 

A further reason for believing that the consistent histories formalism is at 
best incomplete comes from considering the logical relations among events in 
different consistent sets. We say that two projection operators P and Q are 
complementary if they do not commute: PQ yf QP. We say that they are 
contradictory if they sum to the identity, so that P = 1 — Q and PQ = QP = 

^ With different motivations, Gell-Mann and Hartle have explored a “strong deco- 
herence” criterion which is intended to rednce the nnmber of consistent sets [26] . 
However, under their present definition, every consistent set is strongly decoher- 
ent [10]. 
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0, and that they are contrary if they are orthogonal and not contradictory, 
so that PQ = QP = 0 and P <1 — Q. 

Contradictory inferences are never possible in the consistent histories for- 
malism, but it is easy to find examples of contrary inferences from the same 
data [19]. For instance, consider a quantum system whose hamiltonian is zero 
and whose Hilbert space % has dimension greater than two, prepared in the 
state I a). Suppose that the system is left undisturbed from time 0 until time 
t, when it is observed in the state |c), where 0 < |(a | c)| < 1/3, i.e. a single 
quantum measurement is made, and the outcome probability is less than 1/9. 
In consistent histories language, we have the initial density matrix p = |a)(a| 
and the single datum corresponding to the history P[ = {Pc} from the con- 
sistent set S = {{Pc, 1 — Pc}}, where the projection Pc = |c)(c| is taken at 
time t. 

Now consider consistent extensions of S of the form Sb = {{Pb, 1 — 
Pft}, {Pc, 1 — Pc}}, where Pf, = |6)(6| for some normalised vector \b) with 
the property that 

{c\b){b\a) = (cja). (7) 

It is easy to verify that Sb is consistent and that the conditional probability 
of Pb given iJ is 1. It is also easy to see that there are at least two mu- 
tually orthogonal vectors \b) satisfying (7) . For example, let jui), [^ 2 ), jus) 
be orthonormal vectors and take |a) = jui) and |c) = Ajui) -I- /i|u 2 ), where 

+ ImP = 1- Then the vectors 

It,) = (lAP + M!)-./.(A|„) + ± I.,)) (8) 

X XX 

both satisfy (7) and are orthogonal if x is real and = (x — 2)(1 — |Ap), 

which has solutions for |A| < 1/3. Thus the consistent sets Sb^ give contrary 
probability on retrodictions. 

Some brief historical remarks are in order. The existence of contrary infer- 
ences in the consistent histories formalism, though easy to show, was noticed 
only quite recently. In particular, it was not known to the formalism’s orig- 
inal developers.^ It was first explicitly pointed out, and its implications for 
the consistent histories formalism were first examined, in Ref. [19]. Further 
discussion can be found in Refs. [27,28]. However, a noteworthy earlier consis- 
tent histories analysis of an example in which contrary inferences arise can be 
found in a critique by Cohen [29] of Aharonov and Vaidman’s interpretation 
[30] of one of their intriguing examples of pre- and post-selection."* As noted 
in Ref. [19], Cohen’s analysis miscontrues the consistency criterion: however, 
this error does not affect its derivation of contrary inferences. 

® I am grateful to Bob Griffiths, Jim Hartle, and Roland Omnes for helpful corre- 
spondence on this point. 

* I am grateful to Oliver Cohen and Lucien Hardy for drawing this reference to my 
attention. 
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Now, the existence of contrary inferences in the consistent histories for- 
malism needs to be interpreted with care. It is not true that, in any given 
consistent set, two different contrary propositions can be inferred with prob- 
ability one. The inferences made within any given consistent set lead to 
no contradiction. The picture of physics given by any given consistent set 
may or may not be considered natural or plausible — depending on one’s 
intuition and the criteria one uses for naturality — but it is not logically 
self-contradictory. It is however true, as a mathematical statement about the 
properties of the consistent histories formalism, that the propositions inferred 
in the two different sets correspond to contrary projections. The formalism 
makes no physical distinction among different consistent sets, and so requires 
us to conclude that two equally valid pictures of physics can be given, in 
which contrary events take place. 

To put it more formally, the consistent histories approach can be inter- 
preted as setting out rules of reasoning according to which, although physics 
can be described by any of infinitely many equally valid pictures, only one of 
those pictures may be considered in any given argument. Such an interpre- 
tation ensures — tautologically — that no logical contradiction arises, even 
when the pictures contain contrary inferences. 

The consistent histories formalism, in other words, gives a set of rules 
for producing possible pictures of physics within quantum theory, and these 
rules themselves lead to no logical inconsistency. However, consistent histo- 
rians claim much more, arguing that the formalism defines a natural and 
scientifically unproblematic interpretation of quantum theory. Indeed, the 
consistent histories literature tends to suggest that the descriptions of physics 
given by consistent historians are simply and evidently the correct descrip- 
tions which emerge from quantum theory, so that, in querying them, one 
necessarily queries quantum theory itself.® 

This seems patently false. The most basic premise of consistent historians 
— that quantum theory is correctly interpreted by some sort of many-picture 
scheme — leads to such trouble in explaining which particular picture we see, 
and why, that cautious scepticism seems only appropriate. Even if the premise 
were accepted, it would be essential to ask, of any particular many-picture 
scheme, whether its assumptions are natural and whether the descriptions of 
nature it produces are physically plausible or scientifically useful. The partic- 
ular equations used to define consistent sets are, after all, simply interesting 
guesses: there is no compelling theoretical justification for them, and indeed, 
several different definitions of consistency have been proposed [1,4,32]. 

My own view is that there are a number of compelling reasons for regard- 
ing the consistent histories interpretation, as it is presently understood, as 
scientifically unsatisfactory. However, as these questions have been explored 

® For example, the consistent histories interpretation of quantum mechanics has 
been referred to as “the interpretation of quantum mechanics” [16] and even as 
simply “quantum mechanics” [31]. 
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in some detail elsewhere [9,10,18-23,27,28], I here comment only on two spe- 
cific problems raised by contrary inferences. 

First, the fact that we are to take as equally valid and correct pictures of 
physics which include contrary inferences goes against many well developed 
intuitions. No argument based on intuition alone can be conclusive, but I 
think it must be granted that this one has some force. What use, it may 
reasonably be asked, is there in saying that in one picture of reality a particle 
genuinely went, with probability one, through slit A, and that in another 
picture the particle went, also with probability one, through the disjoint slit 
B? Why should we take either picture seriously, given the other? 

On this point, it is worth noting that one of the advertised merits of the 
consistent histories formalism [1,16,17] is that it, unlike the Copenhagen in- 
terpretation, accommodates some (arguably) plausible intuitions about the 
behaviour of microsystems in between observations. For example, the formal- 
ism allows us to say — albeit only as one of an infinite number of incompatible 
descriptions — that a particle observed at a particular detector was travelling 
towards that detector before the observation, and that a particle measured 
to have spin component = Sx had that spin before the measurement took 
place. As Griffiths and Omnes note [1,16], informal discussions of experiments 
are often framed in terms which, if taken literally, suggest that we can make 
this sort of statement about microsystems before a measurement is carried 
out. (“Was the beam correctly aligned going into the second interferometer?”, 
or “Do you think something crazy in the electronics might have triggered [de- 
tector] number 3 just before the particle got there?” [1], for example.) Their 
intuition is that a good interpretation of quantum theory ought to give a 
way of allowing us to take such statements literally — a criterion which, they 
suggest, the consistent histories approach satisfies. 

The intuition is, of course, controversial. A counter-intuition, which most 
interpretations of quantum theory support, is that any description of a mi- 
crosystem before a measurement is carried out should be independent of the 
result of that measurement. 

In any case, to the extent that any intuition is offered as a justification 
of the formalism, it seems reasonable to consider the fact that the formalism 
violates other strongly held intuitions. Few experimenters, after all, can ever 
have intuitively concluded that the entire flux of their beam can sensibly be 
thought of as having followed any of several macroscopically distinct paths 
through the apparatus. Yet this is what the above example, translated into 
an interference experiment, implies. 

The second, and probably deeper, problem is that it seems very hard to 
justify the distinction, which consistent historians are forced to draw, be- 
tween contradictory inferences, which are regarded as a priori unacceptable, 
and contrary inferences, which are regarded as unproblematic. Some justifi- 
cation seems called for, since the distinction is not an accidental feature: it 
is not that the formalism, for unrelated reasons, simply happens to exclude 
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one type of inference and include the other. The definition of consistency is 
motivated precisely by the notion that, when two different sets allow a calcu- 
lation of the probability of the same event (belonging, in the simplest case, 
to a single history in one set and a combination of two histories in the other), 
the calculations should agree. This requires in particular that contradictory 
propositions P and (1 — P) can never be inferred, since if the probability of P 
is one in any set, it must be one in all sets, and so the probability of (1 — P) 
must be zero in all sets. 

Now this last requirement is not absolutely essential, sensible though it 
may seem. No logical contradiction arises in an interpretation of quantum 
theory which follows the basic interpretational ideas of the consistent histo- 
ries formalism but which accepts all complete sets of disjoint quantum his- 
tories, whether consistent or not, as defining valid pictures of physics [9,10]. 
In this “inconsistent histories” interpretation, contradictory inferences can 
generally be made by using different pictures. This possibility is excluded by 
a deliberate theoretical choice. 

It seems natural, then, having made this choice, to look for ways in which 
contrary inferences can similarly be excluded. This is the line of thought pur- 
sued below. Note, however, that the problem of contrary inferences is not the 
only motivation for the ideas introduced below. Whatever one’s view of the 
consistent histories formalism, it is interesting that an alternative formalism 
can be defined relatively simply. It seems fruitful to ask which, if either, is to 
be preferred, and why. And, as we will see, ordered consistent sets raise in- 
dependently interesting questions about the quasiclassical world we actually 
observe. 

5 Relation of Contrary Inferences and Subspace 
Implications 

A contrary inference arises when there exist two consistent sets, 5i and ^ 2 , 
both containing a history H, with the property that there are orthogonal 
propositions Pi and P 2 which are implied by H in the respective sets, so that 
— temporarily adding set suffices for clarity — we have 

PsAPi\H)=PsAP2\H) = 1. (9) 

Now P52((1 ~ ^ 2 )!-^) = 0, and since the probabilities are set-independent 
and ps^{Pi\H) is nonzero, we cannot have Pi = 1 — P 2 . Hence, since Pi 
and P 2 are orthogonal, we have that Pi < 1 — P 2 . Since p{Pi\H) = 1 and 
p((l — P 2 )\H) = 0, this pair of projections violates the subspace implication 
Pi 1 — P 2 . That is, a contrary inference implies the existence of consistent 
histories H and H', belonging to different consistent sets and agreeing on 
all but one projector, such that H has non-zero probability, H' has zero 
probability, and H < H' : in the example of the last section, for instance, we 
have H = {Pb+,PJ and P' = {(1 - Ph_),Pj. 
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Clearly, according to the consistent histories formalism, an observation 
of the datum cannot be taken to imply an observation of the strictly 
larger projector (1 — Pb_). To make that inference would lead directly to a 
contradiction, in the form of the realisation of a probability zero history, if 
Pc were subsequently observed. 

Now it is easy to produce real world examples of contrary inferences, so 
long as those inferences are of unobserved quantities. As we have seen, it 
requires only a three-dimensional quantum system, prepared in one state, 
isolated, and then observed in another state — hardly a taxing experiment. 

It is not obvious, though, that we can produce examples where subspace 
implications fail in a realistic consistent histories description of observations 
of laboratory experiments, or more generally of macroscopic quasiclassical 
physics. That general consistent histories violate subspace implications need 
not necessarily imply that the particular consistent histories used to recover 
standard descriptions of real world physics do so. Both in order to address 
this question, and for its own sake, it is interesting to ask whether there 
might be any alternative treatment of quantum theory within the consistent 
histories framework which respects subspace implications, at least when they 
relate two consistent histories. The next section suggests such a treatment. 

6 Ordered Consistent Sets of Histories 

We have already seen that there is a natural partial ordering for each of the 
standard representations of quantum histories in the consistent histories ap- 
proach. The probability weight defines a second partial ordering: H ^ H' ii 
p{H) < p{H'). The violation of subspace implications reflects the disagree- 
ment between these two partial orderings in the consistent histories formal- 
ism: we can have both H < H' and H >- P[' . The aim of this section is to 
develop a history-based interpretation which restricts attention to collections 
of quantum histories on which the two orderings do not disagree. 

We begin at the level of individual quantum histories, defining an ordered 
consistent history, H, to be a consistent history with the properties that: 

(i) for all consistent histories H' with H' > H we have that p{H') > p{H); 

(ii) for all consistent histories H' with H' < H we have that p{H') < p{H). 

Recall that a consistent history is any quantum history which belongs to some 
consistent set of histories. Properties (i) and (ii) hold trivially for histories H 
and H' which belong to the same consistent set: it is the comparison across 
different sets which makes them useful constraints. 

We now define an ordered consistent set of histories to be a complete set 
of exclusive alternative histories, each of which is ordered consistent. We can 
then define an ordered consistent histories approach to quantum theory in 
precise analogy to the consistent histories approach, using the same definition 
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of probability weight and the same interpretation, simply declaring by fiat 
that only ordered consistent sets of histories are to be considered.® 

The following lemmas show that, within the projection operator formula- 
tion, ordered consistent sets of histories do indeed exist. ^ 

Lemma 1: Any consistent history H = {Pi,...,P„} defined by a 

series of projections which include a minimal projection Pj, so that Pi > Pj 
for all i, is an ordered consistent history. 

Proof: Suppose H' = {P{, . . . ,Pf} is a larger history than H, and 

write P' = Pj + Qi- We have that 

p{H')=Tr{Pf...P[pP[...Pf) 

= Tr{{Pj + Qn) ■ . . {Pj + Qi)p{Pj + Qi) . . . {Pj + Qn)) 

= Pr{Pjp) + Tr(Q„ . . . QipQi . . . Qn) 

> Tr{Pjp) 

= p{H). (10) 

Now suppose that H' is smaller than H. Then in particular Pj < Pj and we 
have that 



p{H') < Tr(Pjp) 

< Tr(P,p) 

= p{H). (11) 



Lemma 2: Let S = {ci, . . . , ct„} be a consistent set of histories defined 

by a series of projective decompositions, of which one, Oj, has the property 
that for each projection P in aj, and for every i, we have that there is precisely 
one projection Q in ai with the property that Q > P (so that all of the other 
projections in ai are contrary to P). Then S is an ordered consistent set of 
histories. 

Proof: Each of the histories of non-zero probability in S satisfies the 

conditions of Lemma 1 and so is ordered consistent. Each of the histories of 
zero probability in S is of the form H ={..., P, ... ,Q, .. .}, where P and Q 
are contrary projections. Now any consistent history smaller than H therefore 
also contains a pair of contrary projections P' < P and Q' < Q. By the 

® Note that there are other collections of quantum histories on which the orderings 
do not disagree. For example, if all the consistent histories H that violate (i) 
are eliminated, the remainder form a collection on which the orderings do not 
disagree and which is not obviously identical to the ordered consistent histories, 
and similarly for (ii). It might be interesting to explore such alternatives, but we 
restrict attention to the ordered consistent histories here. 

I am grateful to Bob Griffiths for suggesting a slight extension of Lemma 1 . 
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consistency axioms, its probability is less than or equal to Tv{Q' P' pP' Q') = 
0, and thus must also be zero. Hence PI is ordered consistent, since any 
consistent history larger than PI has probability greater than or equal to 
zero. 

7 Ordered Consistent Sets and Quasiclassicality 

A formalism based on ordered consistent sets of histories obviously defines 
a more strongly predictive version of quantum theory than that defined by 
the existing consistent histories framework, since it allows strictly fewer sets 
as possible descriptions of physics. But can it describe our empirical obser- 
vations? 

The question subdivides. Are quasiclassical domains generally ordered 
consistent sets of histories? Is our own quasiclassical domain one? If not, are 
its histories generally ordered when compared to histories belonging to other 
consistent sets defined by projections onto ranges of the same quasiclassical 
variables? For example, can we show that consistent histories defined by 
projections onto ranges of densities for chemical species in small volumes are 
generally ordered with respect to one another? If so, then the type of subspace 
implication which is generally used in analyses of observations could still be 
justified. Finally, if either of the previous two properties fail to hold, it would 
be useful to quantify the extent to which they fail. 

Answering any of these questions definitively may require — and, it might 
be hoped, help to develop — a deeper understanding of quasiclassicality than 
is available to us at present. I at any rate do not know the answers, and 
can only offer the questions as interesting ones whose resolution would have 
significant implications. At least ordered consistency does not seem to fall at 
the first hurdle: ordered consistent sets are shown below to be adequate to 
describe quasiclassicality in simple models. 

As a simple example, consider the following model of a series of successive 
measurements of the spin of a spin- 1/2 particle about various axes. We use a 
vector notation for the particle states, so that if u is a unit vector in the 
eigenstates of cr • u are represented by | ± u) . With the analogy of a pointer 
state in mind, we use the basis {| t)fe; I i)fc} to represent the environment 
particle state, together with the linear combinations |±)fc = (| t)fci| i)fc)/'\/2- 
We compactify the notation by writing environment states as single kets, so 
that for example | t)i ® ‘ ‘ ‘ ® I t)n is written as | f • ■ • t)> £^nd we take the 
initial state |V'(0)) to be |v) 0 | t ■ • • t)- 

The interaction between the system and the fc**' environment particle is 
chosen so that it corresponds to a measurement of the system spin along the 
Uk direction, so that the states evolve as follows: 



|Ufc) 0 I t)fc |Ufc) ® I t)fe , 
|-Ufc) 0 I t)fc |-Ufe) ® U)fc . 



( 12 ) 
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A simple unitary operator that generates this evolution is 

Uk{t) = P(ufe) ®h + P(-Ufc) (g) exp(-i6>fc(t)i^fc) , (13) 

where P(x) = |x)(x| and Fk = i| i)fc(t |fc-i| t)fcU U- Here 9k{t) is a function 
defined for each particle k, which varies from 0 to tt/2 and represents how 
far the interaction has progressed. We define Pk{F) = |±)fc(±|fc, so that 

Fk = Pk{+) - Pk{-)- 

The Hamiltonian for this interaction is thus 

Hk{t) = iilkmlit) = 0k{t)P{-Uk) 0 Fk , (14) 

in both the Schrodinger and Heisenberg pictures. We write the extension of 
Uk to the total Hilbert space as 

Vfe = P(ufc) (g) Ji (g) • • • (g) /„ (15) 

+P(-Ufc) (g) 7i (g) • • • (g) Jfc_i (g) exp{-i6k{t)Fk) (g) Ik+i <g) • • • <g) /« . 

We take the system particle to interact initially with particle 1 and then 
with consecutively numbered ones, and there is no interaction between en- 
vironment particles, so that the evolution operator for the complete system 
is 

U{t)=Vn{t)...Vi{t), (16) 

with each factor affecting only the Hilbert spaces of the system and one of 
the environment spins. 

We suppose, finally, that the interactions take place in disjoint time in- 
tervals and that the first interaction begins at t = 0, so that the total Hamil- 
tonian is simply 

n 

H{t) = Y^Hk{t), (17) 

k=l 

and we have that 9i{t) > 0 for t > 0 and that, if 0 < 9k{t) < tt/2, then 
9i{t) = tt/2 for alH < A: and 9i{t) = 0 for all i > k. 

This model has been used elsewhere [25,33] in order to explore algorithms 
which might select a single physically natural consistent set when the physics 
is determined by the simplest type of system-environment interaction. It is 
particularly well suited to such an analysis, since the dynamics are chosen so 
as to allow a simple and quite elegant classification [33] of all the consistent 
sets built from projections onto subspaces defined by the Schmidt decom- 
position. Apart from this, though, the model is unexceptional — one of the 
simpler variants among the many models used in the literature to investigate 
the decoherence of system states by measurement-type interactions with an 
environment. 

To give a physical interpretation of the model, we take it that the en- 
vironment “pointer” variables assume definite values after their respective 
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interactions with the system. That is, after the interaction, the fc**' en- 
vironment particle is in one of the states | t)fe and | i)fc: the probabilities of 
each of these outcomes depend on the outcome of the previous measurement 
(or, in the case of the first measurement, on the initial state) via the standard 
quantum mechanical expressions. 

This description can be recovered from the consistent histories formalism 
by choosing the consistent set 5i, defined by the decompositions 

{ / (g) |ei)i(ei|i (g) / 0 • • • (g) / : Cl = t or |} at time , ( 18 ) 

{ / (g) |ei)i(ei|i (g) |e2)2(£2|2 (g> • • • (g> -f : ei, C2 = t or |} at time ^2 , 



{ 1(g) |ei)i(ei|i (g) |e 2 ) 2 (£ 2|2 (g> • • • (g> |en)™(e„|„ : ei,C 2 , . . . ,e„ = t or |} 

at time . 



Clearly, the histories of non-zero probability in 5i take the form 



={/(g)/(g)/(g)---(g)/, 

I ® |ei)i(ei|i (g) / (g) • • • (g) J, 

I (g) |ei)i(ei|i (g) 162)2(6212 (g) / (g) • • • (g) /, 

. . . , 

I (g) |6l)l(6i|i (g) 162)2(6212 (g) • • • (g) |6„)„(6„|„} , (19) 

for sequences { 61 ,..., 6 „}, each element of which takes the value f or 
Their probabilities, defined by the decoherence functional, are precisely those 
which would be obtained from standard quantum theory by treating each 
interaction as a measurement: 



, s .l-l-a1V.U1l-l-a2U1.U2, , 1 -|- a„u„„i.Un , , , 

p(i?£i... .,£„) = ( ^ )( s )■■■( ^ ), (20) 



where, letting 60 =t> we define a^ = 1 if 6 i and 6 j_i take the same value, and 
Qi = —1 otherwise. 

Now iSi is defined by a nested sequence of increasingly refined projec- 
tive decompositions, all of whose projections commute — a relation which 
is unaltered by moving to the Heisenberg picture. It therefore satisfies the 
conditions of Lemma 2 above, and so is ordered consistent. 

This argument clearly generalizes: in any situation in which Hilbert space 
factorizes into system and environment degrees of freedom, where the self- 
interactions of the latter are negligible, any consistent set defined by nested 
commuting projections onto the environment variables is ordered consistent. 
The model considered above is a particularly crude example: more sophis- 
ticated, and phenomenologically somewhat more plausible, examples of this 
type are analysed in, for example. Refs. [26,32]. 

No sweeping conclusion can be drawn from this, since it is generally agreed 
that familiar quasiclassical physics is not well described in general — at least 
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in any obvious way — by models of this type. (Again, a detailed discussion of 
the limitations of such models can be found in Refs. [26,32].) In other words, 
while it would be hard to defend the hypothesis that familiar quasiclassical 
sets are generally ordered consistent if sets of the type 5i were not, the fact 
that they are is certainly not sufficient evidence. It would be good to find 
sharper tests of the hypothesis, perhaps for example by developing further 
the phenomenological investigations of quasiclassicality pursued in Ref. [26] . 
Meanwhile, the questions raised earlier in this section remain unresolved. 

On the other hand, it would be difficult to make a watertight case that or- 
dered consistent sets are definitely inadequate to describe real-world physics, 
for the following reason. First, it seems hard to exclude the possibility that 
the initial state is pure, so let us temporarily suppose that it is: p = [■i/')(V'l- 
As Gell-Mann and Hartle point out [32] , we can then associate to every con- 
sistent set of histories, S, a nested set of commuting projections defining what 
they term generalized records. The consistent set defines a resolution of the 
initial state into history vectors, 

iv>)= 5] p^...priV’), (21) 

which are guaranteed to be orthogonal by the consistency condition 3. We can 
thus find at least one set of orthogonal projection operators {i?/}, indexed 
by sets of the form I = {ii . . . in}, which project onto the history vectors and 
sum to the identity: 



J2Ri =I, (22) 

/ 

RiRj = SjjRj . 



We can thus [32] construct a set S' , with the same history vectors and the 
same probabilities as S, built from a nested sequence of commuting projec- 
tions defined by sums of the Rj. And, as we have seen, sets of this type are 
ordered consistent. 

There is no reason to expect the projections defining the set S' to be 
closely related to those defining S. In particular, the fact that 5 is a quasi- 
classical domain certainly does not imply that S' is likely to be: its projections 
are not generally likely to be interpretable in terms of familiar variables. But, 
as we have already noted, we have no theoretical criterion which identifies a 
particular consistent set, or a particular type of variable, as fundamentally 
correct for representing the events we observe. The set S' correctly identifies 
the history vectors and predicts their probabilities, and we thus could not 
say for certain (given that we presently have no theory of set selection) that 
its description of physics is fundamentally incorrect, while that given by S is 
fundamentally correct. 
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To make this observation is merely to point out a logical possibility. In 
fact, it would be extremely puzzling if the more complicated and apparently 
derivative set S' were in some sense more fundamental than the associated 
quasiclassical domain S. And even if this were somehow understood to be 
true in principle, we would still need to understand the relationship between 
between ordered consistency and quasiclassicality in order to say whether or 
not standardly used subspace inferences are in fact justifiable. 

8 Ordering and Ordering Violations: Interpretation 

However, one of the main points of this paper — and the main reason for 
taking a particular interest in the properties of ordered consistent sets — is 
that either answer leads to an interesting conclusion. 

If our empirical observations can be accounted for by the predictions of an 
ordered consistent set, then the ordered consistent sets formalism supersedes 
the current versions of the consistent histories formalism as a predictive the- 
ory. Alternatively, if the predictions of ordered consistent histories quantum 
theory are false, then either the consistent histories framework uses entirely 
the wrong language to describe histories of events in quantum theory, or we 
cannot generally rely on subspace implications in analysing our observations. 

Either of these last two possibilities would have far-reaching implications 
for our understanding of nature. It is true that other ways of representing 
quantum histories are known than those used in the consistent histories for- 
malism, but they arise either in non-standard versions of quantum theory, 
such as de Broglie-Bohm theory, or in alternative theories. It is also true 
that there is no way of logically excluding the possibility that subspace im- 
plications generally fail to hold. Any clear violation would, however, lead to 
radical changes in our representation of the world, and in particular to our 
understanding of the relation between theory and empirical observation. 

I would suggest, however, that any version of closed system quantum the- 
ory in which the two orderings disagree leads to radical new interpret ational 
problems. The fundamental problem is that, supposing that the world we 
experience is described by one particular realised quantum history, we never 
know — no matter how precise we try to make our observations — exactly 
what form that history takes. This is not only because we can never com- 
pletely eliminate imprecision from our experimental observations. A deeper 
problem is that we have no theoretical understanding of how, precisely, an 
observation should be represented within quantum theory. We do not know 
precisely when and where any given observation takes place. Nor do we know 
whether is fundamentally correct to represent quantum events by projection 
operators, by quantum effects, by statements associated to space-time regions 
in path integral quantum theory, or in some other way — let alone precisely 
which operator, effect, or statement correctly represents any given event. As a 
result, we are always forced into guesswork and approximation. We are forced 
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to assume, at least as a working hypothesis, that we can find sensible bounds 
on our observations. Roughly speaking, we assume that we can say, at least, 
that a photon hit our photographic plate within a certain region, that the 
observed flux from a distant star was in a certain range, and so forth.® We 
assume also that the probability of the actual observations — whose precise 
form we do not know — is bounded by the probability of the observations as 
we approximately represent them. These assumptions ultimately rely on the 
agreement of two orderings just mentioned: when those orderings disagree, 
we therefore run into new problems. 

It is easy to see, in particular, that this sort of problem arises in any 
careful consistent histories treatment of quantum cosmology. Suppose, for 
example, that we have a sequence of cosmological events which we wish to 
represent theoretically, in order to calculate their probability, given some 
theory of the boundary conditions. Assuming that the basic principles of the 
consistent histories formalism are correct, we know that these events should 
be represented by some history H belonging to some consistent set S. We do 
not, however, know the precise form oi H or of 5: the events are given to us 
as empirical observations rather than as mathematical constructs. 

The best we can then do, following the general principles of the consistent 
histories formalism, is to choose some plausible consistent set S' containing 
histories i/min and i/max which we guess to have the property i?min < H < 
Hmax- in particular, thus, we choose i?min < ^^max- Since i?min and i/max 
belong to the same set, we have that p(i?min) < p(^fmax)- It might naively be 
hoped that we can derive that p(i?min) < p{H) < p(i?max), but since H in 
general will belong to a different consistent set from iJmin and this does 

not generally follow. There is no way to bound p{H), except (in principle) 
by performing the enormous task of explicitly calculating the probabilities 
of all consistent histories bounded by i?min and i?max, and there is no way 
to justify the type of subspace implication — relating observations and true 
data — that we generally take for granted. 

This is not to say that the disagreement of the two orderings necessarily 
leads to logical contradiction. Versions of quantum theory in which the or- 
derings disagree need not be inconsistent, or even impossible to test precisely. 
They do, though, generally seem to require us to identify precisely the correct 
representation of our observations in quantum theory. This is generally a far 
from trivial problem: how are we to tell, a priori, exactly which projection 
operators represent the results of a series of quantum measurements? It is 
not impossible to imagine that theoretical criteria could be found which solve 
the problem, but we certainly do not have such criteria at present. 



In fact such statements are generally made within statistical conhdence limits. 
To consider statistical statements would complicate the discussion a little, but 
does not alter the underlying point. 
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9 Conclusions 

Though the criterion of ordering seems mathematically natural, both in the 
consistent histories approach to quantum theory and in other possible treat- 
ments of quantum histories, it raises very unconventional questions. It seems, 
though, that these questions cannot be avoided in any precise formulation of 
the quantum theory of a closed system which involves a standard represen- 
tation of quantum events and which gives a historical account. 

There seem to be three possibilities, each of which is interesting. The first 
is that the representations of quantum histories discussed here, though stan- 
dard, are not those chosen by nature. Clearly this is a possibility: there are, 
for example, well known non-standard versions of quantum theory [6], and 
related theories [7], in which histories are defined by trajectories or other aux- 
iliary variables, and in which subspace implications follow just as in classical 
physics. 

The second possibility is that our quasiclassical domain can be shown 
to be an ordered consistent set. If so, then the ordered consistent histories 
approach is both predictively stronger than the standard consistent histories 
approach — since there are fewer ordered consistent sets — and compatible 
with empirical observation, and hence superior as a predictive theory. If it is 
compatible with our observations, the ordered consistent histories approach 
would seem at least as natural as the consistent histories approach. 

Even if so, I would not suggest that the ordered consistent histories formal- 
ism is the “right” interpretation of quantum theory, and the consistent his- 
tories approach the “wrong” one. The ordered consistent histories approach 
seems almost certain to suffer from many of the same defects as the consis- 
tent histories approach, since there are still far too many ordered consistent 
sets. The aim here is thus not to propose the ordered consistent histories 
approach as a plausible fundamental interpretation of quantum theory, but 
to suggest that the range of natural and useful mathematical definitions of 
types of quantum history is wider than previously understood. This range 
includes, at least, Goldstein and Page’s criterion of linear positivity [34], the 
various consistency criteria [1,26,32] in the literature, and the criterion of 
ordered consistency introduced here. It seems to me hard to justify taking 
any of these criteria as defining the fundamentally correct interpretation of 
quantum theory. On the one hand, physically interesting quantum histories 
might possibly satisfy any one, or none, of them; on the other hand, most 
quantum histories satisfying any given criterion seem unlikely to be physically 
interesting — and precisely which criteria are useful in which circumstances 
largely remains to be understood. 

The third possibility is that our quasiclassical domain is not an ordered 
consistent set. This would have intriguing theoretical implications. We would 
have, at least in principle, to abandon subspace inferences, and we would 
ultimately need to understand precisely how to characterise the quantum 
events which constitute the history we observe. This would raise profound and 




Quantum Histories and Their Implications 



111 



not easily answerable questions about how we can tell what, precisely, are our 
empirical observations. It might also, depending on the way in which ordered 
consistency was violated, and the extent of any violation, raise significant 
practical problems in the analysis of those observations. 

No compelling argument in favour of any one of these possibilities has 
been given here: it has been shown only that, if quasiclassical sets generally 
fail to be ordered consistent, they do so in a way too subtle to be displayed 
in the simplest models. 

Another caveat is that the above discussion applied the criterion of or- 
dered consistency only to the simplest representation of quantum histories, 
in which individual events are represented by projections at a single time. 
Other representations need to be considered case by case, and our conclu- 
sions might not necessarily generalize. For example, the fact that a consistent 
history built from single time projections is ordered when compared to con- 
sistent histories of the same type does not necessarily imply that it is ordered 
when compared to consistent histories defined by composite events. 

Still, the criterion of ordered consistency defines a new version of the 
consistent histories formulation of quantum theory, which avoids the problems 
caused by contrary inferences. Its other properties and implications largely 
remain to be understood. 
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Appendix: Ordering and Decoherence Functionals 

This Appendix describes a noteworthy earlier discussion of quantum history 
orderings, given by Isham and Linden in Sec. IV of Ref. [14], and its relation 
to the ideas discussed here. 

Isham and Linden abstract the basic ideas of the consistent histories for- 
malism in the following way. First, the space UV of history propositions is 
taken to be a mathematical structure — an orthoalgebra — with a series of 
operations and relations obeying certain axioms. In particular, they propose 
that a partial ordering < and an orthogonality relation T should be defined 
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on UV and should obey natural rules, and that UV should include an identity 
history 1 . 

They then introduce a space T> of decoherence functionals, defined to be 
maps from UV x UV to the complex numbers satisfying certain axioms, and 
go on to consider whether the axioms defining decoherence functionals should 
include axioms relating to the ordering in UV. 

In the language of standard quantum mechanics, UV corresponds to the 
space of all the quantum histories (not only the consistent histories) for a 
given system, whose Hilbert space and hamiltonian are fixed. Any of several 
representations of quantum histories could be considered: the relevant part of 
Isham and Linden’s discussion uses the simplest representation of quantum 
histories, as sequences of projection operators. 

The standard quantum mechanical decoherence functional (as appears 
on the left hand side of (3)) is a member of the space T> in the minimal 
axiom system Isham and Linden eventually choose. As they remark, though, 
it would not be a member of T> if the extra ordering axioms they discuss were 
imposed. Isham and Linden nonetheless consider imposing these ordering 
axioms, since their aim in the relevant discussion is to investigate generalised 
algebraic and logical schemes rather than to propose a formalism applicable 
to standard quantum theory. (They suggest, at the end of section IV, that 
standard quantum theory might perhaps emerge from some such generalised 
scheme in an appropriate limit.) 

Isham and Linden were, as far as I am aware, the first to investigate pos- 
sible uses of orderings in developing the consistent histories formalism. It is 
worth stressing, though, to avoid any possible confusion, that their sugges- 
tions pursue the exploration of orderings in a direction orthogonal to the one 
considered in the present paper. In this paper we restrict attention to stan- 
dard quantum theory, and propose an alternative histories formalism within 
that theory, using the standard quantum theoretic decoherence functional 
throughout. We note also that the subspace implications which underlie our 
basic scientific worldviews depend for their justification on the assumption 
that the quasiclassical set describing the physics we observe is an ordered con- 
sistent set. Isham and Linden’s proposed ordering axioms, on the other hand, 
exclude standard quantum theory and the standard decoherence functional: 
they are possible postulates which might be imposed on non-standard gen- 
eralised decoherence functionals in non-standard generalisations of quantum 
theory. 

Isham and Linden give a minimal set of postulated properties for gener- 
alised decoherence functionals: 

d(0, a) = 0 for all a G UV ; 
d\a, (3) = d{(3, a)* for all a, /3 G UV ; 
d{a, a) > 0 for all a G UV ; 

if a T /3 then, for all 7 , d{a © /3, 7 ) = d{a, 7 ) + d{(3, 7 ) ; 
d(l,l) = 1 . 



( 23 ) 
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They then consider imposing new postulates on decoherence functionals. 
The first of these — their posited inequality 1 — is that: 

for all d G P and for all a, /3 with a < /3 we have d{a, a) < d{(3, (3) . (24) 

As Isham and Linden go on to point out, there are familiar examples 
in standard quantum theory in which (24) is violated for a pair of histories 
a < P in which one of the histories (in their case a) is inconsistent.® 

Two further postulates on generalised decoherence functionals are also 
posited: 

a ± P implies d{a, a) + d{P, P) < I for all d G 22 ; (25) 

and 

for all d G 22 and all 7 G UV we have d(y, 7 ) < 1 . (26) 

Isham and Linden give examples to show that, in standard quantum the- 
ory, with the standard decoherence functional, inconsistent histories do not 
necessarily respect these inequalities either. 

Again, the difference from the examples considered in the present paper 
is worth emphasizing. All of the examples Isham and Linden consider involve 
inconsistent histories — these are all they require in order to investigate possi- 
ble properties of decoherence functionals applied to arbitrary, not necessarily 
consistent, quantum histories. These examples are not problematic for the 
consistent histories approach to quantum theory, according to which the in- 
consistent histories have no physical significance, and they do not give rise 
to new interpretational questions in any conventional quantum histories ap- 
proach, for essentially the same reason. The discussion in the present paper, 
on the other hand, looks at the properties of consistent histories in stan- 
dard quantum theory: we have argued that their failure to respect ordering 
relations is problematic and explained that it does raise new questions. 

Suppose now that we set aside Isham and Linden’s motivations, and alter 
their ordering postulates so that they apply, not to generalised decoherence 
functionals applied to all quantum histories in an abstract generalisation 
of quantum theory, but to the standard decoherence functional applied to 
ordered consistent histories in standard quantum theory. We then obtain the 
following: 



for all a, P with a < P we have d{a, a) < d{P, P ) ; (27) 

a 3- P implies d{a, a) + d{P, P) < 1; (28) 

and 

for all 7 we have ^( 7 , 7 ) < 1 . (29) 

® The suggestion that inequality 1 is true when applied to sequences of projectors 
onto subsets of configuration space in a path-integral quantum theory is thus 
misleading: it is easy to find confignration space analogues of these examples. I 
am grateful to Chris Isham and Noah Linden for discnssions of this point. 
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Here d is the standard decoherence functional, a, l3 and 7 are now taken to 
be ordered consistent histories, and a _L /3 means that a and [3 are disjoint — 
i.e., there is at least one time at which their respective events are represented 
by contrary projections. 

The first of these equations holds by the definition of an ordered consistent 
history, but it might perhaps be hoped that the others could restrict the class 
of histories further. However, the second equation also holds for all ordered 
consistent histories. To see this, note that a T /3 implies that j3 < (1 — a), 
and that if a is a consistent history then (1 — a) is too. The fact that (3 is 
ordered consistent thus implies that 

p{!3) < p{l - a) = I - p{a) . (30) 

The third equation, moreover, holds for all consistent histories, ordered or 
otherwise. It seems, then, that ordered consistency may be the strongest 
natural criterion that can be defined using the basic ingredients of consistency 
and ordering. 
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Abstract. We discuss the role of non-locality in the problem of determining the 
state of a quantum system, one of the most basic problems in quantum mechanics. 



1 Introduction 

In 1964 J. Bell [1] introduced the idea of non-local correlations. He consid- 
ered the situation in which measurements are performed on pairs of particles 
which were prepared in entangled states and the two members of the pair 
are separated in space, and he showed that the results of the measurement 
are correlated in a way which cannot be explained by any local models. Since 
then, non-locality has been recognized as one of the most important aspects 
of quantum mechanics. 

Following Bell, most studies of non-locality focused on situations when 
two or more particles are separated in space and are prepared in entangled 
states. Entanglement appeared to be a sine-qua-non requirement for non- 
locality. Indeed, all entangled states were shown to be non-local [2] while all 
direct product states lead to purely local correlations. 

Surprisingly enough, however, it turns out that non-locality plays an es- 
sential role in problems in which there is no entanglement whatsoever and 
apparently everything is local. One such problem - probably one of the most 
basic problems in quantum mechanics - is that of determining by measure- 
ments the state of a quantum system. That is, given a quantum system in 
some state 'P unknown to us, how can we determine what S' is? As is well- 
known, because quantum measurements yield probabilistic outcomes, and 
due to the fact that measuring an observable we disturb all the other non- 
commuting observables, we cannot determine with certainty the unknown 
state iF out of measurements on a single particle. If we have a large number 
of identically prepared particles - a so called “quantum ensemble” - by per- 
forming measurements on the different particles we can accumulate enough 
results so that from their statistics we could determine <F. But anything less 
than an infinite number of particles limits our ability to determine the state 
with certainty. 

Now, although when we have only a finite ensemble of identically pre- 
pared particles we cannot determine their state with absolute precision, the 
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question remains of how well can we determine the state and which is the 
the measurement we should perform. It is in this context that, unexpectedly, 
non-locality creeps in. 

The story of interest here begins in 1991 when Peres and Wootters [3] 
considered the following particular question of quantum state estimation. 
Suppose that we have a quantum ensemble which consists of only two iden- 
tically prepared spin 1/2 particles. In other words, we are given two spin 1/2 
particles, and we are told that the spins are parallel, but we are not told 
the polarization direction. What is the best way to find the direction? One 
possibility would involve measurements carried out on each spin separately, 
and trying to infer the polarization direction based on the results of the two 
measurements. Another possibility would be again to measure the two spins 
separately, but to measure first one spin and then choose some appropriate 
measurement on the second spin depending on the result of the measurement 
on the first spin. Finally, they conjectured, an even better way might be to 
measure the two spins together, not separately. That is, to perform a mea- 
surement of an operator whose eigenstates are entangled states of the two 
spins. The Peres- Wootters conjecture was subsequently proven by Massar 
and Popescu [4]. 

What we see here is a surprising manifestation of non-locality. Indeed, 
suppose that the two spins are separate in space. They are prepared in direct- 
product states (they are each polarized along the same direction) . Neverthe- 
less, by simply performing local measurements one cannot extract all infor- 
mation from the spins. The optimal measurement must introduce non-locality 
by projecting them onto entangled states. 

The above effect has opened a new direction of research, and its implica- 
tions are just starting to be uncovered. In this paper we will discuss a new 
surprising effect [5] whose existence is made possible by the fact that optimal 
measurements imply non-locality. 

2 Measurements on 2-Particle Systems 
with Parallel or Anti-Parallel Spins 

The problem we consider here is the following. Suppose Alice wants to com- 
municate Bob a space direction n. She may do that in two ways. In the first 
case, Alice sends Bob two spin 1/2 particles polarized along n, i.e. \n,n). 
When Bob receives the spins, he performs some measurement on them and 
then guesses a direction rig which has to be as close as possible to the true 
direction n. The second method is almost identical to the first, with the dif- 
ference that Alice sends |n, — n), i.e. the first spin is polarized along n but 
the second one is polarized in the opposite direction. The question is whether 
these two methods are equally good or, if not, which is better^. 

^ It might be the case that there is some other method which is better for trans- 
mitting directions than the two ones mentioned above. However, we are not 
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For a better perspective, consider first a simpler problem. Suppose Alice 
wants to communicate Bob a space direction n and she may do that by one 
of the following two strategies. In the first case, Alice sends Bob a single spin 
1/2 particle polarized along n, i.e. \n). The second strategy is identical to 
the first, with the difference that when Alice wants to communicate Bob the 
direction n she sends him a single spin 1/2 particle polarized in the opposite 
direction, i.e. \ — n). Which of these two strategies is better? 

Obviously, if the particles would be classical spins then, both methods 
would be equally good, as an arrow defines equally well both the direction 
in which it points and the opposite direction. Is the quantum situation the 
same? 

First, we should note that in general, by sending a single spin 1/2 particle, 
Alice cannot communicate to Bob the direction n with absolute precision. 
Nevertheless, it is still obviously true that the two strategies are equally good. 
Indeed, all Bob has to do in the second case is to perform exactly the same 
measurements as he would do in the first case, only that when his results are 
such that in the first case he would guess rig, in the second case he guesses 

-rig. 

One is thus tempted to think that, similar to the classical case, for the 
purpose of defining a direction n, a quantum mechanical spin polarized along 
n is as good as a spin polarized in the opposite direction: in particular, the 
two two-spin states \n,rn) and \n,—rn) should be equally good. Surprisingly 
however this is not true. 

That there could be any difference between communicating a direction by 
two parallel spins or two anti-parallel spins seems, at first sight, extremely 
surprising. After all, by simply flipping one of the spins we could change 
one case into the other. For example, if Bob knows that Alice indicates the 
direction by two anti-parallel spins he only has to flip the second spin and 
then apply all the measurements as in the case in which Alice sends from the 
beginning parallel spins. Thus, apparently, the two methods are bound to be 
equally good. 

The problem is that one cannot flip a spin of unknown polarization. In- 
deed, the flip operator V defined as 

V\n) = \-n) (I) 

is not unitary but anti-unitary. To see this we note that in Heisenberg rep- 
resentation flipping the spin means changing the sign of all spin operators, 
cr ^ — (T. But this transformation does not preserve the commutation rela- 
tions, i.e. 

[<^x,CTy\ = iCTz -)> [(Jx,<7y\ = -icT^, (2) 



interested here in finding an optimal method for communicating directions; we 
are only interested in comparing the parallel and anti-parallel spins methods. 
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while all unitary transformations leave the commutation relations unchanged. 
Thus there is no physical operation which could implement such a transfor- 
mation. 

A couple of questions arise. First, why is it still the case that a single spin 
polarized along n defines the direction as well as a single spin polarized in 
the opposite direction? The reason is that although Bob cannot implement 
an active transformation, i.e. cannot flip the spin, he can implement a passive 
transformation: he can flip his measuring devices. Indeed, there is no problem 
for Bob in flipping all his Stern-Gerlach apparatuses, or, even simpler than 
that, to merely rename the outputs of each Stern-Gerlach “up”— “down” 
and “down” — >■ “up” . 

At this point it is natural to ask why can’t Bob solve the problem of 
two spins in the same way, namely by performing a passive transformation 
on the apparatuses used to measure the second spin? The problem is the 
entanglement. Indeed, if the optimal strategy for finding the polarization 
direction would involve separate measurements on the two spins then two 
parallel spins would be equivalent to two anti-parallel spins. (This would be 
true even if which measurement is to be performed on the second spin depends 
on the result of the measurement on the first spin.) But, as explained in 
the Introduction, the optimal measurement is not a measurement performed 
separately on the two spins but a measurement which acts on both spins 
simultaneously, that is, the measurement of an operator whose eigenstates 
are entangled states of the two spins. For such a measurement there is no 
way of associating different parts of the measuring device with the different 
spins, and thus there is no way to make a passive flip associated to the second 
spin. Gonsequently there is no way, neither active nor passive to implement 
an equivalence between the parallel and anti-parallel spin cases. 

After understanding that there is indeed room for the two direction com- 
munication methods to be different, let us now investigate them in detail. 

To start with, we have to define some figure of merit which will tell us 
how successful a communication protocol is. For concreteness, let us define 
Bob’s measure of success as the fidelity 



( 3 ) 

where nUg is the scalar product in between the true and the guessed direc- 
tions, the integral is over the different directions n and dn represents the 
a priori probability that a state associated to the direction n, i.e. \n, n) or 
|n, — n) respectively, is emitted by the source; P{g\n) is the probability of 
guessing Ug when the true direction is n. In other words, for each trial Bob 
gets a score which is a (linear) function of the scalar product between the true 
and the guessed direction, and the final score is the average of the individual 



scores. 
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When the different directions n are randomly and uniformly distributed 
over the unit sphere, an optimal measurement for pairs of parallel spins ip = 
\n,n) has been found by Massar and Popescu [4]. Bob has to measure an 
operator A whose eigenvectors 4>j, j = 1, . . . , 4, are 

( 4 ) 



where jf/' ) denotes the singlet state and the Bloch vectors rij point to the 4 
vertices of the tetrahedron: 



ni = (0,0, 1), 



n ^ 

«2 = ( — , 0 .--), 



«3 = ( 



-V2 



3’ 



«4 = ( 



-V2 



3’ 3'^' 



( 5 ) 



The phases used in the definition of \rij) are such that the 4 states (j)j are 
mutually orthogonal. The exact values of the eigenvalues corresponding to 
the above eigenvectors are irrelevant; all that is important is that they are 
different from each other, so that each eigenvector can be unambiguously 
associated to a different outcome of the measurement. If the measurement 
results corresponds to 4>j, then the guessed direction is rij. The corresponding 
optimal fidelity is 3/4 [4]. 

A related case is when the directions n are a priori on the vertices of 
the tetrahedron, with equal probability 1/4. Then the above measurement 
provides a fidelity of 5/6« 0.833, conjectured to be optimal. 

Let us now consider pairs of anti-parallel spins, \tp >= \n, —n), and the 
measurement whose eigenstates are 



6j = a\nj,-rij) - |nfc,-nfc) 



(6) 



with a = ^ 1-095 and (3 = ~ 0.129. The corresponding 

fidelity for uniformly distributed n is F = ~ 0.789; and for n lying 

on the tetrahedron F = ~ 0.955. In both cases the fidelity obtained 

for pairs of anti-parallel spins is larger than for pairs of parallel spins! 

It is useful now to investigate in more detail what is going on. We have 
claimed above that when we perform a measurement of an operator whose 
eigenstates are entangled states of the two spins, there is no way of making a 
passive flip associated with the second spin. We would like to comment more 
about this point. 

It is clear that in the case of a measuring device corresponding to an 
operator whose eigenstates are entangled states of the two spins, we cannot 
identify one part of the apparatus as acting solely on one spin and another 
part of the apparatus as acting on the second spin. Thus we cannot simply 
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isolate a part of the measuring device and rename its outcomes. But perhaps 
one could make such a passive transformation at a mathematical level, that 
is, in the mathematical description of the operator associated to the mea- 
surement and then physically construct an apparatus which corresponds to 
the new operator. 

In the case of two parallel spins the optimal measurement is described 
by a nondegenerate operator whose eigenstates are given by (4) and 
(5). It is convenient to consider the projectors = \4’j){4’j \ associated with 
the eigenstates. As is well-known, any unit-trace hermitian operator, and in 
particular any projector, can be written as 

po = 1 (/ + )). (7) 

with some appropriate coefficients , (3^ and R\i- (The upper indexes on 
the spin operators mean “particle 1” or “2” and I denotes the identity). 
Why then couldn’t we simply make the passive spin flip by considering a 
measurement described by the projectors 

= i(/ + (8) 

obtained by the flip of the operators associated second spin, — >• — 

The reason is that the transformed operators are no longer projectors! 
Indeed, each projection operator P^ could also be viewed as a density matrix 
/A = P^ = The passive spin flip which leads from Eq. (7) to Eq. (8) 

is nothing more than the partial transpose of the density matrices fP with re- 
spect to the second spin. But each density matrix p> is non-separahle (because 
it describes the entangled state \4>j))- However, according to the well-known 
result of the Horodeckis [6,7] the partial transpose of a non-separable density 
matrix of two spin 1/2 particles has a negative eigenvalue and thus it cannot 
represent a projector anymore. On the other hand, if the optimal measure- 
ment would have consisted of independent measurements on the two spins, 
each projector would have been a direct product density matrix and the spin 
flip would have transformed them into new projectors, and thus led to a valid 
new measurement. 

The above analysis of encoding directions by parallel or anti-parallel spins 
shows a most important aspect of the problem of state estimation. Consider 
the two sets of states, that of parallel spins and that of anti-parallel spins. 
The distance in between any two states in the first set is equal to the distance 
in between the corresponding pair of states in the second set. That is, 

|(n, njm, m)p = |(n, — njm, — m) p. (9) 

Nevertheless, as a whole, the anti-parallel spin states are further apart than 
the parallel ones! Indeed, the anti-parallel spin states span the entire 4- 
dimensional Hilbert space of the two spin 1/2, while the parallel spin states 
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span only the 3-dimensional subspace of symmetric states. This is similar to 
a 3-spin example discovered by R. Jozsa and J. Schlienz [8]. 

Furthermore, suppose we consider a simpler problem in which Alice has 
to communicate Bob one out of only two possible directions, n and m. Then, 
since \ {n, n|m, m)p = |(n, — n|m, — m)p the parallel and anti-parallel spins 
methods would be equally good. The two methods are also equally good in the 
case when Alice has to communicate to Bob an arbitrary direction in a plane. 
Indeed, suppose that the directions Alice has to communicate are restricted to 
the x—y plane. Then a rotation of the second spin by 180 degrees around the z 
axis can transform any state of parallel spins into a state of anti-parallel ones 
and vice-versa. It is only when Alice has to communicate directions which 
do not lie in the same plane, that the two methods become different. This 
shows that the problem of state estimation depends on the global structure 
of the set of states under investigation and cannot be reduced to the problem 
of pairwise distinguishibility.^ 

3 Conclusions 

To conclude, we have shown that non-locality plays a fundamental role in 
what is probably the most basic quantum mechanical problem - determining 
the state of a quantum system. The link between non-locality and measure- 
ment theory is completely unexpected, and without any doubt, it will lead to 
new insights into the very nature of quantum mechanics. Furthermore, it is 
already leading to possible practical applications in quantum communication 
[9]. 
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Abstract. The loss of coherence of a quantum system coupled to a heat bath as 
expressed by the reduced density matrix is shown to lead to the mis-characterization 
of some systems as being incoherent when they are not. The spin boson problem and 
the harmonic oscillator with massive scalar field heat baths are given as examples 
of reduced incoherent density matrices which nevertheless still represent perfectly 
coherent systems. 



1 Massive Field Heat Bath and a Two Level System 

How does an environment affect the quantum nature of a system? The stan- 
dard technique is to look at the reduced density matrix, in which one has 
traced out the environment variables. If this changes from a pure state to a 
mixed state (entropy — Trpln p not equal to zero) one argues that the system 
has lost quantum coherence, and quantum interference effects are suppressed. 
However this criterion is too strong. There are couplings to the environment 
which are such that this reduced density matrix has a high entropy, while the 
system alone retains virtually all of its original quantum coherence in certain 
experiments. 

The key idea is that the external environment can be different for different 
states of the system. There is a strong correlation between the system and the 
environment. As usual, such correlations lead to decoherence in the reduced 
density matrix. However, the environment in these cases is actually tied to the 
system, and is adiabatically dragged along by the system. Thus although the 
state of the environment is different for the two states, one can manipulate 
the system alone so as to cause these apparently incoherent states to interfere 
with each other. One simply causes a sufficiently slow change in the system 
so as to drag the environment variables into common states so the quantum 
interference of the system can again manifest itself. 

An example is if one looks at an electron with its attached electromagnetic 
field. Consider the electron at two different positions. The static Coulomb 
field of the two charges differ, and thus the states of the electromagnetic field 
differ with the electron in the two positions. These differences can be sufficient 
to cause the reduced electron wave function loose coherence for a state which 
is a coherent sum of states located at these two positions. However, if one 
causes the system to evolve so as to cause the electron in those two positions 
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to come together (e.g., by having a force field such that the electron in both 
positions to be brought together at some central point for example), those 
two apparently incoherent states will interfere, demonstrating that the loss 
of coherence was not real. 

Another example is light propagating through a slab of glass. If one simply 
looks at the electromagnetic field, and traces out over the states of the atoms 
in the glass, the light beams traveling through two separate regions of the 
glass will clearly decohere- the reduced density matrix for the electromagnetic 
field will lose coherence in position space- but those two beams of light will 
also clearly interfere when they exit the glass or even when they are within 
the glass. 

The above is not to be taken as proof, but as a motivation for the further 
investigation of the problem. The primary example I will take will be of a 
spin-i particle (or other two- level system). I will also examine a harmonic 
oscillator as the system of interest. In both cases, the heat bath will be a 
massive one dimensional scalar field. This heat bath is of the general Caldeira- 
Leggett type [1] (and in fact is entirely equivalent to that model in general). 
The mass of the scalar field will be taken to be larger than the inverse time 
scale of the dynamical behaviour of the system. This is not to be taken as 
an attempt to model some real heat bath, but to display the phenomenon in 
its clearest form. Realistic heat baths will in general also have low frequency 
excitations which will introduce other phenomena like damping and genuine 
loss of coherence into the problem. 



2 Spin- 1 System 



Let us take as our first example that of a spin-^ system coupled to an external 
environment. We will take this external environment to be a one-dimensional 
massive scalar field. The coupling to the spin system will be via purely the 
3-component of the spin. I will use the velocity coupling which I have used 
elsewhere as a simple example of an environment (which for a massless field 
is completely equivalent to the Caldeira-Leggett model). The Lagrangian is 



21 



which gives the Hamiltonian 
1 



(^(x))^ — 4>'{xY — + ‘2€<p{x)h{x)a3 



da;, 



H = 



{tt{x) — eh{x)a3)^ + (j)'{x)‘^ + m^(j){xY da;. 



( 1 ) 

(2) 



h{x) is the interaction range function, and its Fourier transform is related to 
the spectral response function of Leggett and Caldeira. 

This system is easily solvable. I will look at this system in the following 
way. Start initially with the field in its free (e = 0) vacuum state, and the 
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system is in the +1 eigenstate of cti. I will start with the coupling e initially 
zero and gradually increase it to some large value. I will look at the reduced 
density matrix for the system, and show that it reduces to one which is almost 
the identity matrix (the maximally incoherent density matrix) for strong 
coupling. Now I let e slowly drop to zero again. At the end of the procedure, 
the state of the system will again be found to be in the original eigenstate 
of (Ti. The intermediate maximally incoherent density matrix would seem 
to imply that the system no longer has any quantum coherence. However, 
this lack of coherence is illusionary. Slowly decoupling the system from the 
environment should in the usual course simply maintain the incoherence of 
the system. Yet here, as if by magic, an almost completely incoherent density 
matrix magically becomes coherent when the system is decoupled from the 
environment. 

In analyzing the system, I will look at the states of the field corresponding 
to the two possible eigenstates of the system. These two states of the field 
are almost orthogonal for strong coupling. However they correspond to fields 
tightly bound to the spin system. As the coupling is reduced, the two states 
of the field adiabatically come closer and closer together until finally they 
coincide when e is again zero. The two states of the environment are now the 
same, there is no correlation between the environment and the system, and 
the system regains its coherence. 

The density matrix for the spin system can always be written as 

pW = ^(1 + pW'O-) (3) 

where 

p(t) = Tr(crp(t)). (4) 

We have 

p(t) = Tr(^^r[e-‘/o"‘i‘] l(l + p(0).^)i?or[e-‘/^‘i‘]^^ , (5) 

where i?o is the initial density matrix for the field (assumed to be the vac- 
uum), and 'T is the time-ordering operator. Because e and thus H is time- 
dependent, the iL’s at different times do not commute. This leads to the re- 
quirement for the time-ordering in the expression. As usual, the time-ordered 
integral is the way of writing the time ordered product = 

Let us first calculate ps(t). We have 

P3 (t) = Tr ^(JaT e“*-^o^di ^ ^ 

= Tr + p(0) • ^)i?or[e-‘/^'^*]t^ 
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= Tr ^<73^(1 + p(0) • ( t ) Ro ^ 

= P3(0) (6) 

because commutes with H{t) for all t. We now define 

<^+ = \{<^i + i<^2) = \+){-\, <J_=a\. (7) 



Using (T+CT3 = — CT+ and CT3CT+ = (t+ we have 



Tr (^a+r 

= Tr. fr 









it 



i j {HQ — e{i) f 7r{x)h{x)dx)dt 



lt 



(8) 



r 



^ — i J f 7r(x)h(x)dx)dt 



H-(i + p(o).^)|+) 



— (pi(0) + *P2(0)) J(t), 

where Hq is the Hamiltonian with e = 0, i.e., the free Hamiltonian for the 
scalar field and 



J{t) = 

Tr. fr 



i f (Ho — €(t) f 7r(x)h(x)dx)dt 



lt 



r 



i f (Ho+e(t) J 7r{x)h{x)dx)dt 




Breaking up the time ordered product in the standard way into a large num- 
ber of small time steps, using the fact that exp[— ie(t) f h(x)^(x)dxj is the 
displacement operator for the field momentum through a distance of e{t)h{x), 
and commuting the free field Hamiltonian terms through, this can be written 
as 



J{t) = Tr^ 



( t/dt 
n^l 



t/dt 

n—1 



gi£(0)<f(0)^ 




(10) 



where = ndt and dt is a very small value, <?(t) = f h{x)(j){t,x)dx and 
(po{t, x) is the free field Heisenberg field operator. Using the Campbell-Baker- 
Hausdorff formula, realizing that the commutators of the <l>s are c-numbers, 
and noticing that these c-numbers cancel between the two products, we finally 
get 



J(t) = Tr 0 



( 11 ) 
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from which we get 



ln( J(t)) = -2Tr^ fjio (^e(i)m ~ <0)<P(0) + ^ i(t')^(t')dt'^ ) . (12) 



I will assume that e(0) = 0, and that i(t) is very small, and that it can be 
neglected. (The neglected terms are of the form 

J J e^i$(t')^t''))dt'dt" « ehT{<P(0)^) 

which for a massive scalar field has the coherence time scale t « 1/m. Thus, 
as we let e go to zero these terms go to zero.) 

We finally have 

ln(J(t)) = -2e(t)2 < <P(t)^ > 

= -2e(t)^ I \h(k)\^-^=^dk. (13) 

Choosing h(k) = we finally get 

ln(J(t)) = -4 / e(t)^ (14) 

JO y'(k^ + m^) 

This goes roughly as In(Tm) for small Fm, (which I will assume is true). For 
r sufficiently small, this makes J very small, and the density matrix reduces 
to essentially diagonal form (pz{t) « Py(t) « 0, pz(t) = Pz(0).) 

However it is clear that if e(t) is now lowered slowly to zero, the decoher- 
ence factor J goes back to unity, since it depends only on e(t). The density 
matrix now has exactly its initial form again. The loss of coherence at the 
intermediate times was illusionary. By decoupling the system from the envi- 
ronment after the coherence had been lost, the coherence is restored. This is 
in contrast with the naive expectation in which the loss of coherence comes 
about because of the correlations between the system and the environment. 
Decoupling the system from the environment should not in itself destroy that 
correlation, and should not reestablish the coherence. 

The above approach, while giving the correct results, is not very trans- 
parent in explaining what is happening. Let us therefore take a different ap- 
proach. Let us solve the Heisenberg equations of motion for the field <j>(t,x). 
The equations are (after eliminating tt) 

df(j)(t, x) — d^(f)(t, x) + m?'(j)(t, x) = —k(t)azh(x), (15) 

’n(t,x) = (j){t,x) + e(t)h(x)a^. (16) 

If e is slowly varying in time, we can solve this approximately by 

(j){t,x) = 4>o(t,x) + k(t) y l/i(a;')da;'(J 3 -I- V’(i,a;)e(0)(J3, (17) 

7r(t, x) = 4’o(i, x) + e(t)h{x)a 3 + ip{t, x)e(0)a3, (18) 
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where 4>o{t,x) and TTo{t,x) are free field solution to the equations of motion 
in absence of the coupling, with the same initial conditions 



0o(O,a;) = 7 t(0,x), (19) 

0o(O,a;) = </>(0,a;), (20) 

while t/’ is also a solution of the free field equations but with initial conditions 

V-(0,cr) = 0, (21) 

^(0, cc) = — /i(x). (22) 



If we examine this for the two possible eigenstates of (T 3 , we find the two 
solutions 

(j)±{t,x) Ki 4>o{t,x) ± (^e(t) J ^e"™l"=“"^'l/i(x')da:' + V'(t,a;)^ , (23) 

7T±(t, x) « (po{t, x) + 0(e) ± {e{t)h{x) + e{0)ip{t, x)). (24) 

These solutions neglect terms of higher derivatives in e. The state of the field 
is the vacuum state of (pOiT^o- 4>± and 7 t± are equal to this initial field plus 
c-number fields. Thus in terms of the (j>± and tt±, the state is a coherent state 
with non-trivial displacement from the vacuum. Writing the fields in terms 
of their creation and annihilation operators, 

+ ^ (25) 

7T±{t, x)=ij Ak±{t)e''^^ - 4±e”‘"" (26) 

we find that we can write Ak± in terms of the initial operators A^o as 

Au±{t) « Afcoe-‘“'=* ± - e{0)e-^^^*){h{k)/^k + 0{e{t))), (27) 

where = \/k'^ + w? ■ Again I will neglect the terms of order e in comparison 
with the e terms. Since the state is the vacuum state with respect to the initial 
operators A^o, it will be a coherent state with respect to the operators Ak±, 
the annihilation operators for the field at time t. We thus have two possible 
coherent states for the field, depending on whether the spin is in the upper 
or lower eigenstate of But these two coherent states will have a small 
overlap. If A\a) = a\a) then we have 

|a) = e“^^-l“l'/2|o). (28) 

Furthermore, if we have two coherent states ja) and \a'), then the overlap is 
given by 

{a\a') = (0|e“*^-l“l"/^e^^^-l^l'/2|o) = e“*/5-(l“l"+l/3l')/2. 



(29) 
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In our case, taking the two states these correspond to coherent states 
with 

« = -«'= ^i(e(t) - e(0)e-‘“'“‘) = ]^it{t)h{k) / ^. (30) 

Thus we have 






-t{ty\h(k)\^ /{k^+m^ 



— e{ty f dk /oi \ 

= e ' ^ = J(t). (31) 



Let us assume that we began with the state of the spin as :^(|+) + |— ))• 
The state of the system at time t in the Schrodinger representation is 

^(l+)l+0 W) + l-)l-0)) 

and the reduced density matrix is 

P = ^(l+)(+l + l-)H + r{t)\+){-\ + J(t)|-)(+|). ( 32 ) 



The off diagonal terms of the density matrix are suppressed by the function 
J(t). J(t) however depends only on e{t) and thus , as long as we keep i small, 
the loss of coherence represented by J can be reversed simply by decoupling 
the system from the environment slowly. 

The apparent decoherence comes about precisely because the system in 
either the two eigenstates of 173 drives the field into two different coherent 
states. For large e, these two states have small overlap. However, this distor- 
tion of the state of the field is tied to the system, tt changes only locally, and 
the changes in the field caused by the system do not radiate away. As e slowly 
changes, this bound state of the field also slowly changes in concert. However 
if one examines only the system, one sees a loss of coherence because the field 
states have only a small overlap with each other. 

The behaviour is very different if the system or the interaction changes 
rapidly. In that case the decoherence can become real. As an example, con- 
sider the above case in which e(t) suddenly is reduced to zero. In that case, 
the field is left as a free field, but a free field whose state ( the coherent state) 
depends on the state of the system. In this case the field radiates away as real 
(not bound) excitations of the scalar field. The correlations with the system 
are carried away, and even if the coupling were again turned on, the loss of 
coherence would be permanent. 



3 Oscillator 

For the harmonic oscillator coupled to a heat bath, the Hamiltonian can be 
taken as 

H J [{TT{x)-e{t)q{t)h{x))^ + {dx4'ix))‘^ + m‘^(j){t,xf]dx+^{p^ + P2‘^q‘^). 

(33) 
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Let us assume that m is much larger than 17 or the inverse timescale of change 
of e. The solution for the field is given by 

, 'YYt I 31 I 

4>{t,x) K. 4>Q{t,x) +'il}{t,x)e{Q)q{Q) - e{t)q{t) I h{x')dx', (34) 



7r(t, x) « x) + Tpit, x)e{0)q{0) 

,, p 777- 1 31 31 I 

-e{t)q{t) J h{x')dx' + e(t)q{t)h{x), (35) 

where again (j>o is the free field operator, is a free field solution with ip{0) = 
0, ^/’(O) = —h{x). Retaining terms only of the lowest order in e, 

(j>{t,x) ^ (j)o{t,x), (36) 

7r(t, x) « x) + e{t)q{t)h{x). (37) 

The equation of motion for q is 

q{t)=p{t), (38) 

p{t) = - q + e{t)<i>{t) , (39) 



where <?(t) = f h(x)c/>(t, x)dx. Substitution in the expression for (/>, we get 

p— m|3i— 31^1 



g(t) + f2^g(t) ^ e(t)^o(i) ~ e(t)e(t)(?(i) 



h{x)h{x') da;dxL (40) 



Neglecting the derivatives of e (i.e., assuming that e changes slowly even on 
the time scale of 1/17), this becomes 



p p jyq I I \ 

I + e(t)^ / h{x)h{x ')~ — ^ dxdx'\ q+ fl:^q = (41) 



The interaction with the field thus renormalizes the mass of the oscillator to 

p— m|3i— 31^1 

J\ ^ 



M=l + e 



r f p-m\x-x I 

(t)^ / / h{x)h{x') ^ dxdx'. 



The solution for q is thus 

q{t) « q(0) cos J7(t)dt^ + i sin J7(t)dt^ p(0) 

+ i J sin n{t)dt^ dt{e{t')e{t)<Po{t')dt' , (42) 

where 17(f) « Q j ^ M(£) . 

The important point is that the forcing term dependent on is a rapidly 
oscillating term of frequency at least m. Thus if we look for example at ( 9 ^), 
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the deviation from the free evolution of the oscillator (with the renormalized 
mass) is of the order of 

J sin(l7t — t') sin(a;(t — 



But is a rapidly oscillating function of frequency at least m, 

while the rest of the integrand is a slowly varying function with frequency 
much less than m. Thus this integral will be very small (at least f2/m but 
typically much smaller than this depending on the time dependence of e). 
Thus the deviation of q{t) from the free motion will in general be very very 
small, and I will neglect it. 

Let us now look at the field. The field is put into a coherent state which 
depends on the value of q, because Tr{t,x) « 4>o(t,x) + e{t)q{t)h{x). Thus, 

Ak{t) « aofee"‘“'“‘ + \hi{k)e{t)q{t)/u>k. (43) 

The overlap integral for these coherent states with various values of q is 

n(i^^(A)e(t)gM|i^^(A:)e(t)g'M) (44) 

k 

The density matrix for the Harmonic oscillator is thus 

= (45) 

where po is the density matrix for a free harmonic oscillator (with the renor- 
malized mass). 

We see a strong loss of coherence of the off diagonal terms of the density 
matrix. However this loss of coherence is false. If we take the initial state for 
example with two packets widely separated in space, these two packets will 
loose their coherence. However, as time proceeds, the natural evolution of the 
Harmonic oscillator will bring those two packets together {q — q' small across 
the wave packet). For the free evolution they would then interfere. They still 
do. The loss of coherence which was apparent when the two packets were 
widely separated disappears, and the two packets interfere just as if there 
were no coupling to the environment. The effect of the particular environment 
used is thus to renormalise the mass, and to make the density matrix appear 
to loose coherence. 

4 Spin Boson Problem 

Let us now complicate the spin problem in the first section by introducing 
into the system a free Hamiltonian for the spin as well as the coupling to 
the environment. Following the example of the spin boson problem, let me 
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introduce a free Hamiltonian for the spin of the form |l7cri, whose effect is 
to rotate the 0-3 states (or to rotate the vector p in the 2 — 3 plane) with 
frequency 17. 

The Hamiltonian now is 

H i [(7r(t, x) - e(t)h(x)a 3 y + (d^<p(x)y + x)'^]dx + , 

. . 

where again e(t) is a slowly varying function of time. We will solve this in the 
manner of the second part of section 2. 

If we let 17 be zero, then the eigenstates of cr^ are eigenstates of the 
Hamiltonian. The field Hamiltonian (for constant e) is given by 

H± = ^ y [(tt - {±e{t)h{x))f + {d^(l)f]dx. (47) 

Defining tt = tt — (±/i(x)), tt has the same commutation relations with tt 
and (j) as does tt. Thus in terms of tt we just have the Hamiltonian for the 
free scalar field. The instantaneous minimum energy state is therefore the 
ground state energy for the free scalar field for both H±. Thus the two states 
are degenerate in energy. In terms of the operators tt and </>, these ground 
states are coherent states with respect to the vacuum state of the original 
uncoupled (e = 0) free field, with the displacement of each mode given by 

afc|±) = ±ie(i)^^y=|±), (48) 

or 

i±) = n I ^ 

k 

where the ja^) are coherent states for the modes with coherence param- 
eter Qffe = ie(t)^l, and the states |±)<t 3 are the two eigenstates of 0-3. (In 
the following I will eliminate the symbol.) The energy to the next excited 
state in each case is just m, the mass of the free field. 

We now introduce the fiux as a perturbation parameter. The two lowest 
states (and in fact the excited states) are two-fold degenerate. Using degen- 
erate perturbation theory to find the new lowest energy eigenstates, we must 
calculate the overlap integral of the perturbation between the original degen- 
erate states and must then diagonalise the resultant matrix to lowest order 
in 17. The perturbation is |l7cri . All terms between the same states are 
zero, because of the (±|cT3cri|±)cr3 = 0. Thus the only terms that survive for 
determining the lowest order correction to the lowest energy eigenvalues are 

y+|l2ai|-) = y-|I7ai|+)* 

k k 



(50) 

(51) 
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= = ^nj{t). (52) 

The eigenstates of energy thus have energy of E{t)± = Eq± ^17 J(t), and the 
eigenstates are \J\{\+) ± |— ))• If e varies slowly enough, the instantaneous 
energy eigenstates will be the actual adiabatic eigenstates at all times, and 
the time evolution of the system will just be in terms of these instantaneous 
energy eigenstates. Thus the system will evolve as 

|V.(t)) = ((c+ + c_)e-'/ + |-)) 

+ (c_ - c+)e+‘/ - I-))) , (53) 



where the c_|_ and c_ are the initial amplitudes for the |+)ct 3 and |— )o-3 states. 
The reduced density matrix for the spin system in the basis can now be 
written as 

p{t) = {J{t)poi{t),J{t)po 2 {t),pos{t)), (54) 

where po{t) is the density matrix that one would obtain for a free spin half 
particle moving under the Hamiltonian J{t)Qa\, 



Poiii) — Pi(0), 

P02(t) = P2(0) cos ^17 J J(t')df'^ + P3(0) sin ^17 J , (55) 

Po3(^) = P3(0) cos ^17 J -P2(0)sin ^17 J J{t')d1^ . 



Thus, if J(t) is very small (i.e., e large), we have a renormalized frequency 
for the spin system, and the the off diagonal terms (in the representation) 
of the density matrix are strongly suppressed by a factor of J{t). Thus if we 
begin in an eigenstate of CT3 the density matrix will begin with the vector p 
as a unit vector pointing in the 3 direction. As time goes on the 3 component 
gradually decreases to zero, but the 2 component increases only to the small 
value of J{t). The system looks almost like a completely incoherent state, 
with almost the maximal entropy that the spin system could have. However, 
as we wait longer, the 3 component of the density vector reappears and grows 
back to its full unit value in the opposite direction, and the entropy drop to 
zero again. This cycle repeats itself endlessly with the entropy oscillating 
between its minimum and maximum value forever. 

The decoherence of the density matrix (the small off diagonal terms) ob- 
viously represent a false loss of coherence. It represents a strong correlation 
between the system and the environment. However the environment is bound 
to the system, and essentially forms a part of the system itself, at least as 
long as the system moves slowly. However the reduced density matrix makes 
no distinction between whether or not the correlations between the system 
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and the environment are in some sense bound to the system, or are correla- 
tions between the system and a freely propagating modes of the medium in 
which case the correlations can be extremely difficult to recover, and certainly 
cannot be recovered purely by manipulations of the system alone. 



5 Instantaneous Change 

In the above I have assumed throughout that the system moves slowly with 
respect to the excitations of the heat bath. Let us now look at what happens 
in the spin system if we rapidly change the spin of the system. In particular 
I will assume that the system is as in section 1, a spin coupled only to the 
massive heat bath via the component of the spin. Then at a time to, I 
instantly rotate the spin through some angle 6 about the 1 axis. In this case 
we will find that the environment cannot adjust rapidly enough, and at least 
a part of the loss of coherence becomes real, becomes unrecoverable purely 
through manipulations of the spin alone. 

The Hamiltonian is 

H = ^ y [(7r(t, x) - e{t)h{x)a 3 )‘^ + x))'^ + x)]dx 

+9 /25{t — to)<Ji- (56) 

Until the time to o’a is a constant of the motion, and similarly afterward. 
Before the time to, the energy eigenstates state of the system are as in the 
last section given by 

l±)^) = {l+)<T3|afe(i)) or {|-)a3l -ak{t))}. (57) 

An arbitrary state for the spin-environment system is given by 

IV’) = c+|-k) -kc_|-). (58) 

Now, at time to, the rotation carries this to 

|(/)(to)) = c+(cos(6»/2)|-h)<,3 -kisin(6»/2)|-)<,3)|afc(t)) 

-k c_(cos( 6»/2)|-),^3 -hisin(6»/2)|-h),^3)| - afc(t)) 

= cos(6l/2)(c+|-k)-kc_|-)) (59) 

+ isin(6l/2)(c+|-),^3|afc(t)) -c_|-k),^3| -ak{t)). 

The first term is still a simple sum of eigenvectors of the Hamiltonian after 
the interaction. The second term, however, is not. We thus need to follow 
the evolution of the two states |— )<t 3 |afc(to)) and |-I-)(T3| — cek{to))- Since ao 
is a constant of the motion after the interaction again, the evolution takes 
place completely in the field sector. Let us look at the first state first. (The 
evolution of the second can be derived easily from that for the first because 
of the symmetry of the problem.) 
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I will again work in the Heisenberg representation. The field obeys 



= TT_{t,x) + e{t)h{x), (60) 

fr_(t, x) = x) — a;). (61) 

At the time to the field is in the coherent state \ak)- This can be represented 
by taking the field operator to be of the form 

(j)-{to,x) = (l)o{to,x), (62) 

7T_{to,x) = ^o{to,x) + e{to)h{x), (63) 



where the state \ak) is the vacuum state for the free field 4>q. We can now 
solve the equations of motion for <j)_ and obtain (again assuming that e(t) is 
slowly varying) 



(j)-{t,x) = (j)o{t,x) + 2f/'(t,x)e(to), (64) 

7T_(t, x) = <po{t, x) + x)e{to) — e{t)h{x), (65) 

where V’(to; a;) = 0 and ipito, x) = h{x). Thus again, the field is in a coherent 
state set by both 2e(to)V' and e{t)h{x). The field ip propagates away from the 
interaction region determined by h{x), and I will assume that I am interested 
in times t a long time after the time tg- At these times I will assume that 
f h(x)'ip(t, x)dx = 0. (This overlap dies out as The calculations can 

be carried out for times nearer to as well — the expressions are just messier 
and not particularly informative.) 

Let me define the new coherent state as | — ak{t) + Pk{t)), where at is as 
before and 

(3k{t) = 2e{to)uJk-ip{t, k) = 2ie(to)e‘“'‘*/i(fc)/wfc. (66) 

(The assumption regarding the overlap of h{x) and ip{t) corresponds to the 
assumption that / a^(t)/3fc(t)dfc = 0). Thus the state |— )<T 3 |Q!fc) evolves to 
the state |— )ct 3 I — + f3k{t))- Similarly, the state |+)ct 3 | — ctk) evolves to 

|+)<T3|afe - Pk(t)).) 

We now calculate the overlaps of the various states of interest. 



{ak\ak ± Pk) = (— afcl — afc ± Pk) = e~ = J(to), (67) 

{—ak\ak ± Pk) = (afcl — afc ± Pk) = J{t)J{to), (68) 

(-Ofc + Pk\ak - Pk) = {-Oik - PkWk + Pk) = J{t)J{tQY. (69) 

The density matrix becomes 

P3 = cos(6»)po3 + sin( 6 ») J(to)po 2 , (70) 

Pi = J{t) (cos( 6 ») + J‘‘(to)sin( 6 »)) poi, (71) 

P 2 {t) = J{t) (- sin( 6 »)po 3 + (cos(6»/2) - J^itP) sin( 6 >))po 2 ) , (72) 
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where 



P03=^(|c+p-|c_n, 


(73) 


Poi = Re(c+c)l), 


(74) 


P 02 = Im(c+c)l). 


(75) 



If we now let e{t) go slowly to zero again ( to find the ‘real’ loss of coherence), 
we find that unless poi = P 02 = 0 the system has really lost coherence during 
the sudden transition. The maximum real loss of coherence occurs if the 
rotation is a spin flip (6* = tt) and po 3 was zero. In that case the density 
vector dropped to J(to)^ of its original value. If the density matrix was in an 
eigenstate of on the other hand, the density matrix remained a coherent 
density matrix, but the environment was still excited by the spin. 

We can use the models of a fast or a slow spin flip interaction to discuss the 
problem of the tunneling time. As Leggett et al. argue [3], the spin system is 
a good model for the consideration of the behaviour of a particle in two wells, 
with a tunneling barrier between the two wells. One view of the transition 
from one well to the other is that the particle sits in one well for a long 
time. Then at some random time it suddenly jumps through the barrier to 
the other side. An alternative view would be to see the particle as if it were 
a fluid, with a narrow pipe connecting it to the other well- the fluid slowly 
sloshing between the two wells. The former is supported by the fact that if 
one periodically observes which of the two wells the particle is in, one sees 
it staying in one well for a long time, and then between two observations, 
suddenly finding it in the other well. This would, if one regarded it as a 
classical particle imply that the whole tunneling must have occurred between 
the two observations. It is as if the system were in an eigenstate and at some 
random time an interaction dipped the particle from one well to the other. 
However, this is not a good picture. The environment is continually observing 
the system. It is really moved rapidly from one to the other, the environment 
would see the rapid change, and would radiate. Instead, left on its own, the 
environment in this problem ( with a mass much greater than the frequency 
of transition of the system) simply adjust continually to the changes in the 
system. The tunneling thus seems to take place continually and slowly. 

6 Discussion 

The high frequency modes of the environment lead to a loss of coherence 
(decay of the off-diagonal terms in the density matrix) of the system, but as 
long as the changes in the system are slow enough this decoherence is false- 
it does not prevent the quantum interference of the system. The reason is 
that the changes in the environment caused by these modes are essentially 
tied to the system, they are adiabatic changes to the environment which can 
easily be adiabatically reversed. Loosely one can say that coherence is lost by 
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the transfer of information (coherence) from the system to the environment. 
However in order for this information to be truly lost, it must be carried 
away by the environment, separated from the system by some mechanism 
or another so that it cannot come back into the system. In the environment 
above, this occurs when the information travels off to infinity. Thus the loss 
of coherence as represented by the reduced density matrix is in some sense 
the maximum loss of coherence of the system. Rapid changes to the system, 
or rapid decoupling of the system from the environment, will make this a 
true decoherence. However, gradual changes in the system or in the coupling 
to the external world can cause the environment to adiabatically track the 
system and restore the coherence apparently lost. 

This is of special importance to understanding the effects of the envi- 
ronmental cutoff in many environments [3]. For “ohmic” or “superohmic” 
environments (where h does not fall off for large arguments), one has to in- 
troduce a cutoff into the calculation for the reduced density matrix. This 
cutoff has always been a bit mysterious, especially as the loss of coherence 
depends sensitively on the value of this cutoff. If one imagines the environ- 
ment to include say the electromagnetic field, what is the right value for this 
cutoff? Choosing the Plank scale seems silly, but what is the proper value? 
The arguments of this paper suggest that in fact the cutoff is unnecessary 
except in renormalising the dynamics of the system. The behaviour of the 
environment at frequencies much higher than the inverse time scale of the 
system leads to a false loss of coherence, a loss of coherence which does not 
affect the actual coherence (ability to interfere with itself) of the system. 
Thus the true coherence is independent of cutoff. 

As far as the system itself is concerned, one should regard it as “dressed” 
with a polarization of the high frequency components of the environment. One 
should regard not the system itself as important for the quantum coherence, 
but a combination of variables of the system plus the environment .What is 
difficult is the question as to which degrees of freedom of the environment are 
simply dressing and which degrees of freedom can lead to loss of coherence. 
This question depends crucially on the motion and the interactions of the 
system itself. They are history dependent, not simply state dependent. This 
makes it very difficult to simply find some transformation which will express 
the system plus environment in terms of variables which are genuinely inde- 
pendent, in the sense that if the new variable loose coherence, then that loss 
is real. 

These observations emphasise the importance of not making too rapid 
conclusions from the decoherence of the system. This is especially true in 
cosmology, where high frequency modes of the cosmological system are used 
to decohere low frequency quantum modes of the universe. Those high fre- 
quency modes are likely to behave adiabatically with respect to the low fre- 
quency behaviour of the universe. Thus, although they will lead to a reduced 
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density matrix for the low frequency modes which is apparently incoherent, 
that incoherence is likely to be a false loss of coherence. 
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